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ABSTRACT 


This dissertation deals with theoretical aspects of 
the determination of the magnitude and the rate of decay 
of hydrostatic excess pore pressures caused by applied ex- 
Mernalmstcesses in diifass vor sonl. Numerical mectiods sor 
the solution of the problems are developed. 

The existing theories of primary consolidation used 
to obtain the compression characteristics of clay are criti- 
cally reviewed. The techniques of solution are mentioned 
briefly and the types of boundary value problems that may be 
encountered are summarized. 

A general theory of one-dimensional primary consolida- 
tion accounting for the self weight of the soil grains is 
presented. This general theory is based on nonlinear re- 
lationships between the coefficient of permeability, the 
coefficient of compressibility and the effective stress. 
Terzaghi's linear theory of consolidation can be considered 
a special=ease of this general theory. The general theory 
evolved can be applied to the case of a sedimenting soil. 
The data recorded from the Mississippi Continental Shelf 
is recoverable with appropriate input data. 

This dissertation further develops a numerical method 
in two dimensions for determining the construction pore 
pressures in embankment sections as a function of their 
construction history. The numerical algorithm developed 


is versatile. Problems involving complex boundary conditions, 
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stratified soils and gradual loading can be solved. The 
coefficient of consolidation can also be varied as a function 
of position and time. Case records of six embankments have 
been analyzed and results compared with field observation. 
Generally good agreements were obtained. 

The numerical algorithm developed is extended to treat 
the cases of impeded drainage (retarded consolidation) at 
the junctions of inner cores and outer shells of earth em- 
bankments. The impedance to drainage is measured by a 
parameter called the Impedance Factor. The impedance fac- 
tor depends on the relative lengths of the drainage paths 
and the relative permeabilities of the inner cores and the 
outer shells of the embankment. Results are presented for 
a hypothetical dam wnose dimensions are varied and for 
different values of the impedance factor. 

Terzaghi's classical theory of one-dimensional con- 
solidation is extended to develop a second-order weakly non- 
linear partial differential equation governing the pore 
pressure equalization arising out of the erosion of a fully 
saturated sediment at a prescribed rate. This part of the 
study aims at correlating some of the geological and physical 


parameters involved in valley formation. 
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CHAPTER I 


INTRODUCTION 


[Nel rURPOSE Ore THERRESEARCH 


The study comprising this thesis is based on the 
well-known classical Terzaghi theory. In deriving the 
differential equation governing one-dimensional consolida- 
tion, Terzaghi assumed constant coefficients of permeabi- 
lity and compressibility for the soil together with 
Paaimrclesinalestrains ito realizesa hinearntheory. veSuch 
assumptions characterized early studies in mathematical 
physics for two very important reasons. First, a linear 
theory yields governing equations which are amenable to 
Closed form solution. This requirement was almost essential 
prior to the advent of high speed computers. Second, a 
linear theory enables one to apply the principle of super- 
position. However, the idealizations assumed by Terzaghi 
obviously are at variance with reality. 

Differences in observed and predicted results have 
stimulated a series of modifications to the theory. Non- 
homogeneous soils whose properties vary with depth have 
been treated by Schiffman and Gibson (1964) and Raymond 
(1965). Small strain theories for compressible, normally 
consolidated clay incorporating stress-dependent coefficients 
of compressibility and permeability have been presented by 


Davis and Raymond (1965) and Barden and Berry (1965). Large 
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Strains have been considered by Mikasa (1965) and by Gibson 
eum be (1967) " Gibson et lal. “Cop? cit!) “have ‘incorporated 
the relative velocity between water and soil in applying 
Darcy's law. Many of the above treatments lead to non- 
linear partial differential equations which cannot be 
solved rigorously. Recourse is usually made to numerical 
methods. 

Under three-dimensional conditions, three-dimensional 
strains as well as three-dimensional flow of water take 
place. The three-dimensional treatments have all assumed 
Hinear elastic characteristics, as anything more realistic 
leads to extreme mathematical complexity. Pseudo three- 
dimensional cases which involve three-dimensional flow of 
water but only one-dimensional vertical strain are compara- 
tively simple to treat if the bulk stress remains constant. 
The solutions for sand drain design (Barron, 1948) are of 
this type, as also is the more complex analysis for sand 
drains in layered clays (Horne, 1964). 


All of the studies mentioned above neglect the true 


self-weight of the soil. 


Ve cee COREOPs THE wo TUDY 


The scope of this study includes a summary of all of 
the existing theories of one- and three-dimensional consoli- 
dation. The assumptions involved in each theory are 


examined and criticized. A general, nonlinear theory of 
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one-dimensional consolidation is developed incorporating 
muerse lf welgnt of the soil, the relative velocity Ons Oule 
anadawater in Darcy's law, and a variation of the coeffi- 
cients of permeability and compressibility with effective 
eeress Which in turn is a function of depth and time. 

A two-dimensional plane strain case for consolidation 
is developed on the basis of Terzaghi's assumptions, which 
involves two-dimensional dissipation of water and one- 
dimensional vertical compression. The development of the 
equation is general and incorporates a moving boundary 
involving the thickness of the clay layer increasing with 
mame. Ihis equation is used to predict the construction 
pore pressures and the effect of dissipation during con- 
Struction shutdown. Six case histories have been studied 
and the calculated and recorded pore pressures are compared. 

Arising out of the solution of the two-dimensional 
equation, rectangular clay cores having side drains are 
examined. The impedance of the side drains to the normal 
Tiow Of Water due to dissipation is studied. 

Woaparteote ant investigation antoephocesses associated 
with valley formation the erosion of a fully saturated soil 
mass at a specified rate has been studied and the resulting 
pore pressure isochrones have been computed. This is a 
corollary to the one-dimensional sedimentation case (Gibson, 
1958). The equation governing one-dimensional erosion is 


developed and results are obtained for typical cases. 
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TSS8 ORGANEZAPIONG@ES THE THESIS 


This study can be divided into the following main 
Dain S-. | 

Chapter Lil deals with a brief critical review of 
me existing theories of consolidation of clay layers. 

Chapter III comprises a brief treatment of the 
boundary conditions that may arise in consolidation problems. 

Chapter IV includes a summary of the numerical tech- 
niques involved and a detailed discussion of the methods 
employed. 

Chapter V is a discussion of the Terzaghi two-dimen- 
Sional consolidation case. Also an extension of the Davis- 
Raymond theory to two dimensions is included. A complete 
section is devoted to the development and discussion of 
impeded drainage in two-dimensional dissipation. 

Chapter VI deals with the analysis and techniques 
of solution of a two-dimensional consolidation equation 
which yields the construction pore pressures in a growing 
wedge. A summary of the theoretical methods for predicting 
construction pore pressures is provided. Six case histories 
are presented and comparisons are made between the numeri- 
cal and recorded pore pressure results. 

Chapter VII is the development of analysis and 
solution of a general, nonlinear, finite strain, one-dimen- 
sional consolidation theory for clay layers increasing in 


thickness with time. The development accounts for the 
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Selr weight of the soil grains?*? The material property 
behavior is based on observed data. A case history is 
cited. 

Chapter VIII refers to the analysis of the one- 
dimensional equation governing the deficient pore pressure 
dissipation arising out>of the depositionalvcand erosional 
balance in valley~-formation. Results are provided for few 


general cases. 
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CRAPTER A 


REVIEW OF CONSOLIDATION THEORIES 


2elerGENERAL 


When a load is applied to a saturated soil mass, 
initially the entire load is taken up by the water in the 
pores of the soil skeleton. fhe pore water thus develops 
a pressure excess (in excess of the equilibrium pore water 
pressure) equal to the applied load at the instant of load- 
mig, If the duration of load application to the saturated 
Somemeiiess 1S sShOrt,.;1.e,5 Short» compared to. the dissipation 
time of excess pore pressure, water begins to flow due to 
the hydraulic gradient caused by the excess pore pressure, 
and the soil changes in volume. The process involving a 
decrease of the water content of a saturated soil without 
replacement by air is called the ‘Process of Consolidation’. 
The mathematical theory describing the dissipation of excess 
pore pressures and associated deformation of the soil is 
called 'Consolidation Theory'. The amount of compression 
thacenassoccurrced in the Soil ‘skeleton at any time is ‘rela- 
ted not only to the applied load but also to the amount of 
Sumess transmitted at tne usOll, particie contacts. .1.e., to 
the difference between the applied stress and the pore pres- 


Sure. this differences called stha;'Effective Stress’. 
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2.2 BIOT'S THEORY OF THREE-DIMENSIONAL CONSOLIDATION 


Biot (1941) has presented a rigorous and complete 
Preaument of tne theory of, consolidation. The following 
basic properties of the saturated soil were assumed: 

(a) isotropy of the material 

(b) reversibility of stress-strain relations under 

final equilibrium conditions 

(c) linearity of stress-strain relations 

(d)® small strains 

(e) the water contained in the pores is incompressible 
and (f) the water flows through the porous soil skeleton 

PCCOrdiINnGs LOMVareyns Faw: 
The strains are assumed to be related to the effective 


Stress changes by equations such as 


Xian : 5 a i 2 eta 
where Evy denotes strain in the x-direction 
Cae denotes effective stress in the x-direction 
. (ory - u) 
oy denotes total stress in the x-direction 
U denotes pore water pressure 


u denotes Poisson's ratio of the soil skeleton 
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were > is called the first effective stress imvariant, and 
eeenemtTiysit total stress. invariant. Similar’ relations: such 
as 2.la can be written for strains in the y- and z-direc- 
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Equation 2.4 yields 


EA 


i Sees): 


substituting 2.]lc in the above equation leads to 


A EA 
6 = 3u re Enea 
ie) 
iL bea 
or 3 0 = uU'= Schmealt): 


purcher, Darcy *s law’ is assumed” valid for ‘the ~flow 


of pore water through the porous soil skeleton. Therefore, 


; k du 
Vo = =~ kdl, =-= 206 
» % Gy aX 
where aS denotes the velocity of fluid flow in the 
x-direction 
i denotes the hydraulic gradient in the x-direction 


k denotes the coefficient of permeability of the 
S-O: lake 


On differentiation, equation 2.6 becomes 
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Equation 2.10 is the general equation for consolidation 
in three-dimensions. An unique solution to equation 2.10 


depends on the specified boundary and initial conditions. 


Boundary and initial conditions are discussed in a subseque 


section. 

Of the basic assumptions made, (b) and (c) are most 
subject to criticism. Biot (1941) contends that it can be 
imagined that the grains comprising the soil are held to- 
gether by surface tension forces and tend to assume a con- 
figuration of minimum potential energy; this Biot believes 
is essentially true for the colloidal particles constitu- 


ting clay. Biot (op. cit.) assumed that for small strains, 
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when the soil grain pattern is not too much disturbed, the 
assumption of reversibility will be applicable. The assump- 
tion (a) of isotropy is not essential and anisotropy can 
easily be introduced as a refinement. 

In Biot's presentation it is emphasized that all 
equilibrium, compatibiiity and stress-strain relationships 
mwvempeen satisfied. Lt is; evident that, the formal solution 
of the equation 2.10a is rather complex; only few special 
cases have been solved. For example, De Jong (1957) pre- 
senved the solution for a circular uniformly Joaded area 
Supported on a compressible soil of infinite extent and 
Gibson and McNamee (1957) solved the case of a rectangular 
uniformly loaded area on the surface of a semi infinite, 
POuropic, homogeneous, and fully saturated porous. elastic 
medium with a fully permeable upper surface’ he analysis 
was restricted to a medium with a Poisson's Ratio equal to 


ZEW . 


2.3 RENDULIC'S THEORY OF THREE-DIMENSIONAL CONSOLIDATIO 


Rendulic (1936) presented a theory for three-dimen- 
ional consolidations, Assuming that, the soil) as isotropic 
and homogeneous, the three-dimensional consolidation equa- 


tion may be written as 
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where a, =. Ove ay tee eta eee © 12 @) Ee) 


and Cus __kE 
Tenens 


This equation was obtained by assuming a linear elastic 
stress-strain law, and no attempt was made to refer to 
strain compatibility. This theory is called the pseudo- 
three-dimensional theory of consolidation. 

GUationg2s Li isesvdenticalewi tiie. 10, far an 
castre skeleton; but it is “important to, realise that 8, 
in 2.1] is the sum of the externally applied total stresses, 
whereas 6 in 2.10a is the sum of the total stresses at the 


point wonder consideration. The value of oG is determined 


ot 
not only by the time rate of change of external stresses 
Che) 
(as is srl)» but also by the change in volume as the con- 


solidation process progresses (Schiffman, 1967). 


2.4 (COMPARISON .OF “FHEORTES*OR=bBi OT AND RENDULIC 


re) : : 
Generally the term y_ 1S Non-zero, since stress re- 


distribution does take place within the soil mass during 
consolidation, even when the applied load is maintained 

constant. The term oul will be zero when the externally 
applied load is constant. When, as a first approximation 


“¢ is taken equal to zero (as in Terzaghi theory, to be 


followed), then the governing equation for consolidation is 
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As pointed out earlier, under three~dimensional 
Seredan ‘conditions, 2.12 assumes the form 
2 2 ks 
o,(25 + Sh + 28) = ou Zi 
0X Oy dZ . 
where One ae a 
ii 3 Sil = 2u)y 
Under two-dimensional Strain conditions, the strain 
in one of the directions (say y-direction) is considered to 
be zero (plane strain case) and it can be shown that equa- 
tion 2.12 assumes the form 
2 2 
c, (2-4 ot) = + 2.14 
OX 0Z 
it kE 
where Copal eoeeih (le eal Te 
and is the coefficient of consolidation under two-dimen- 
Sonal strains 
For one-dimensional strain condition, the equation 
becomes 
2 
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and is the coefficient of consolidation under one-dimen- 
2 EL Strat. 

Thus it may be seen that the form of the equation 
meicsdenends: on the tTlow (or drainage)” conditions, while 
the value of the coefficient of consolidation depends on 


erems crain COndTcions.9 Ltewild beemoted that 
: ; Src W llamas 
Cy 2(1 U)C, Si )c, 


Oma Valueof q =" 0.5 Cy = Coes C. 
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POrua Value of 1 Cy Co 3 
The difference between these values of c can be important. 
inner value /of.c in the-definition of time factor Torta tri- 
exdaletest 1S C25 whereas that for the oedometer test is Cc). 
The three-dimensional theory of primary consolida- 
tion (poroelasticity) due to Biot makes no assumptions as 
to the time-dependent variables. The state of total stress 
is a function of the excess pore pressure and the compati- 
Dit yerelationsnwps TOM wmte son skeleton. Ibis. theory 
is completely self-conststent and does not call for a sepa- 
ration of the magnitude and progress of consolidation. The 
pseudo-three-dimensional theory maintains a lack of coupling 


between magnitude and progress of consolidation. 
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Comparing these two theories, Schiffman (1967) considered 
Piemprovlem of avsitrip ‘oad (plane! strain) on the surface 
Oe Seni int initerbody.. It was) found that the total stress 
is not constant with time. Furthermore, excessive pore pres- 
sures, instead of decreasing with time, increase at first, 
and then decrease. This effect, called te Mandel-Cryer 
effect, appears only in the three-dimensional theory, and 
is absent from the pseudo-three-dimensional theory. 

The solution of three-dimensional problems and sub- 
sequent application to practical cases is seriously limited 
because of the lack of techniques of solution and also be- 


Cause accurate values of E and » ane difficult to obtain. 


2eoue TERZAGH LS CLASSICAL ONE =D MENS LOMA Wai G CORY. 


The one-dimensional consolidation theory (Terzaghi, 
1923) has been applied widely to many field problems. The 
results obtained are usually acceptable and, Wain Hacc 17 
cal limits, quite reliable. The assumptions involved in 
ine aerzaghi theory ares 

(i) homogeneous soil 
(ii) incompressible soil grains and the voids fully 
saturated by an incompressible fluid 
(aii) small strains 
(iv) void ratio is linearly related to the effective 


stress 


' | b ete . x ¢ bal Ga ¥¥ . 
By pees (foe Ma Re Mt alee tig 2 ER “piney WHF ma i 
: 7 ow Cal’? sa. nt ‘ , BY 


hy “bait PET 
“lS tied eae : 
wide hy... en, Hi Bau saison 
eet ms 


4 , d ‘ 
ey ee ; et ete | 
To) ste 7, : ap 
| | : Me 
‘ erik alice be pear te ee ih: 
bn wot Ieee gts, 
(Pi) ig 


» 


bed ag on veg 


5 ihe ey? r 
4 | “sol ara Lb er ; 
. 
ton We ie A Lee 29m 
y vi } 


vV qi 7 


j ‘ , a 


ys} PUP ead Nhe IRE 2, a 
. - a if : 


: > a? 
4 AR Si viigaelé gotbett joecen hehe Feeenen one aT cP: 
{ ’ he « a ey i (et 


we} oa idee ae ied nin. os ile hia? ‘Sar fawe ve al | 


} 


veto hil thw aed gros, Uhsiew.aa er ait 


Serna CAR? Veqtrt 3 Gsire Aili sfunt lan ¢ wie 9 abit | 


on . ms pes + dare, Gi GAs,’ ian 
oe Ser ee Schoentteg atom i) a0 


ro 


a 


VE COTRBRVOL SHochay Batero) bos 9 dees pee ARy | ae a 


. phat | saver tdigsn' ie Ye besntesea a E 


y . | nian cs id tare waa 


cy ieetvaan we hoteles” yl soon F..2 ofp wr, hte | “tat 2 


(vy) ‘the’ fluid flow through the soil skeleton is 
governed by Darcy's law. 
and (vi) permeability and compressibility are constant 
throughout the process of consolidation. 
This theory considers only the diffusion of water 
through a porous medium. For an isotropic soil in which 
water flowsonly in one direction and compression occurs only 


in the same direction, the governing equation may be written as. 


The one-dimensionality of consolidation may be due 
either to a surface loading of uniform magnitude infinitely 
extended in all directions or due to a lateral physical 
Restraint (fas in c&he oedometer test 

The application of one-dimensional theory, which 
depends on unreal assumptions as to the direction of con- 
pression and iaterial properties, paradoxically, yields 
acceptable results and permits a large variety of solutions 
to be obtained. Recent developments in numerical analysis 
and in computer technology provide improved capability in 
handling one-dimensional consolidation equation for an 
arbitrary loading history and for variations in soil pro- 
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Pb gov pl CALLONS (.0y TERZAGHIAS THEORY 


From time to time, many of the assumptions made 
in Terzagni's theory have been modified. 

It is known that the Terzaghi theory does not hold 
good for the case of soft clay whose permeability k and 
compressibility my, change their values during the consolida- 
tion process. Small strain theories for compressible 
normally consolidated clay incorporating a decreasing my, 
and k have been presented by Davis and Raymond (1965). 
Non-homogeneous soils whose properties vary spatially and 
with time have been treated by Schiffman and Gibson (1964) 
and Raymond (1965). Layered soils can also be considered. 

The assumption of the walidity of Darcy's law: 1s 
generally acceptable but a modification may be introduced. 
The seepage velocity of flow is replaced by the relative 
Verocuty, On flow oF fluldowith wespects vos tat of the so lid 
Guomnsw(schneldedger, 1960)"  sinis extends Dancy 4s lawe toa 
more presentable and accurate form. 

Tie) imitation on SWall os trains, was specen, Overcome Dy 
making use of Lagrangian coordinates (Gibson et al., 1967). 
Also Mikasa (1965) considers large strains. 

Many of the above modifications in combination or 
singly lead invariably to nonlinear partial differential 
equations governing the process of consolidation, for which 
analytical (closed form) solutions are difficult to obtain. 


In actual practice, the total load on a clay layer 
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is usually applied over a period of time, th. i ties iS 
SeneunatvOnly a Smali, proportion of the total consolidation 
occurs during this period then tne load may be considered to 
be applied instantaneously. On the other hand, if a suffi- 
ciently large proportion of the excess pore pressure has 
been dissipated during the period te then the problem may 
be treated in either of the following ways: 

(i) The time-dependency of loading may be considered 
Fomoe. tne Variation in loading With time, ice., the rate of 
mtcmaoohication 1S ‘arbitrary, @.9.-, COMms tructiIOn rate of 
loading. 

(ii) The total load on the element of soil under 
consideration may also change because of the addition of 
Supsequent layers of So1l On top of ‘the element. “Such a case 
is also one of time-dependent loading except that this case 
involves the physical movement of the top boundary. Tnis 


is a moving boundary problem. 
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BOUNDARY CONDITIONS - A BRIEF TREATMENT 


821 GENERAL 


The basic governing equation for the dissipation 
Sreexcess pore pressure than ideal soil according to 


Merzagni' S theory 1s 
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The solution of this differential equation for various 
geometrical configurations and boundary conditions* has 
been attempted by many since the advent of Fourier mathe- 
focics. “Equation 3.1 will have numerous solutions unless a 
SounOte boundary CONdTtLONS andwans initial condition: are pres 
Scribed. The conditions prescribed at the boundary surfaces 
of the region may be linear or non-linear. 

The general problem may be formulated as follows: 

It is required to determine the distribution of excess 
pore pressure in a homogeneous and isotropic (soil) body 


——s 


* The term ‘boundary condition', sometimes, may include the 
womcalled “initial concicions ... This means, according: co 
Some. schools of thought, that. “bounds are considered to 
be limitations in time as well as in space. 
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at any given instant when the following is given. 

(a) The pore pressure excess distribution at any 
Ser einstance, mostly ‘at van earlier instance (e.g., at 
eiicewe —- 05 1.e., initial pore pressure distribution). 

(b) The influence of the surroundings of the body 
Omeresesurtace. For -example, it’ may be assumed that by 
some means a definite excess pore pressure distribution can 
be imposed on the surface of the body by external forces; 
this distribution may be constant or time-dependent. 

The conditions under (a) and (b) above are called 
the boundary conditions, the first being with respect to 


time and the second with space. 


3.2 BOUNDARY VALUE PROBLEM 


In the equation 


Ul wat (Ocuyoaeze te) Geez 


u may be interpreted as given by a function f within a domain 
of a four dimensional space - time region and to solve for 

u uniquely mathematics demand certain values (called condi- 
tions) at the boundary of this domain. In mathematical 
terminology, these conditions which must be satisfied in 
addition to the governing differential equation, are known 
as boundary conditions. And the problem(s) thus posed are 


called BOUNDARY VALUE PROBLEM(S). 
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A. BOUNDARY CONDITION FOR TIME 


The boundary condition for time may be given by a 


function 


Uae) ase 


which represents the distribution of excess pore pressure 
Sony Given instant. “this distribution can be chosen, 
mebmerariby,, to be continuous or discontinuous. In most 
problems, excess pore pressure at a later instant is of 
interest. This is obtained using the governing equation in 


conjunction with the boundary conditions for space. 


B. BOUNDARY CONDITIONS IN SPACE 


In the study of boundary value problems of heat con- 
duction in mathematical physics, the conditions at the 
boundary surfaces of the region may be linear or nonlinear. 
The same holds good for excess pore pressure dissipation in 
a (soil) body as the similarity between pore pressure dissi- 
pation and temperature transmission has been well established 
(Terzaghi, 1923). In the following pages, only linear boun- 
dary conditions are considered. 

For convenience the linear boundary conditions may 


be presented as follows: 
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(i) BOUNDARY CONDITION OF THE FIRST KIND 


Excess pore pressure is prescribed at every point 
on the boundary surface and for the general case it is 


merunction Of both time and, position, i.e. ; 
u = fF, (xsy5z,t) oad 


the function itself may be arbitrary, continuous or discon- 
tinuous with respect to time and position. 

Special cases include excess pore pressure at the 
boundary surface as a function of position only, or a func- 
tion of time only, or a constant. If the excess pore pres- 


Sure u at the boundary surface vanishes, then 


This special case is called the homogeneous boundary condi- 
Pionvat the first kind.» A boundary surface which 1s (kept 
@cea constant pore pressure U.,; also satisfies the homogeneous 
Boundary condition of the Tirst Kind If the excess pore 


pressure is measured in excess of U. 


(ii) BOUNDARY CONDITION OF THE SECOND KIND 


The second kind of boundary condition consists in 
expressing the amount of water exchanged (in case of thermo- 


dynamics, hte iseutnenuguantity of heat exchange at the surface, 
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Svenieass ab runctionaos time! and» posi tion.© «This! function 
may be quite arbitrary. 

The normal derivative of excess pore pressure is 
prescribed at the boundary surface, e.g., in an orthogonal 
coordinate system it is expressed in the form 


ao 
Peper Fo(Xs¥5Z,t) 3.5a 


where d/dn denotes the differentiation along the outward 
drawn normal at the boundary surface. This boundary condi- 
tion is equivalent to that of prescribing the magnitude of 
the quantity of water exchanged along the boundary surface. 
Special cases of the above equation include the 

normal derivative of excess pore pressure at the boundary 
surtaces toubesa funetion oF POSTTiOnionlyj;oor’a function of 
time only, or a constant. If the normal derivative of excess 


pore pressure at the boundary surface vanishes, then 


This special case is called the homogeneous boundary condi- 


tion of the second kind. An impervious boundary (an insu- 


OS 


lated boundary, in thermodynamics) satisfies this condition. 
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(iii) BOUNDARY CONDITION OF THE THIRD KIND 


The third kind of boundary condition gives a linear 
combination of the excess pore pressure and its normal 
derivative at the boundary Sia cel OvnenewOnas 7. tNuS 
type consists in defining the excess pore pressure at the 
Surface together with a law for the water exchange between 
the surface of the (soil) body and its surroundings. 


In mathematical terms this may be expressed as, 


ou + = oF 
A, ah Anu f3(x,ysZ,t) 3. 6a 


meus this type tS "a combination of ‘the firstotwo types tof 
boundary conditions and each of the first two voundary condi- 
tions can be obtained by appropriately choosing Ay or Ay 

to vanish. 


A’specital’ case’ of "the above equation is, 


which is called the homogeneous boundary condition of the 
third kind. Such a situation can be realized in reality 
through the dissipation of water into a surrounding which 
impedes drainage for example, the core of a dam dissipating 
excess pore pressure into a side drain whose efficiency is 


not unity 
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A graphical representation may be useful to illust- 
rate the nature of the three types of boundary conditions 
miececussed above. -In Figure: 3,1.ds represénts an element 
on the surface of a (soil) body and n the outward drawn 
(positive direction) normal. 

In the first type of boundary value problem the 
excess pore pressure Uy at the surface is known. The hydrau- 
lic gradient and the amount of water exchanged at the boun- 
dary surface are the unknowns to be determined (Figure 3.la). 

In the second type of boundary condition problem, 
the conditions are opposite to those of the first type. 

In this type, the value of the excess pore pressure at the 
surface is the unknown, while the hydraulic gradient and/or 
the amount of water exchanged at the boundary surface is 
known (Figure 3.1b). 

In the third type of boundary value problem a point 
on tne outward drawn normal is given through which all the 
possible tangents to. thevexcess pore pressurescurvevatythe 
Simface mustipesss:. Thisepointi known pase POLES) |i essatra 
Gas tercenk/ae=chi from thetsurfaces«' (k issthe coefficient 
of permeability and a is the ‘water exchange coefficient’). 
The unknowns in this type are the excess pore pressure at 
Surface and the quantity of water exchanged. (the proof 
is given by Grober et al., 1963). 

The three types of boundary conditions described 
above cover most cases of practical interest. Boundary 


conditions associated with powers of excess pore pressure 
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other than unity are nonlinear. Also boundary conditions 
associated with changes in material content of the (soil) 
body (e.g., erosion or deposition) are nonlinear boundary 
conditions in that the boundary conditions are moving with 


respect to time. 


3.3 MOVING BOUNDARY CONDITIONS 


As mentioned in the previous paragraph, an important 
class of problems in pore pressure dissipation (as in heat 
conduction) deals with the determination of pore pressure 
(temperature) in a body whose boundaries are not fixed in 
Space. Two types of problems may possibly be distinguished: 

| (a) those’ in which the motion of the’boutndary is 
due to rigid body motion of the entire body, and 

(b) those in which the motion of the boundary is 
due to local conditions near the boundary. 

Case (a) is usually referred to as moving body pro- 
blem(s). Solutions for several problems of this type for 
bodies moving with a constant velocity are given by Carslaw 
and Jaeger (1959). 

Case (b) may be divided into two cases, namely those 
in which the motion of the boundary is prescribed, and those 
in which it must be determined as part of the solution of 
the problem. The first of these cases is simpler than the 
second, since it is a linear problem although with variable 


coefficients. An example of this is the problem solved by 
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Gibson (1958). The latter case is sometimes referred to as 
‘floating - boundary' problem. 

However, very few exact (closed form) solutions are 
monleblertor such proplems, and these problems indeeed 
tend to become rather eeoree except in few ycases. As 
a consequence approximate and/or numerical methods of solu- 
tion are very important. 

The importance of floating-boundary problems in pore 
Pressure dissipation arises primarily from the ‘problems of 
Geposttion “and erosion, iniwhich the position of the inter- 
face between the deposit and the freshly deposited material 
is not known beforehand and must be determined in the solu- 
Ponce lneds typical problem, Ghesconditivons coebessatis sped 
on the moving front are that the pore pressure should be 
the same at the interface in both media and the continuity 
of flow must be maintained. 

To simplify the analysis of problems such as’ erosion, 
it is assumed that the eroded material is completely removed 
from the surface as soon as it is eroded. © There will there- 
fore be no eroded phase on the surface of erosion: “In 
reality a certain thickness of eroded arena will always 
pemin veontact With the surface of ‘erusdon.. Sucnecases! “are 
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CHAPTER AY 


TECHNIQUES OF SOLUTION - A BRIEF SUMMARY 


4.1 GENERAL 


Wich is; chapter an examination OF TUNE erinitve dit rer. 
ence approximation of the boundary value problems of excess 
pore pressure dissipation (similar to problems of heat con- 
duction) is presented. Digital and analog computers are of 
great value for solving problems that cannot ordinarily be 
handled analytically because of complicated geometry and/or 
boundary conditions, or because the numerical evaluation of 
the analytical solution becomes too laborious. Digital com- 
puters are frequently put to use because of the fact that 
high-speed precise and versatile computers are available. 

By analogy with heat conduction theory it is possible 
eomdistinguishn threeranalytical methods for ther solubion of 


the equation (Abbott, 1960), 


Boog oallnagou 4.1 


These are: 
(a) Classical methods based on Fourier mathematics, 
(b) Green's function, 
(c) Integral transformation, such as Laplace's. 


The classical methods are useful for calculations of 
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pore pressure distribution in homogeneous soils and implicit 
in these is the underlying assumption of linearity, whereby 
the principle of superposition may be empioyed. However, 
most soils exhibit a nonlinear behavior, hence the principle 
of Superposition is no longer applicable. There are thus 
considerable objections to further employment of classical 
methods of Fourier mathematics. The same objections Jimit 
the application of Green's function. 

Integral transformations offer excellent means for 
the solution of the above equation except that computation 
OF a Set of results presents considerable algebraic and 
arithmetical problems. Generally, integral transformations 
are limited to solutions of tinear differential equations. 

Thus closed-form solutions of the above equation for 
consolidation problems of non-homogeneous and nonlinear soils 
are difficult to obtain and alternate methods must often be 
employed. 

An alternative to analytical methods is numerical 
analysis, wherein the operation of solving the equation is 


reduced to solving a series of algebraic equations. 


4.2 FUNDAMENTAL CONCEPTS IN FINITE DEEPER ENCE 
APPROXIMATIONS 


A fundamental concept in the finite difference appro- 


ximation of a differential equation is the expansion of the 
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function by Taylor's series. Once the finite difference 
approximation is obtained, the problem reduces to the solu- 
tion of a set of algebraic equations. These equations may 
be translated into a set of commands understood and per- 
formed by a digital computer. 

immune tN rOCesSS. Of =SOlLUtITON 101 diferencia lsequations:. 
eymors are introdireed at each step of the calculation due 
to approximations involved in finite difference and numeri- 
Pelecalcilations es) Cumulative effects of such errors on the 
final result and stability of the solution are very import- 
ant. An exhaustive treatment of the above aspects can be 


found in Fox (1962), Richtmeyer (1957), and others. 


423) FORMS OF FINITE DIFFERENCE APPROXIMATION 


Various schemes are available to express the differen- 
tial equation in a finite difference form. Richtmeyer (1957) 
lists thirteen different schemes, ranging from the explicit 
fouinecOrascwltiy Impiielt, form for finite di irerencingsor tne 
one-dimensional time-dependent heat conduction equation. 

An explicit scheme provides a noniterative ‘marching' 
process for obtaining the solution at each present nodal 
point in terms of known preceding and boundary points. 
Muestians OF StabilitVyeareamost Critical for explicit «schemes - 
Implicit procedures are usually iterative simultaneous cal- 
culations of many known present values together with pre- 


ceding values and boundary conditions. 
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Each of these difference schemes has 


himitations as explained in 4.5. 


3) 


ints merits: and 


When the boundary condi- 


tions involve derivatives, a forward, backward, or central 


differencing scheme may be used in expressing the boundary 


eondici10n in finite difference 


4.4 


The differential equation 


form. 


GENERAL FORM OF FINITE DIFFERENCE SCHEME 


iene nequon 0 < x < Il for 0 < fe< 2 can be putin finite 


difference form by expanding the derivatives. 


derivatives may be expanded in terms of first order central 


The space 


differences, and the time derivative as a first order forward 


Giifernence (Figure 4.1). 
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The factor A may take an arbitrary value between zero and 
Unity. 
Hires vallvevor siUts rican be found in terms of ui.) 


maeuisktl), byteulinear interpolation over 1. as 
u(i,kta) = au(i,kt1) + (1-a) u(i,k) 4.3d 


The boundary and initial conditions can also be ex- 
pressed in finite difference notation. 
From equations 4.3, equation 4.2 may be obtained in 


a general form as 


Moauci-1,k+1) - Bu(i,k+1)} 4.4 


(4-1) {A,u(4-1,k) Te CIA) ste ale 


u 


where A, B, Als By» Ci> and C, are quantities involving Az 
OnGmNteand as such are constants. 


Choosing a value for A will yield a particular method 


of computation. 
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45 os EXPLICIT, AND ‘IMPLICIT SCHEMES 


Pe EAL COI SiGH EME 


Let X» be chosen equal to zero. 


Equation 4.4 becomes 


Ci sktl) = C,u(i-1,k) it C,u(i,k) ue Cu(itl,k) ED (Ginak eee 


where Cy» Cys C., and D are constants involving Ax and AT. 
ius ethis. scheme calculates u at: point 1. in space 

emdekt!) inetime, in terms of the values of u at points 

Git). i. (i+]) at the preceding time step k. | Hence this 


method ‘marches forward' in time. 


Ba IMPIICIT SCHEME 


betes bDeeequal- Tosuuat ye 


Equation 4.4 yields 


Cju(i-1,k+1) : Cou(i,k+1) + Cau(it],k+T) 


+20 Vik ine uiek) 4.6 


This scheme involves the unknown values of u at time 
Moe iveror points (i-l),.1, (itl) to calculate the value of 


u at time k for space point i. In other words, the value of 
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u at any time is implicitly expressed in terms of unknown 
values of u at a following time step hence the name 
maplacdt. 

lines terative solue@ion fom Dhis scheme may besobtaimed 
by solving a set of (ntl) linear simultaneous equations for 


meng eunknowns u(i,k#l) for 1=0,192,...<n2 


C2 CRANK = NICHOLSON SCHEME 


Let A» be set equal to one-half. 


Another implicit scheme results and is given by 


By u(i-1,k+1) + Bou(i,k+1) + Bau(itl,k+1) + BACi,k+1) 


= A,u(i-1,k) + A,u(i,k) a A,u(itl,k) + Aa (aieaks) 


where BL B. > B. By» Al > Ay» As Ay LAVO 1 Ve™ AXe-ond elk and 
hence constants. | 

The right-hand side of equation 4.7 contains only quan- 
tities at cime K and: hence 1tais Simidian conequation 4.6. 


Selutroneor equation 4.7 is similar ~to that of equation 4.6: 


D. COMPARISON OF METHODS 


An explicit method of finite difference representation 
of the differential equation 4.2 results in a set of simple 


algebraic relations which could easily be solved with a digi- 
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tal computer. However, stability considerations restrict 
mre sizeof the time step AT for a given value of Ax. If 
the space step Ax is to be chosen small for improved accur- 
acy, then AT correspondingly also has to be chosen small. 
mas results in a large number of time steps and computation 
Oeecones unwieldy. In such cases an implicit method is 
usually employed. 

iMemaimnplicre methodeddes not tresitrich tnew value 
Giant and has no’ restriction on iteration rate. The time 
taken for each iteration may be high, but total time taken 
poreant! vverations is small; smaller than the time saeen Ber 
the explicit method. However, it is more difficult to set 
Ppeene calculation procedure’ for an implicit method= than for 


the explicit method. 


4.6 ALTERNATING - DIRECTION IMPLICIT METHOD 


As noted in a previous paragraph, the implicit method 
has the advantage that it is unconditionally stable for all 
values of time step AT. On the other hand, the computational 
problems become enormous when two- or three-dimensional 
pore pressure dissipation (heat conduction) equations are 
to be solved over a region requiring a large number of sub- 
divisions. For example, for a three-dimensional problem 
with N interior points in each direction there is a total 

3 3 


Ort n3 interior points, hence an NX x N° matrix must be solved 


for each time increment. The procedure becomes obviously 
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ampractical if N exceeds, say, about 10. 

Peaceman and Rachford (1955) introduced an alternat- 
migecirection implicit (ADI) method for use in problems 
mivolving a large number of internal nodal points. In this 
method the size of the matrix to be solved at each time step 
is reduced at the expense of solving the reduced matrix 
many times for each time step. For example, referring to 
the above problem of a three-dimensional region the ADI 
method transforms the problem to solving an N x N matrix ne 
times for each time step, which is easier than solving an 
Matrix ata tyne. 

An alternating-direction implicit procedure requires 

a line-by-line solution of small sets of simultaneous equa- 
tions that can be solved by a direct, non-iterative method. 
Analysis of the procedure shows it to be stable for any size 
time step and to require much less work than other methods 
that have been studied. As a practical test, the new pro- 
cedure was used (Peaceman and Rachford, 1955) to solve the 
heat flow equation with boundary conditions for which a for- 
mal solution is known-* The two solutions were in good agree- 
ment. Also a rapid convergence for the solution of Laplace's 
equation in a square was obtained using a suitable set of 
iteration parameters which were easily calculated. An analy- 
sis was presented (Peaceman and Rachford, op. cit.) that 
showedthe method to require about (2 log N)/N as many cal- 
culations as the best previously known iterative procedure 


for solving Laplace's equation, where Ne is the number of 
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points for which the solution is computed. 


In obtaining the usual implicit differenc equation 


for 


both the second derivatives are replaced by second differ- 


ences evaluated in terms of the unknown values of u. 
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Brea mesh”is’choséen' such that Ax’ =' Az: 

Thus the unknown excess pore pressure at the time step 
k may be solved for implicitly by the above equation, pro- 
vided the excess pore pressure at the time step (k+l) is 
known. 

Large sets of simultaneous equations are formed, 
mercn can.be solved practically, .only,by.iteration... If, 
however, only one of the second derivatives, say au pus 


dZ 
replaced by a second difference evaluated in terms of the 


2 
unknown values of u, while the other derivative 2 a *1iS 

3 
replaced by a second difference evaluated in Reaee of the 
known values of u, sets of Siiultaneous equations are formed 
that can be solved easily without iteration. These equations 
ere implicit in the z-direction. If the procedure is then 
repeated for a second time step of equal size, with differ- 
ence equations implicit in the x-direction, the overall pro- 
cedure for the two time steps is proved to be stable for any 


Size time step. Thus two difference equations are used, one 


mor the first’ time step, the other for the second time step. 
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These equations may be arranged in the following form 


Suitable for calculation. 
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Use of each of the above equations at each time step 
Weads* to) Ne setsvof, N. simultaneous’ equations... The’ general 
Bracedure for solving these equations, the criteria for stabi- 
lity and convergence and related phenomena are fully discussed 


by Peaceman and Rachford (op. cit.) and Douglas (1955). 
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42 
4,7 BARAKAT AND CLARK METHOD 


Barakat and Clark (1966) described an explicit differ- 
ence scheme which is unconditionally stable for the solution 
of the multi-dimensional, time-dependent heat~-conduction 
equation. The method possesses the advantages of the implicit 
scheme (i.e., there is no severe limitation on the size of 
time increment), and the simplicity of the explicit scheme. 

Barakat and Clark (op. cit.) examined the finite 
difference representation of the following boundary value 


problem of heat conduction, 


[AA 
ed 
JE 
a 
DO 


[AA 


Ue = Oe > 
2 POA Te 


Subjected to prescribed temperature at the boundaries and 
Pond. prescribed. initial, condiLion. 

To solve the above problem with finite difference, they 
considered two auxiliary functions V(i,j,k) and W(i,j,k) 
satisfying the following difference equations, 
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and, a Simi far expression for W(i,j,k) with subscripts k and 
(k+1) interchanged on the Teftehandysidecofreqdation <4. 13: 
Assumimngenedy Ff ok tandaudisy, kK) cfanctians also satisfy 
tie boundary and initial conditions for the heat eonduction 
problem, the temperature U(i,j,k) at any time level k may be 


maceneas (the artthmetic “average sof Vii pp sk) and WCTSp, Ry Pei.e., 


CI ieee Oi tia Kk) re WOT 5 Kn Bald 


If finite. difference equations. for V{(i,j,k) and 
W(i,j,k) functions are stable the solution of the equation 
mez will be stable. _ It has been found that the finite 
Mmrrerence solution for 4.12 is unconditionally stable for 
mimevalues of At, Ax, Az. 

Equation 4.12 was solved for the following boun- 
dary and initial conditions employing both the present method 
and the ADI method to study the effect of the form of the 
difference equation on the relative accuracy. Tne initial 


and boundary conditions are; 
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mre’ results” are giver in Table 4.4% “Both the numeri - 
cal metnods give results of comparable accuracy. ss boTcn, tne 
methods are unconditionally stable for any time step inter- 
Gal ADI method requires’ the solution of a ‘tridiagonal 
matrix whose solution can be obtained by an algorithm 
derived from the Guassian elimination technique. The pre- 
menu method possesses the simplicity of an explicit method 
and employs the same "marching' type technique of solution. 
Thus the present method is simpler to formulate and easier 
momprogram. fable 4.2 gives a summary of the comparison 
between the two methods. 

It is interesting to compare the machine time required 
by the present and the ADI methods for two dimensional pro- 
blems. The ADI method needs to solve N sets of N simultaneous 
equations for’ N x N nodal points at each time step. And the 
present method requires the solution of N explicit equations 
N times for the same number of nodal points at each time 
Step. As illustrated in Table 4.2, the present method takes 
less time than that taken by ADI. 

The one serious drawback of the present method is 
that it requires at least twice the storage requirement of 
fre ADI (Table 4.2). This is inherent in the method of 
Mamoron since a single function U(i,j.k) is) split. into 
two components V(i,j,k) and W(i,j,k) and both components 
are to be stored in memory in order to evaluate U(i,j,k) 
by equation 4.14. For many of the problems considered in 


determining the excess pore pressure distribution the storage 
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requirement by the present method will be excessive for a 
large number of time steps. 

It may be mentioned here that the present method can 
be extended to handle problems having three dimensions which 


is not convenient for'the ADI method (Barakat and Clark, 


op. cit.). 
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Table 4,.1la Comparison between ADI and Bara- 
Kat and Clark Methods: Percentage 
pore pressure values at the im- 
pervious boundary on the centre 


Tine. 
7A Bae Tos ORAS | MX =e2 NZ = 2 
Time Factor ADI Barakat & 
Clark 
0 160.0 100.0 
OPs0i 98.1 99.8 
0202 Shoe! 92.9 
O03 88.6 ERS 
0.04 Soe). Sard 
O05 192 (Ono 
0.06 Lee 74.7 
0.07 jl | Talece 
0.08 68.5 68.0 
0.09 Oats Gores 
Oe i0 Onl rey 


Table 4.1b Comparison between ADI and Bara- 
kat and Clark Methods: Values of 
average degree of consolidation 
in percentage. 
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Comparison between ADI and 
Barakat and Clark Methods 
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ADI AND BARAKAT AND CLARK METHODS 
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CHAPTER V 


TWO-DIMENSIONAL CONSOLIDATION 


A. GENERAL 


When a soil consolidates under an applied load drain- 
age and compression frequently take place in three dimensions 
(as opposed to the assumption of one-dimensional drainage 
and compression in the classical theory of one-dimensional 
consolidation due to Terzaghi). The present state of know- 
ledge of the compression characteristics of soil in three 
dimensions, however, is meagre; and the displacement of 
Principal importance. In structural applications is» the 
Meortical.compressdon. (With; due recognition; of) this’ fact, 
the effect of compression in directions other than the ver- 
tical may often be neglected. However, the process of drain- 
age of water in three dimensions does have a considerable 
Siuect on the yate of settlement of a soil loaded at tne 
surface. 

Ine actual rates at which foundations on clay settle 
(published evidence) are generally faster than those predicted 
by the one-dimensional Terzaghi consolidation theory. In 
Some instances the differences are very large. [In many 
Mractical cases, it is obvious that the geometric conditions 
are far from being one-dimensional and that horizontal dissi- 
pation of pore pressure (not accounted for in the Classical 


theory) must make the actual rate of settlement (i .6.5 rate 
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of dissipation of pore pressure) more rapid. Rowe (1968) 
has demonstrated that the real drainage behavior (and hence 
mene pore pressure dissipation) of a deposit as a whole 
depends on the geological details of its formation. Veins 
fmeotlt along fissures, or organic inclusions can affect 
the permeability of the mass to a large extent and hence 


Pies drainage characteristics: 


9.A.1 TERZAGHI - RENDULIC . TWO-DIMENSIONAL THEORY 


The simple consolidation theory generally attributed 
to Terzaghi (1925) and Rendulic (1937) considers the three- 
dimensional diffusion of water through a porous medium. 
bom an isotropic soil in which water flows in.all three 
eoordinate directions, such a cansideration leads to a single 


e@uation 2.12 whichis reproduced as 
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Equation 5.1 is a half-way solution between the classi- 
cal one-dimensional theory and a true three-dimensional 
theory. This equation retains all the characteristics of 
one-dimensional theory and is based exactly on the same 
assumptions, except that it is more complex than the classi- 
cal one-dimensional theory. Equation 5.1 provides a means 


for dissipation of excess pore pressure in all directions 
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though the compression is restricted only to the vertical 
aarect ion. 
Under™pTane™ Strain conditions, equation 5.1 reduces 


mo equation 2.l4-which is 
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fens equatton has all’ the characteristics’ of the classical 

One-dimensional theory and yet provides a means for account- 

ing the dissipation of excess pore pressure in the horizon- 

pan direction. 
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Though equation 5.2 is more complex than the one- 
dimensional consolidation equation 2.15 and less complex 
mam €GUatTIONn. 5.1... it has. noti lost. its versatility andimay 
be put to use for all types of problems solved using classi- 
Cal one-dimensional theory. Equation 5.2. is called the 
Terzaghi- equation for. consolidation in two dimensions. The 
Solution of this equation for given boundary and initial 


conditions is presented in a subsequent section. 


OAc ane Nol ON OF DAVIS RAYMOND THEORY T0 
TWO DIMENSIONS 


Davis and Raymond (1965) proposed a modification to 
the classical Terzaghi theory for overcoming the inaccuracies 


of the following assumptions: 


ae 


=) 
- 
J 
a) & 
“ 
= 
Ae 
SS 
pea 
’ 
' La, 
tat 
Me @ 


(POG 


Od) Tihvopi 


iS 


pareg2 avb- a) Vioeas ineanet teat venta 
ee We - an eney. 
Tal Fag fort, pitiwol tot ait: 
YF ee a 


n 
rt fight) ee 
if yf 
4 
. 
ae | 
A & 
a 
& 
= 
! 
fi 
< 
ty 
i ' 
4 


EI oh 


2 


reat 


Pk 


- CIVAG 30 HOT ERS Va ‘Sh | f 


en i ‘a a 
ea eas 
i; YY 
% 
ut 
rod 
=i sf 
; 
aL 
* F 
‘ 
be ae * ND ok. 
¢ Se m ’ 
i 
~~ _— 


e002 ante ae 


>. ae ‘<) odie PS We 


oan tie | 


ie 
cee i é st ‘i hth 7 . 4 


cies 1hene uot seu oF 


iil \vehaatige ito eae ed 


ww 
i ann vert ; seh anipa 
re its ata ¥ ey ad 4 not 


wt babooteng a engl 


enn rananty. OW Te 
a a 


oY ‘o wae 


| 
(qamyed ef) Wom tna mie hahg 2, 


(a) J thecoeffiicient of permeability k is a constant 
during consolidation under a given stress incre- 
ment; and 

(b) the coefficient of compressibility my Tsai S Ona 
constant during consolidation under a given stress 
increment. 

For a normally consolidated clay, Davis and Raymond 
nave pointed out that, during a stress increment the coeffi- 
cient of consolidation Cy is usually observed to change imuch 
fess than the “coefficient of compressibility. They assumed 


Cy GOmDe “al Cons tant = Yt se). 5 


Cie ener constant. Do 


This is seen to be equivalent to assuming that, as 
the soil particles move closer together, the decrease in 
eperficient of permeability kK 1s proportional to the décrease 
in volume compressibility my 

Further, instead of assuming that ay iowa CONSLanicr, 


they adopted an empirical relationship between the void 


Maelo @ and the effective stress co” as 
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where A is a constant whicn depends on the initial void ratio 
eg and tne compressibility index Coe 

With these modifications, Davis and Raymond (op. cit.) 
have developed a nonlinear theory of one-dimensional con- 
solidation. In the following section a three-dimensional 
theory is developed on the basis of the Terzaghi-Rendulic 
assumption with the modifications suggested by Davis and 
Raymond (op. cit.). The development is exactly on the lines 
adopted by Davis and Raymond (op. cit.). 


The coefficient of compressibility Mm, 1S given by 


ieee 
ey Oe, 307 9.4a 
i +o! +g! 
where 8, eee - ae 5. 4b 
since all round compress ony lis) considered. 
Substituting equation 5.3b in equation 5.4a yields 
f 0.434 Cc. nas are 
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where A = AGETRINE Q. 
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in che eqtation S50 thesterm (Ite) tts assumed) con= 
meant.) Ihus the change in woid ratio e is considered too small 
to affect the total volume (Ite). 

Darcy's law is assumed to be applicable and hence in 


general 
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where v., denotes the velocity of flow in the x-direction 


kK, denotes the coefficient of permeability in the 


x-direction 
1, denotes the hydraulic gradient in the x-direction 
u. denotes the excess pore pressure 
and y, denotes the unit weight of water. 


Puastt curving equdtion Dasa m-equation 5.6 end different tate 


ing the resulting equation with respect to x yields 
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Similar expressions such as 5.7a can be written in the y- 


and z-directions. They are 
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The soil is assumed to be completely saturated and 
the pore water and soil particles are incompressible relative 
to the soil skeleton. Therefore the effective stress at any 


point may be written as 


The vertical strain developed in the soil skeleton 


is e whicn may be expressed as 
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where ey denotes the void) ratio. corresponding to zero strain 


and effective stress Og: 


Using equation 5.3b, the above equation may be written 
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Differentiating equation’ 5.8b with respect. to 


and making use of equation 5.5a yields 
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The equation of.continuity of flow demands that the 


rate of water lost through an element dxdydz should equal 


the rate of volume decrease within that element. 
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Tnerefore, 


Making use of equations 5.7 and 5.8c, the above equation 


may be reduced to 
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Equation 5.9 is the general three-dimensional con- 


Solidation equation for the assumptions made in Davis and 
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moymondytneory, jAs in their original theory, no assumptions 


have been made as to the type or rate of loading. 


Special Cases 
CASE 1 


To derive the form of the original Davis-Raymond one- 
. Z Z 
dimensional theory, we may put 24 = Q0 = z : Pen OND Cm LOW 
aX oy 
in the x- and y-directions is assumed to be absent (and hence 


ee 0 = ay 


ae ay Then the equation 5.9b reduces to 


which is same as given by Davis and Raymond (1965). 


CASE 2 CONSTANT LOAD 


Let the applied load be constant with depth and as 
well with time (as in case of oedometer or thin clay deposits 
when load is maintained the same). Hence (o,, + Gaye Tl 
is constant in the following term 


Which on differentiation yields 
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which is the governing equation for a normally consolidated 


soul under, a constant-load. 


For the plane strain case (the strain in the y- 


direction may be assumed zero), equation 5.11 becomes 
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depending on the specified boundary and initial conditions. 
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mresapovesequation-is identical in-form to,that-of. the ,Terzaghi 
two-dimensional equation (5.2) and can be solved using the 


ADI method described in Chapter IV. 


9.A.3 SOLUTION OF THE TWO-DIMENSIONAL EQUATION 


Equation 5.2 may be solved by any of the following 
means: 

(a) closed form (analytical) method 

(b) electric analogue method, and 

(c) numerical analysis method. 

The solution obtained depends primarily on the boundary 
conditions specified. 

Rigorous analytical solutions obtained for one-dimen- 
sional classical theory (equation 2.15) may be extended to 
this Squepione EGualion 5.2 1s mdentical With) tiesdvimeren-— 
tial equation of heat. flow (Carslaw and Jaeger, 1959). If 
u is interpreted as temperature and c as the diffusivity 
constant, the solutions of the equation of heat conduction 
Sayebe utilized in the theory of consolidation, Several 
solutions of heat conduction problems nave in fact provided 
immediate answers to problems in consolidation. 

Carillo (1942) has demonstrated how the familiar one- 


dimensional solutions of the differential equation of consoli- 
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dation (of the analogous heat flow) can easily be combined 
to furnish simple solutions for several important two- and 
enree-dimensional problems. In finding the solution for 
Sqeation 9.2, Caritio (op. cit.) made use of the following 
theorem: 

oA ih Tee Fi (x,t) Isa solution of the Minear -flow 


equation 
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and Une Fo (yst) is a solution of linear flow equation 


then u = UUs is necessarily a solution of the two-dimensional 
flow equation 


2 2 
3 ' 
ae ere r) 4.) — ou 


OX oy ot 

A variety of other analytical techniques such as inte- 
gral transforms, Green's functions may be used to obtain 
desired results. Integral transformation offers, in prin- 
Ciple, a means for the solution of the equation 5.2. However, 
the mere computation of a set of results presents considerable 


algebraic and arithmetical problems (Abbott, 1960). 
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Another technique that may be employed to obtain 
mer xcal tyesults* forthe equation 5 e2is .the "analogue tech- 
nique. Bernell and Nilsson (1957) have described an electric 
analogue method to solve the equation 5.2. However, it will 
be shown that the results furnished (Bernell, 1958) are 
Mecorrect in section 5.A.5. 

An alternative to analytical and analogue methods 
Mast; described is numerical analysis (Gibson and Lumb, 1953). 
By this the process of integrating the differential equations 
is reduced to solving a series of simultaneous algebraic equa- 
tions. The Alternating Direction Implicit method (ADI, as 
described in Chapter IV) can be made use of for the advant- 


ages cited herein. 


5.A.4 COMPARISON OF NUMERICAL SOLUTION 
WITH CLOSED FORM SOLUTION 


For a cnosen rectangular soil mass and prescribed 
boundary and initial conditions, equation 5.2 was solved 
both analytically and by the ADI method. 

For the analytical technique the solution (Art. 2-8, 
per od, Ozisik, 1968) for a finite rectangle (Figure 5.1) 
is made use of. For the boundary conditions of the third 
kind specified therein, the solution for equation 5.2 in the 


megion.0 sx < a, 0 < 2 < b, t > 0 is 
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5 -c(p¢ Pee )it 
MXe Zt iam e ny + e m n Ki) ae Keeeeze) 
m=1 n=] AL 1 
ik t 2 2 
c a : 
ECE ey, aie e (8, Vee AL Broa t dts} 5.15 
where 
K(B, sx) bea 
ACB sv, st \ia= eee ar teuG ee Rivas) Fah aU) dzs 
SEE Hee G piacere cs 
t - ¥ Ki ines Zs EE Merete) oe 
Ko a Bvt n 2 
K(v sz) . 
ae ana + GAM EKNS a italy ae) de 
3 x'=0 
Snick: f? K(aaxt) tplx' ot!) 
+ ' K( Bre aes xe AD ede 
Ky Bal ah acieey, m 4 
and 
7 = a b 1 i} i I 1 tis 
F(Bys%q) = L2 JP kCBysx!) K(vgs2") FO 2 "Dax! a 


Where summation is taken over all eigenvalues B and ies The 


Kernels K(B,»x) and K(v sz) and the corresponding eigenvalues 
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a and Vecre obtained: trom wable 2.1% pp. 60-51, cf. (Ozisik, 
1968) for various combinations of boundary Zieh ine 

Figure 5.2 represents a rectangular section with the 
boundary conditions specified thereon. For these conditions 


eauvacion ©.15 ~educes to 


© foe) -(min® ' nny, AU 
UE rZet ener et CAT ve wl ie D 
m=1 n=] mnt 
sin(8_x) sin(v,z) cos(mm-1) cos(n7-1) 5.16 


Equation 5.16 yields the excess pore pressure Ui fat. ah tpodmnt 
(x,z) at a given time t, provided the value Ug» the excess 
pore pressure at initial time (t=0) is specified. Equation 
9.16 has been made use of to calculate the values of excess 
Borespressure at points (a/2, b/2) and (a/2, 3b/4): for 
Wectanguiar sections of height to width ratio of 1.5 and 
ue he Values are plotted tn Figure 5.3. 

Equation 5.2 has also been solved by numerical analy- 
Sis using ADI method for the boundary conditions specified 
pierioure 5.2. The initial conditiTonhis*theesame as that 
used for the analytical solution. Results have been obtained 
forthe same points K ‘and P Cmigune (52 30ee tie tresaits 
obtained by both techniques are in very close agreement with 
€ach other as is evident from Tables 5.1 and 5.2. Thus either 


of the methods, analytical or numerical, may be employed in 
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the solution of equation 5.2. 

The analytical method yields a rigorous and complete 
solution, while the numerical method only approximates, to 
ee large degree, the actual.valuess -But then thejanalytical 
Hechod"iselimited insits versatility, because of. the%con- 
Siderable algebraic and arithmetical complexities which arise. 
On the other’ hand, the: numerical method is applicable to. all 
sorts of problems and can theoretically handle any number 
of variations in the parameters of the continuum, limited 
only by the capacity of the computer which is used. Also 
with the advances in computer science the numerical solutions 
are obtained at a faster and cheaper rate. Thus the numeri- 
cal technique has a decided advantage over the analytical 
(closed form) method. This clearly demonstrates the relia- 
bility of the numerical technique that can be employed in 


mroesoluUciOneomesuchnequations. ass5.2. 


beA.5) COMPARESON WITH BERNELL™S VALUES 


Equation 5.2 was also solved using an electric ana- 
logue method described by Bernell and Nilsson (1957). Bernell 
(1958) made use of an idealized dam cross-section (Figure 
5.4) and developed a series of pore pressure time curves. 
Pyeciesewresults, Bernel!]-(opsecit.) demonstrated» therdisst- 
pation of initial construction pore pressures as a function 
of the slope of the core of the dam. The idealized dam had 


an impervious core with varying side slopes. The core was 
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edrroundedtby a»spervious fill with side slopes of 1] on 2. 
The dam was founded on an impermeable base and had a height 
meecusmetersawithvalcrestiwidth ofy4 meters..,The coefficient 
of consolidation of the core material was assumed equal to 
O51 cm@/sec. The core was assigned sideislopes of 10 on 1, 
eeeon, 1,72 -oned}iand 1)on 2... Two cases. were considered: 

Faye checdammhause compbeteds “inelayadrull height. of 

10 meters, and 

(by) tne> danypraiseds tos therfull, height. of 20,meters. 

The initial pore pressure was assumed equal to 100 
per cent, which is often the condition experienced with the 
wet fill technique of dam construction (especially in Sweden). 
The effect of core thickness is easily discernible (Figure 
Omi) % 

The same hypothetical dam with the prescribed boundary 
and initial conditions was considered for solution by the 
ADI numerical technique. Both the cases (a) and (b) mentioned 
in the previous paragraph were solved and the pore pressure- 
time curves were drawn for the same assumptions of material 
Pipert ies: As “lisw-eyv iden tjonrrometatqune: 5258 mative 
respective pore pressure time curves do not have any agree- 
menceatiia ldereiBerne lidsy curves are flat and are almost straight 
econ semis logamiuthmic. splot Jacl tsa st Never whe case ani prac- 
tice. » [m almost’ all cases, theidissipation of pore pressure 
is rapid at the initial stages because of the development 
of enormously high hydraulic gradients. As dissipation 


continues, the value of the hydraulic gradient decreases, 
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and hence the rate of pore pressure dissipation decreases 

and ultimately when the excess pore pressures are very small 
eieerate Of®dissipatton ‘also will be small, or in other words, 
the time taken will be large. Thus the shape for pore pres- 
Sure-time plots anticipated are of inverted S-type, the 
familiar consolidation curves. Such S-type plots have been 
obtained by the numerical technique employed. 

The’ Bernell’s curves are almost straight lines and 
depict that the hydraulic gradient is almost constant through- 
out the dissipation of excess hydrostatic pressure over the 
days. This indicates that at earlier stages Bernell's curves 
yield a lower excess pore pressure than those which actually 
exist in the dam while at later stages the pore pressure 
Prcess Values are higher. This is particularly true when 
mewn! heignt of 20 meters is built up (Figure 5.5b). 
nus the greater the height of the dam, larger will be the 
error involved in the Bernell's values. This is amplified 
mle the more in the case of thick cores (side Slopes 5) o0n 1, 
Or more) for shorter times and for longer times in case of 
thinner cores. Thus, the pore pressure time curves obtained 
by Bernell through electric analogue method are misleading 
and yield erroneous results. 

The ADI method used has been tested for its reliabi- 
lity in the previous section. And as such the results 
obtained are trustworthy. At this stage it may be pointed 
that this technique is versatile and can handle any set of 


geometric configurations and boundary conditions. Also it 
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is capable of handling any type of initial excess pore pres- 
Sure distribution in the soil mass. It may be mentioned, 
further, that the algorithm developed can handle any non- 
homogenity within the soil mass and also any arbitrary soil 
Pomavior With time. | That is to say any variation in the co- 
efficient of consolidation Cy (to mention only one soil pro- 
perty) with space and/or time can very easily be handled by 
the program. This ADI technique is the basis for all the 
numerical work discussed in the following sections and 


chapters dealing with two-dimensional problems. 
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B. IMPEDED DRAINAGE 


Op oad elt RODUCT ION 


Embankments of granular or noncohesive materials are 
more stable than those made of cohesive soils because granu- 
lar materials have a higher frictional resistance and because 
their greater permeability permits rapid dissipation of pore 
water pressures. Embankments of homogeneous materials of 
relatively low permeability generally have flat slopes. 

Zoned embankments, which have free-draining outer zones 
Supporting inner zones Of relatively impervious materials, 
are made use of where 

(a) homogeneous materials are not available in large 

measure; 

(b) sufficient pervious material is available; and 

(c) the section of the dam cannot be made too large 

because of cost considerations. 

The two principal requirements for a satisfactory 
free-draining outer zone are that it must be more pervious 
than the protected soil in order to act as a drain and that 
it must be fine enough to prevent particles of the protected 
Soil from washing into its voids. Where the difference in 
coarseness between the fine and coarse embankment zones is 
too great to meet filter criteria, zones of intermediate grada- 
tion must be provided. These are constructed of sands and 
gravels with special gradation characteristics and are called 


mieters. 
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The purpose of a filter is to provide an easy path 
moununehexcessawater to floWeasyitsis squeezed outsof a soil 
layer during consolidation and during steady seepage. Thus 
Some rnectiverti ltercwill vacceleratesthe process of consoli- 
g@eeion.y® In practice, however, such installations. of«filters 
have met with varied success. Two case histories will be 
cited in support of this statement. 

Valajoskoski Dam (Arhippainen, 1964) is a composite 
Sareagtilpsectionewith accentral*yolled,earthfill oftglacial 
moraine surrounded by transition zones of graded filter 
fiecerialland supportedrby»pervioussearthfili zones:y<ltvhas 
been observed that the pore water pressures measured in cen- 
tGalymorainic fill and the outer pervious earthfill are com- 
parable in magnitude. 

Nanak Sagar Dam (Gupta and Sharma, 1964) is a symmet- 
rical rolled’fill zoned section constructed with silty clay 
material in the core and Sandy material with traces of fines 
in the outer casing. Here again the measured pore water 
pressures long after the construction has been completed in 
the core and the outer shell are of comparable magnitude. 

These two field examples clearly indicate the inade- 
Guacy of the performance of some outer shells as efficient 
Filters. 

To evaluate pore water pressures in the core of a 
Zoned embankment after it has been constructed, generally 
it is believed that the filter zones act with full efficiency. 


This implies that the pore fluid at the interface of the 
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core and outer shells remains at atmospheric pressure, 
although for this to be strictly true the filter should be 
Infinitely-pervious. In other words, the degree of dissi- 
pation of pore pressure at any given time depends on the 
relative permeabilities and dimensions of the core Saimaig 
miter’ zones*inecontact with it. =< This’ problem of retarded 
Consolidation is similar mathematically to the heat «conduc- 
tion problem of a slab initially at a temperature, and cooling 


a 


morius Surtaces accordingito Newton's «Taw “oferadiatton. 
analogous to condition 5.19a (Carslaw and Jaeger, 1959). 

The problem of retarded consolidation (or otherwise 
called as Impeded Drainage) is also experienced in laboratory 
testing. Impeded drainage in testing is concerned with the 
Sirects of iporous stones andvfilter drains ion the degree of 
dissipation. In an oedometer test, the effect of fine porous 
Scone can affect «the tesit results and «the coeffierents derived 
memecneabasis oF testiresulhtsay latectrraxtalitests ethnesemhect 
Sierine porousWend capsivand’ of filtercdrains with *finite per - 
Meabilityscansaffect thescoefficient, the “shear ‘strength, 
and the duration of testing inva drained test.. Newland and 
Allely (1960) have investigated the effect of various types 
of filter disk and different specimen thickness on the observed 
time-settlement behaviour of undisturbed and remoulded speci- 
mens of Whangamarino clay. Their results not only empha- 
size the importance of adequate permeability but also suggest 
that overall disk permeability alone may be an insufficient 


criterion, and that special precautions may be necessary to 
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avoid 
(a) the formation of a thin boundary layer of the 
disk becoming clogged with clay, and 
(b) the head drop at the specimen-disk interface due 
to convergence of flow into the pores of the disk. 
Schiffman (unpublished work) recommends the following to 
have effective drainage in oedometer and triaxial tests: 
(i) ee rorssotteclays: coarse stones with t1icer paper, 
(je forestint clays: Tine stones “to avOtds clogging 
and make necessary corrections thereafter. 
Bishop and Gibson (1963) considered the effects of 
disk permeability on the progress of one-dimensional consoli- 
Gation of a clay layer, assuming that the clay is fully 
Saturated and behaves in accordance with the simplifying 
assumptions adopted by Terzaghi. 
In many cases, the field evidence shows that 
the drainage is far from being one-dimensional. For example, 
Gibson (1958) while comparing the calculated and observed 
pore pressures for Usk Dam remarked that towards the end of 
construction the lateral flow of pore water became important 
fmespite of the fact that all the internal drainage provided 
is by horizontal filters at various levels. In this chapter, 
the dissipation is considered to be two-dimensional and an 
ee enent of the influence of impedance of actual drains on 


dissipation is made. 
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debe SANALYSPSCANDORECGHNIQUE SOF SSOLUTION 


Consider the general problem of two-dimensional 
Tne pore pressure dissi- 


(horizontal and vertical) drainage. 


pation is governed by equation 5.2 which is 


meme ne pore pressure is defined in the region 0) < x. <. 2H 
moayur< Zz < L. The general boundary conditions are, 


du 4 . 
A, So (yas Oe ON tack By u(x ,0,t) 0 
(eee ett) SoBe ek. ee) 20 
2 OZ es 2 fees 
du ‘ 
A, 5x (05z>t) + B. Le Orezeate) 0 
du = 0 
Ay aa (Ot eet.) By Useceezeests) 
These boundary conditions can be 
hGeenucrai ning s «it A. = 0 
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Miepequation 5:2 yields the progress of consolidation in 

two dimensions of a clay layer, assumninguthat sthe 1c lays is 
fully saturated and behaves in accordance with the simplify- 
ing assumptions due to Terzaghi. Equation 5.2 may be solved 
moteea given set of boundary, conditions either analytically 
pearitio, 1942) or by numerical methods. (ADI). 

As mentioned earlier, the filters are usually not 
fully effective. To assess the effectiveness of a filter, 
consider a rectangular SOCCLON With Vertical side wrainsso1 
the same thickness on either side (Figure 5.6a). The height 
to breadth ratio of the section may be varied. 

The rectangular (core) section of the clay is of height 
Peandewidth ZH; theedrain on each@sidecis ofvheight L and 
thickness d. The base of the Peettan is assumed impervious. 
The core has a permeability Ky and its coefficient of two- 
dimensional consolidation is Cy. The drain material is 
assumed to be ideally non-consolidating, i.e., tne drain 
material is an incompressible layer. In addition, it is 
mesumed to be fully saturated. As the drain material 1s in- 
compressible, its two-dimensional compressibility may be taken 
Bsvedual to zero; therefore, its coefficient of two-dimen- 
Sonal consolidation is infinite. 

The clay is assumed to have uniform excess pore pressure 
due to suddenly applied load at time t = 0, and thereafter 
Miewexcess pore pressure.u(x,z,t) in) the clay. layer 1s,gover- 

medeby the equation 5.2. The initial condition (t = 0) is 


Mee 0) = Ao, a constant within the entire clay region. 
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Only half the section (Figure 5.6b) need be considered for 
analysis, because of symmetry. The boundary conditions for 


moe clay layer then, -at any instant t > 0, are 


JU Fs i 
ae OT Rate » 240= 20 Siu Bae 
“FO eat Zk 3) elas 
and SUSE eaten a= 10 iy 
fe esc 


The boundary condition at the internal boundary (x = H) 
between the two materials is to be determined. 
The equation governing the dissipation of excess pore 


pressure in the drain material is 


Since an oo. ' hatatton o-ls Ts ‘the Laplacian equa tionsror 
Steady state condition of flow of water through the drain 


material. 


If the drain is efficient, it is evident that the con- 


peecye=+e-since, the drain material is, assumed incompressible. 
T te assumption ofincompressibility leads to the condition 
that the hydraulic gradient in the drain thickness is linear 
with hydrostatic pressure equalling atmosphere at the outer 
boundary of drain. 
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solidation of the compressible core material will take place 
as a result of the flow of the pore water from the core into 
the drains in the horizontal direction; depending on the 

ratio of the height to the breadth of the core some pore Mater 
Will get dissipated from the top since the base is assumed 
completely impervious 

PLomPRigurees<obsuTtits vseencthat there exis ts an 
interface between the two materials A and B. At this internal 
boundary the conditions to be satisfied are: 

(a) continuity of flow be maintained; and 

(b) the pore pressure be the same between the core 

and the drain. 

Water squeezed out from the consolidating material 
(core), in the horizontal direction, enters the drain through 
maemintcerface. thus to maintain full continuity at tne 
interface the amount of water expelled from the core material 
Should equal the amount of water entering the drain. In 
Other words, the velocities of flow in both tne materials 


Smoutd bersame on either. side of the interface. Therefore, 


dU oun 7 5.19 
Ky aX (He Zot) = Ko ax (Haze eo) te a 
where Ky and Ko respectively are the permeabilities of con- 


Solidating core material A and drain material B. 


Uy u,(H,z,t) ee 
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Q> 


Q 


- ey 


path ea TARE phe eae ape’ oftiaelt | 
Bere re nee pot ie 
nah Snoteeeat servos ey ‘bi a 3 


E i et a : ey i y ~~ - re a: 
sappy ane ee Neha ee i ca 


‘a Ww eo 


i tea . aie th | ay OR : it = YO3 he as VTCF ey 
» : Q mi e\ * ; , } j “ it 


r 7 heel Lan ; hs 
Lapin oaddtioree 
; oe, ee ee ee ee é i * 2! ne ee : ra ’ ai } " " ohowee 


y sneh , Fines : 


r rk 


. Br 2) " Le: as Ry at ' eg ai bee cy, aap gl 4 : 
hebvedein Say pee ae Sei hs tt 
| iy 2 Yi hah Hi? Seay > Roles ish b if han ener sal wid 


Hts hw 3 ot WOME et hades it hone Bee 


ve a ; i 


‘eae nas 


<< 
= 
~a 


te. s4, 1 nip Fath cwhehy re 


’ 4 Pins ‘’ . ‘ q = A arse 7 _ pe : 
5 ait f . Ve 5 ; yy 7 f 
; : nn‘ CPI Bi u \ ] 4 is | re Mer At s * we why pi 5 oa 
pore»: "are aS wi, Be hee aie 
7 si 


78 


MA deriving equation 5:18) the material B is assumed «to ‘be 
incompressible and also that the excess water pressure is 
Zero at the outer boundary (x = H+d) of the drain. 

When consolidation takes place in the vertical direc- 
fron, the water\(of dissipation) travels in that direction. 
As such the water expelled ain Unis dissiveation coes not 
enacer the drain; hence.during vertical dissipation the condi- 
mom OF continuity of flow. (at.the. interface) need .not_be 
considered. 

The value of the pore pressure shall be the same at 
the interface (x = H). This is accomplished by chosing a 
Single value at the internal boundary, for horizontal dissi- 
Pe Tohsn.tOr.both thesmaterials,.i.e., u(H,z,t), = uj (H,z,t) 


and then satisfying the equation 5.17. Thus 


OU . UCHez. C) 2 
Ky aX CHezene) = Ko Tee ee Te 


mopenorizontal dissipation. 

During vertical dissipation, it is assumed that the 
value of pore pressure at the interface remains at the same 
value as that obtained during previous horizontal dissipation. 
This is justified by the fact that the equation 5.16 govern- 
ing the flow of water in the drain is in a steady state con- 
dition. Also, during vertical dissipation water squeezed 
out of the consolidating material does not enter the drain 


and hence does in no way disturb the equilibrium attained 
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during and after horizontal dissipation. Thus this condition 
helps in determining the value of pore pressure at the inter- 


face as, 


u(H,z,t),, =)Pu (Hey Zach at At)y 5.19d 


for vertical dissipation, where 


entHisez Us) is the pore pressure value at time t 


H 

during “erizon tal dissapetions= and 

UH Zt), 15 theapore nresstire: Yaite at utime 

te ete tie rng vertical i155 pauion. 

Thus equations 5.19 yield the boundary condition at 
the interface (x = H) between the two materials. Equation 
9.2 may now be solved for the boundary conditions given by 
Sons orhid, o.1/D, ").17cs and 5.19c (for 5.19d) Pouether 
with the prescribed initial condition. 

the solution to equation 5.2 for the conditions men- 


tioned above may be obtained by numerical methods using ADI. 


Equation 5.2 may be expressed in a dimensionless form as 
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Equation 5.20 may be expressed in finite difference 


form for any general nodal point of the finite difference 
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mesh. During the operation by ADI, the dissipation is allowed 
to take place in either direction at each time step alter- 


natively. An impermeable boundary, signifying no flow across 


the boundary, is denoted by au (or au) =) OP Peto irate ani a 


truncation error of 0(ax)? or 0(Az)*, fictitious points are 
introduced in the grid system and the impermeable boundary 
condition (at i = 1) is satisfied by setting 

UE ait) = Ul, Fre) 
Permeable boundaries (at i = I) are represented by setting 


(tla j tg 870 


The boundary condition given by the equation 5.19c may be 


expressed in finite difference form as 
CQ tt) =u (IS Zz, tj) 7 An = ((k5/k,) CEs toy 
Srercu( lols zetjor a(lsz5t)y = (k,/k,) CHY a) Sx Une et) 
mena u(tsz43t)© =) (1/(1+XAx)) eu(1-1,2z5¢) 5.2] 


Where A = (k5/k,) (H/d); is called the impedance factor. 


Hence knowing the value of excess pore pressure if Cltosl ee Zeayte) 


at the nodal point just previous to the interface, the value 
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@aeexcess "pore pressure u(Ijzyt) at’ the interface can be 
@ecermined. Thus the interface condition is’ satisfied. 

jie wupedance factor A, involves the ratio of the 
permeabilities of the two materials A and B, and also the 
ratio of the relative lengths of flow of water in both the 
fmorerials., hus A may be said to be a true indicator of the 
‘effectiveness of the nultersdrain.. The greater the value of 
meeelther the permeability of the drain material is large 
eompeared to that of clay material or the thickness of the 
drain is very small compared to the half width of the com- 
pressible clay layer. In either case the drain is quite 
efficient. And in the case A’ > ~, the pore pressure excess 
at the interface is zero, the drain is termed perfectly 
Pepicient. The smaller the value of A: either the relative 
permeability is low (i.e., the drain permeability is compara- 
ble to that of the core material) or the length of the seep- 
age path in the drain material is quite large compared to 
half width of the consolidating material; and for such small 
Papies of 4. the drain is impeding, and it is not fully effi- 
@ient. the major effect of the impedance to drainage is on 
the excess pore pressure dissipation. Impeded drainage 
retards the dissipation of excess pore pressures, which may 
be detrimental to the safety of the embankment and extend 


the duration of settlement. 
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Pee aoulsl > CAND SD ESCUSS TON 


Invtnis study, the common case considered is a rect- 
egutar= section of height Land width 2H provided with side 
Peoins Of thickness d each on either side. ~ The entire section 
is founded on an impermeable base. The drains are considered 
to be fully effective (i.e., the impedance factor A > ~) to 
muaty impeded (ive., 4 + 0). The impedance factor Xd is 
varied over a wide range of values differing mainly by an 
PewmOornagmucude. e@.g., 0.0, 0701, 0.1, f.02°10.0 andco 
for given rectangular sections. The cases considered are 
for various values of height/width ratio of the consolidating 
Mrmevrete 9.07 5-0,02.0, 10 0.5,--and) 0.25..." VIhe nesuits are 
Maosented in Figures 5.8 to 5.14. The average degree of 
consolidation is calculated for the x-z plane at each time 
Step... Ihis is accomplished by first integrating the values 
of pore pressure along one axis (say z-axis) and obtaining 
Ordinates at nodal points along the other axis (x-axis). 
These values are then again integrated. Thus tnis double 
integration at yi time step compared to the double inte- 
gration over the entire plane at the initial time yields 
the average degree of consolidation at that time step. The 
integration in either direction is carried out making use 
of the Simpson's rule for the largest even number of consecu- 
tive complete subintervals in the range of integration; the 
remaining subintervals, if any, are evaluated according to 


the trapezoidal rule. This may be termed as the 'two-dimen- 
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Sional integration'. 

Invonesdimenstonal; consolidation: probtems, the dissi- 
pation is allowed to take place only in one direction along 
a chosen axis. At any time step the average degree of con- 
polidationealong that axis is»calculated.. This*may be termed 
esm One-dimensional integration’. . But usually in plane strain 
meoes, tne poreapressure dissipation: is. inetwo directions. 
One way to calculate the average degree of consolidation in 
Such cases is by two-dimensional integration. Another is 
to obtain the values of pore pressure (after two-dimensional 
dissipation) along a chosen axis and calculate the average 
degree of consolidation by one-dimensional integration. By 
miewlatter means the,results could be, compared with» the, pub- 
lished results of one-dimensional consolidation. 

POrnaenieigntotowwidtn ratio of 0725 (FigurerS2ejatne 
eraeceteothA, the impedance factor., isialmost negligible 
This is mainly because of one-dimensional drainage taking 
bDlace in the vertical direction and the impedance at the sides 
does not have any influence on the rate (or amount) of dissi- 
Patrons. Adlsonbecausesof the small height, there,are rela- 
tively high rates of drainage occurring in the early stages 
of the process of consolidation. To fully ascertain the 
effect of impedance, the drainage is allowed to take place in 
two dimensions (horizontal and vertical) while the average 
degree of consolidation is calculated 

(a) only along a vertical section passing through 


the center line; and 
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(b) for—the entire region by two-dimensional inte- 
gracion. 

Figure 5.9 shows that dX does not have any effect on the 
average degree of consolidation calculated at the center line 
by one-dimensional integration. It is obvious from Figure 
ono that the-effect of A on the average degree of consolida- 
tion through two-dimensional integration is small. 

Figure 5.10 gives the average degree of consolidation 
memos the time factor for av height to width ratio of 0.5. 
For this ratio of 0.5, the drainage may be considered two- 
dimensional - horizontal and vertical. The horizontal drain- 
age iS impeded because of the presence of an inefficient 
anadin. In’ these curves, given by Figure 5.10, the effect of 
Bamroeciearty discernible.” Thetsmalier the value of A, greater 
the impedance and consequently the smaller is the degree of 
Giss#pation. However, as the consolidation proceeds, the 
average degree of consolidation almost equals the same value 
merespective of the value of A -at very large time factor. 
myrs* is because, at very large=values of time the hydraulic 
gradient is small; the dissipation is almost complete and 
hence the same value of the degree of consolidation is obtained. 

A height to width ratio of 1.0 yields a family of 
curves (Figure 5.11) for the degree of consolidation against 
time factor. The rectangular section is square, and the 
drainage definitely is two-dimensional. As such it Ts x= 
pected that the impedance (due to the drain) to horizontal 


drainage does have a profound influence on the aout Ob COs 
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Solidation. This is amply evident from Figure sae For 
email values-of A, the rate of dissipation is very slow to 
Start witn compared to that for higher values of A. Similar 
trends may be observed in the cases of height to width ratio 
Bnmec.Ovand 570. 

POG denergntto widthovratiorof 10.0, .theydissipation 
may be expected to be mainly one-dimensional in the horizon- 
mmole ctions To ,ascertain,.1f it.is.so,athe drains «are 
maintained effectively non-impeding (A + ~) and the values 
of the average degree of consolidation (two-dimensional 
integration) are obtained for various values of time factor. 
These values are compared against the classical one-dimen- 
Sional results for the same drainage conditions (Figure 5.14). 
The results agree very closely. Further to make sure that 
the dissipation is mainly one-dimensional, one-dimensional 
integration is carried out at the centre line in the hori- 
zontal direction and the average degree of consolidation is 
Calculated. These values.coincide with those.of the classi- 
ealeone-dimensional.theory (Figure 5.14)...Since it is esta- 
blished that the dissipation is mostly one-dimensional, the 
effect of an impeded drain in the horizontal direction is 
very high (Figure baloae 

From the Figures 5.8 to 5.14 it may be seen that a 
oaluevof 24 =.0.0 or a value .of A = 0,01 or for that matter 
ceyaluesof ).= 0.1 has the same effection,the average,.degree 
of consolidation. For all these values, the average degree 


mweconsolidation versus the time factor curves are almost 
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miessame TOr any value of height to width ratio. This shows 
that for thts range an order of variation in the relative 
permeabilities of the drain and the consolidating material 
does not have much influence on the degree of dissipation, 
ime, tney effectively impede the drainage to the same extent. 
Then again, there is a wide variation in the values of the 
average degree of consolidation for a value of A = 0.1 to a 
value of A = 1.0 to 10.0. The values of the average degree 
Of consolidation obtained for A = 10.0 and those obtained for 
X > are almost identical, which again indicates that an 
order of variation in the relative permeabilities after a 
Percain Value, say »X = 10.0, in the increasing direction 
does not really add to the effectiveness of the filter drain. 
In other words, whether the drain is infinitely permeable 
Seeperteable to an extent such that 4 = 10.0, tne average 
degree of consolidaticn is not affected to any appreciable 
degree. 

Reais or interest to iInvesticace "the value rdf wich 
obtain in some of the existing dam sections. 

McNary Dam in Oregon, U.S.A. (Sherad et al., 1963) 
has a vertical core of total width 50 feet with vertical side 
drains 8 feet wide on either side. The central core is made 
up of fine silt while the side drains are built of sand and 
gravel. Terzaghi and Peck (1967, p. 55) quote the following 
values for the coefficient of permeability of the above 


materials: 
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Veuevaine, silt Sec toaOke cm/sec. 
k for sand and gravel mix - 1 to 107° cm/sec. 


Taking average values of oe? Ci Sec Ovra Ii epscusle sand Os 


emysec for sandy gravel mixture, a value for A, the impedance 


Bactor 1S. Obtained: 


e2 
SUF 2 ule he Feito htogrsesi 2p 


107° 


Boome ce, order of magnitude of A calculated no impedance to 
Orainage should exist and the side drain may be termed as an 
errective filter. 

Nanak Sagar Dam (Gupta and Sharma, 1964) has a central 
core (trapezoidal section) of 20 meters at the base and is 
provided with casing (on either side) having a width of 24 


meters. The values for the coefficient of permeability are 


given 
aroracoresy = eenx 107° cm/sec 
k for casing = 2 x io cm/sec 
Therefore, ARF oO aoe Le 4.17 
| Zea 


From this value of A, it may be inferred that the flow of 


the pore water in the horizontal direction is impeded. The 
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observations from the piezometers located in the casing in- 
dicate that the pore water pressure values in the casing are 
substantial and are of the order of 20% to 25% of the pore 
water pressures in the core. Thus this shows that the para- 
meter X correctly indicates the efficiency of the drain with- 
mie practicalvilimits 

Sherard et al. (1963) state that theoretically a filter 
panepe very, thin. ;However, “from the practicaldsitandpoint, a 
minimum thickness exists depending on the filter material 
available and the method of construction adopted. Also the 
thickness is governed by the fact that there should exist 
weegaps sor iticorporation of segregated material. «Horizontal 
layers can be as thin as possible because they can be easily 
placed. Minimum thickness for horizontal layers are about 
6 inches for sand and 12 inches for gravel, although thicker 
layers are generally specified. For vertical or inclined 
mubters., a minimumewidth.of 8 feet to 10 feet is suggested 
for ease of construction while 12 feet to 15 feet is prefer- 
able. 

Sherard et al. (op. cit.) have noted that the deter-. 
mination of a minimum thickness for the central core of an 
earth dam is not amenable to theoretical treatment. It is 
governed for practical purposes by the following factors: 

(a) the tolerable seepage loss; 

(b) the minimum width which will permit proper con- 

Scythe Caron: 


(c) the type of material available for core and shells; 
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and (d) tne design of proposed filter layers. 
The following criteria represent current field practice: 
(1) Cores with a width of 30% to 50% of the water 
head have proved satisfactory on many dams under 
diverse *conditions: 'Probablytatcorenat this 
width is adequate for any soil type and dam height. 
(2) Cores with a width of 15% to 20% of the water 
-head are considered thin but, if adequately 
designed and constructed filter layers are used, 
they are satisfactory under most Atecunetanoass 
Eianowyetorget anvideatofiwhat should besthe ratios 
Of permeabilities of the filter and the core materials, let 
us choose the following: 
(a) a minimum width of 12 feet for a vertical filter; 
(b) a minimum width (for the core) of 40% of the water 
head, say H; hence the half width of the central 
Corer s- OM2ns 


The expression for id is 


Table 5.3 gives the values for A for various values 
of the head of water impounded (which for practical purposes 
is the height of the dam). It can be seen that for a given 
Neiant of the dam built to conservative standards, the success- 


ful behavior of the dam (and the drain) depends primarily on 
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the ratio of the permeabilities (i.e., for given core and 
pain materials), The greater is this ratio the better is 
the performance of the dam. Further, for a given value of 
the ratio of permeabilities, the value of A is higher for 
high dams when Parad tO, that fOr low dams. ,annus for the 
Same materials a higher dam functions more satisfactorily 
than a low dam. In other words, a more permeable filter 
erain is to. be provided for a low dam when compared to a 

high dam for the same performance. The performance is mea- 
sured in terms of the ease with which water squeezed from the 


fere due to consolidation {is drained through the filter. 
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FIG. 5.1 A RECTANGLE WITH NONHOMOGENEOUS BOUNDARY CONDITIONS 


OF THE THIRD KIND. 
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FIG. 5.2 A RECTANGLE WITH MIXED BOUNDARY CONDITIONS 


OF FIRST AND THIRD KIND. 
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FIG. 5.4 IDEALIZED DAM CROSS-SECTION (Bernell, 1958) 
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Impervious Base 


FIG. 5.6 a RECTANGULAR SECTION WITH SIDE DRAINS 
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FIG. 5.6 b RECTANGULAR SECTION WITH SIDE DRAINS 
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FIG. 5.7 a SPACE GRID FOR MATERIALS FOR A AND B 


FIG. 5.7 b SPACE GRID FOR MATERIAL A ONLY 
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CHAPTER VI 


INFLUENCE OF DISSIPATION ON CONSTRUCTION 
BORE PRESSURES 


Galea GeNERAL 


A theoretical analysis is developed to determine the 
influence of dissipation on construction pore pressures de- 
veloped in earth embankments. Using a multiple lift con- 
Struction procedure and assuming conditions of plane strain, 
embankments with linear, isotropic, and homogeneous material 
properties are studied. The equation 6.4 governing the dissi- 
pation of excess pore pressure developed in partially satu- 
rated soils is solved using the ADI technique which was 
developed and tested in Chapter V. Six field cases are 
Studied and the comparison between the calculated and mea- 
Sured pore pressures is good within practical limits of 
eccuracy.= The Hosts indicate that the construction pore 
pressures are neey sensitive to the assumed construction 
sequence and that in the anlysis the construction stoppages 
and the resulting pore pressure dissipation must be considered 


to make reasonable pore pressure predictions. 


6.2 SUMMARY OF THEORETICAL METHODS FOR PREDICTING 


CONSTRUCTION PORE PRESSURES 


During construction of a rolled embankment temporarily 
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high pore water pressures develop. This development of high 
construction pore pressure presents a major design problem 
from the viewpoints of safety and economy. Numerous methods 
are available to predict the development and consequent 
dissipation of pore pressures during the construction of 
earth dams. Few of them have general application. For this 
reason, many practicing engineers and designers use a rule- 
of-thumb, generally based on past experience, for estimating 
construction pore pressures. Theoretical methods for pre- 
dicting construction pore pressures are based on a minimum 
Of assumptions and these metnods provide an insight into the 
variation of pore pressures with different soil types and 
conditions. In the following the available theoretical 
methods are summarized. 

The first analytical method for dealing with construc- 
tion pore pressures in an embankment was published by the 
U.S. Bureau of Reclamation (1939). The basic development 
Rests on the assumption that the time rate of change of the 
Sum of the volume of moisture and free air in an earth mass 
is caused by: 

(i) flow into the unit volume due to percolation; and 

(ii) change in the free-air volume due to changes in 

pressure and temperature. 
This time rate of change is equated to the time rate of con- 
solidation, and the integration of this equation provides 
the solution for the problem of transient pore pressures. 


Hilf (1948) proposed a solution to the unsteady flow 
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of water with the assumption that no drainage occurs and 
Surface tension can be neglected. Pore pressures in a 
partly saturated cohesive soil mass are caused by the com- 
pression of the two phase pore fluid comprising a mixture 
of air and water as a result of an increase in superimposed 
moad. An expression’has been/derived by Hilf (op. cit.) 
for such pressures by combining Boyle's law for compression 
of gases and Henry's law for solubility of gases and assum- 


ing the soil to be completely non-draining. 
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pore pressure 
De di air pressure after initial compaction very close 


to atmospheric pressure 


A = consolidation: or volume change in percentage of 


initial volume of Soll mass 
V = volume of free air in voids after initial com- 
baction in percentage of initial volume of soil 


mass 


H = Henry's constant of solubility of air in water 


by volume = 0.0196 


Hy = volume of dissolved air in percentage of initial 


volume of soil mass 
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Although the assumptions used in the development of 
Hilf's method limit its applicability, the method eliminated 
comptex*mathematical difficulties: Hilf criticized his own 
assumption of no surface tension. The assumption of no 
drainage is overconservative particularly where internal 
drainage is provided. Despite the drawbacks inherent in 
Hilf's method, it often enables a reasonable estimate of 
the construction pore pressure to be made with the available 
data from consolidation and compaction tests. 

An approximate method of allowing for the various fac- 
eorseiicadingyto therdissipationvofrporerpressures in the field 
has been developed by Bishop (1957). The method comprises 
step computations of pore pressures by Hilf's equation at 
certain arbitrary intervals during which some pressure re- 
lease is permitted so as to simulate the combined action of 
various factors leading to dissipation in the field. The 
various steps in the analysis are shown in Figure 6.1. The 
effective stress-volume change curve is first plotted from 


the laboratory consolidation test. For an arbitrary assumed 
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increment oy in effective stress, the volume change A, is 
obtained. The pore pressure u, corresponding to A, assuming 
no drainage is computed from Hilf's equation and point A, 
obtained. Next, the assumed percentage of dissipation is 
applied, bringing down the pore pressure and increasing the 
Prrective Stress. Let this be represented by By on the curve. 
The pore pressure corresponding to this stage is represented 
by C.. For’ the next stress increment of Tos the correspond- 
ing volume change is Ay. For this stage of loading the con- 
ererons at points By and C, are assumed as initial conditions. 
The pore pressure dissipation is again allowed to take place 
in the same ratio and the process repeated. The points Ci» 
wr -.. @te. thus obtained indicate the pore pressure after 
ailowing for the assumed dissipation rate. Ine ratio in 
which the pressure relief is permitted, known as the dissi- 
Poron factor, depends on soil characteristics, rate of 
construction, and cross-section of the embankment. Even a 
low ae OF this dissipation tactor (as low asel/7) leads 
to considerable lowering of pore pressures: 
Paeudivectivyeds aarestit On diSsipatloneslLUsc ih, 
(b) due to increase in effective stress which leads 

to the application of flatter portion of the con- 

solidation curve for the same range of loading; 

thus corresponding to a greatly reduced com- 

pressibility. 


The step computations are based on the assumptions that: 
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Gillis Ss equatton continues. to be valid throughout 
the step computaticn; 

ied ing dissipation, the mixture of air and water 

drained away has the same proportion by weight 
aS vtiat remaining in 1t, .0r im other words, the 
degree of saturation remains unchanged before 
and after drainage. 

The method has been used for a number of dams having 
different design problems and has provided reasonable esti- 
mates of the construction pore pressures (Snerman and Clough, 
1968). For example in the case of Usk dam (Bishop, 1957) 
this method proved very useful. 

Skempton (1954) proposed a solution for determining 
the pore water pressure set up in a fully saturated soil, 

Pyeomessed in terms.of the major and minor Drincipal Stresses , 
assuming that the undrained soil behaved in accordance with 
elastic theory. To compensate for the discrepancy between 
the elastic behavior and actual behavior of the soil, Skempton 
Fepraced the elastic constants in the expression by a factor 


found in a laboratory triaxial test, the well known A factor. 


Hence, 
MT) eeNees A(Ao, - Ag 3) 
where Au - the change in pore pressure 
Ad. - the change in all round pressure 


Ag, - the change in axial stress. 
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Through development by Skempton and others, a practical 
pore pressure theory for use in predicting embankment pore 
pressure during construction was established. 

Skempton extended his expression for pore water pres- 


Sure in saturated soils to unsaturated soils in which 


Au. = BL Ac, hs A( do, - Ao3)] 


and B is another experimentally determined factor. The value 
of B was found to range from zero for completely dry soils 

eoounity for completely saturated soils. For any given soil, 
the coefficient A was found to vary with stresses and strains. 


Bishop (1954) rearranged Skempton's expression as 


The above equation involves the assumption of no dissipation 
miepore pressures, Since excess pore pressure is a) function 
of applied stress only. Also the factor B was found to be 
not constant. The factor B depends on the pore pressure co- 
efficients A and B which vary with the amount of strain that 
takes place. Further Skempton and Bishop (1954) have shown 
that the value of B varies nonlinearly with the degree of 


Saturation. 


Bishop (1957) noted that the above method gave over- 
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conservative results if allowances were not made for drainage 
in dams. In the case of the Usk Dam, a dam compacted wet 
of optimum with numerous drainage blankets in the core, 
Bishop (op. cit.) found the assumption to be overconservative 
and contended that pore pressure dissipation during construc- 
tion shutdown resulted in a two-fold effect on the value of 
maeeossepore pressure atthe end of-constnuctions. The first 
and most obvious effect is that of relief of pore pressure 
by dissipation. The second effect occurs because of com- 
paction wet of optimum. After a period of dissipation, the 
increase in pore water pressure with increase in stress was 
found to be less than the increase before dissipation. 
Bernell and Nilsson (1957) developed an electric ana- 
logue equipment for the study of transient, two-dimensional 
flow problems in earth dams. The equipment permitted the 
analysis of construction pore pressures at any time during 
the construction of the dam. The pore pressure u at a time 


t can be determined by the following equation for two-dimen- 


Sronal consolidation: 
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where C is the coefficient of consolidation whose value may 
V 

be obtained from a triaxial consolidation test. From the 

Hllustrations presented, it appeared that difficult boundary 


conditions can be adequately reproduced in the equipment. 
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Bernell (1958) reported that the model analogy approach to 
miesproblem transient pore pressure dissipation is less dif- 
Ficult than the complex mathematical approach. However, it 
ieomvpeen demonstrated in section 5.A.5 that the results pre- 
eewived by Bernell (op. cit.) are misleading and not trust- 
worthy. Also with the progress of numerical analysis and 
advances in computer technology more accurate results can 
be obtained. | 

Li (1959) on analysis of pore pressure data from the 
construction of the Quebradona Dam, Colombia, South America, 
concluded that allowance should be made for dissipation of 
pore pressure during construction in predicting pore pres- 
sures. In a subsequent publication Li (1967) made a compara- 
tive study of the influence of various factors on the develop- 
ment of construction pore pressures in three earth dams, 
namely Quebrandona Dam, Troneras Dam and Miraflores Dam (all 
mmcolonbia, South: America). Li (op. cit.) concluded: that 
the pore pressure development is due to the combined effect 
of many factors ctl it would be an oversimplification and 
Subject to serious inaccuracy to estimate construction pore 
pressure without an overall evaluation of all possible fac- 
tors involved. Factors which affect the build-up of construc- 
tion pore pressures are numerous (Sherman and Clough, 1968) 


ema are discussed in section 6.5. 
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6.3 ANALYSIS AND TECHNIQUE OF SOLUTION 


Equation 2.11 represents the consolidation of a soil 
mass when three-dimensional drainage is taking place and 


is reproduced below as 


Q 

<J 

= 
Il 


- 3 aE 6.1 


where 8, is the sum of the externally applied stresses. And 
8, on ay C5 1 02. Herein oF and Oz are respectively major 
and minor principal stresses and C5 is the intermediate 


principal stress. Equation 6.1] may be written as 


for the case of plane strain and where the strain in the 
W-direction is considered zero. The quantity Ao is the 
eiange in the principal strésses at any point. Further, 1f 
meabe asstimed that change: in.principalwstress at. every point 
Bougetne depth of the soi] is equal to the change in tne ‘ex- 
ternally applied load, then it will be inferred that the 

Bare of change in. the external applied load is only time- 


J AC) Ve , 
dependent. Thus the term Ae in equation 6.2 assumes tne 


form ao) . Hence the equation 6.2 reduces to 
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The added weight at the top, which is time-dependent, 
may be due to the variation in the Rageeaes} loading in form 
Of surficial loads... Also load at the top may vary because 
of the addition of supplementary layers and thereby cause an 
increase in the maximum principal stress. In this case the 
top boundary is moving with time (moving boundary problem) 
while in the former case the top boundary does not move. 

The presentation has been restricted to fully saturated 
clays, but the analysis may be extended to partially saturated 
Soils by introducing the pore water pressure parameter B 
(Skempton, 1954). This is necessary to estimate the pore 
pressures set up in an earth embankment during and after 
Sonscnuction. Equation 6.3 1s made use .of to estimate sthe 
construction pore water pressures (and those during post- 


construction conditions) and is expressed as 


io) 
con 
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where the parameter B is the fraction of the water pressure 
set up at any point under conditions of no drainage in the 
partially saturated soil to the water pressure that would be 
set up under the same conditions in a fully saturated soils 


and c' is the coefficient of consolidation in two-dimensional 
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dissipation. The term y is the bulk density of the soil and 
the term dh/dt is the time variation in the thickness of the 
soulelayer that is added at the top. Thus sit is seen that 
equation 6.4 has been derived with the usual assumptions made 
in the classical one-dimensional Terzaghi theory except for 
the following modifications: 

(a) the drainage of the pore water takes place in 

;tWotdirections..ife¥. a leongex-"andirzedimens ions; 
(b) the soil is partly saturated as opposed to fully 
Saturation; and 

oy Pthene “ts addition: of soil layers at theeton,1.e.., 
the load application is time-dependent and the 
ton boundary iS moving. 

It will be seen on comparing equations 6.3 and 6.4 
eiateanyesolution obtained forraypartly saturated®soa) is 
equaliy valid for a fully saturated soil if the term yB is 
replaced by y. 

Equation 6.4 7isvuseful tinspredieting porenuressinres 
eateany instant during’and after constructionsofpan-earth-fill 
dam. This equation yields the excess pore water pressure 
provided certain physical constants such as the properties 
of the material going into the embankment; the rate of 
construction, and the cross-section of the dam together with 
the boundary conditions are known. 

The solution to equation 6.4 may be obtained with the 
use of numerical technique ADI. For this to be feasible, 


the cross-section must be discretized in the x-z plane at 


1) AO a ene eee Be 
? it mes ty Ge Ae | -- ie a of of 


pe chit 24S) rete 

ray sina thal wat 

asat “i weed she i@ O , wintery ah 

si helrerk i « bs. ~w BOTH og Be aa ang aaa awd * eae | 
Peqsiceiss ene ei are ih, ai 


a : 


ea 
- 
= 
a 
te 
4 


ee th, Pau er " a 
a hood ole te) oneyel hoa Vea era fiw et ayo0y) ™ ¥ 
MAG" ink Soorisoqued- sorte “At neh aR vsef er 
uh he simi et oan 
, OITHMEOH on’ +4 Ave vi gk al fr, 
ee Fife bweutedon KbS4ng e- SOF a wf WE bi | 


erin mae 09 98 Pies 2076 sah Se ahha &. not | 


“ha Nak: ri i 
“ay b 


q an 


Portes a’ 2te¢9 en ES a | nays ae fet a at ve ry ‘ett ; 

Pres cddees Reto Heleyevseroe TI TK bias spat sistant 
" 

seeeeta sesan evag’ O20 Ven oy pitahy toh amigas a wt 


ee 
eu , Saye 4c etzZ eh Won ZT e ff Ps yn fevhegag bee 


= 


7 rae eae \andanne sion. gg owt gntog. it 


| athe Vodcqgel meh eds Xe: sass rear, ane ‘bhi 
ie wre + thw oa ‘enond tea 
$42 Giiv Won side OM Yom $e nat kung it rotate 
4 hate «t oF ayat 14s Aah suptaioa? Heh foe 
m ai Whaty <-khieds Ae bas Paes bred! an . 


ii pee 


. , UL 
, 7 ot 


lzs 


any instant of time. In order to achieve this Symmetrical 
cross-sections are most helpful but in practice symmet- 
rical sections are very rarely found. Dams and embankments 
with non-symmetrica]l central cores are more common. In the 
eventsof a non-Symmetrical core, "the «cross-section has to 
be idealed as follows: 

(a) Real boundary conditions expressed in finite 
differences are easy to satisfy if a rectangular mesh is 
used. An idealization (mentioned below) is made use of where 
a rectangular mesh is developed at the expense of real boun- 
Gary econditions : 

The idealization affected is for the upstream and down- 
Stream slopes. This stems from the fact that usually the 
values of upstream and downstream slopes bear a whole number 
relationship. This relationship is used in making up the 
mesh pattern. When the mesh is made up, usually the slopes 
(upstream and/or downstream) are approximated as steps 
(Figure 6.4(b)). This results in widening of the corss- 
See@cron of ‘the dan at icertain sections. “ihe rspace step in 
the vertical direction is so chosen that the height of the 
embankment never exceeds the actual height. Thus the lengths 
of the seepage paths in the horizontal direction is affected 
to a certain extent but care was taken not to lengthen unduly 
the lines of seepage. To assess the effect of this lengthen- 
ing of the seepage path, two different idealized cross-~sec- 
tions for the Seitenoikea Dam (Arhippainen, 1964) were drawn 


and the values of pore pressure were determined. It was 
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found that the pore pressure values agree remarkably well 
(within an order of 0.1%). The line of seepage along the 
vertical direction is not affected and also the weight of 

tie soil above the points of comparison (piezometer locations) 
were maintained the same as in the field. 

(b) The actual rate of construction is usually adhered 
to. However in many cases, the rate of construction over 
few weeks or even few days is not uniform (e.g., Figure 6.3). 
If these rates are taken into calculation, the program be- 
comes unduly unwieldy and may run into computational diffi- 
omnes. “In, order to facilitate computation, the rate of 
construction éver a region of time is chosen so that the 
assumed rate of construction is approximately an average of 
faemactudl rates in) that) interval, of time. Thus instead of 
iewing a large, number of varying construction stages, the 
@ntire construction season is divided into a small number of 
easily manageable rates of construction. The duration of 
construction stoppages is rigidly adhered to. 

(c) One of the most important parameters that must be 
known in calculating the dissipation of excess of pore pres- 
Sures is the coefficient of consolidation Cy of the material 
that goes into the embankment. Alternatively, the coefficient 
of permeability and the compressibility characteris tics sor 
the core material should be made known. Only, occasionally 
are such data available in the literature. If available, 


such information leads directly to the determination of the 
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Frequently, however, only the recorded pore pressure 
readings during and after construction are available, and a 
representative value of c, must be determined. The follow- 
ing procedure has been adopted. 

A particular pore pressure cell is selected with res- 
pect to the boundary conditions. The recorded pore pressure 
readings are considered starting immediately after construc- 
tion is completed. The pore pressure value just at the end 
of construction (i.e., beginning of post-construction period) 
is considered as the datum. Based on this datum, the degree 
Of consolidation at various time intervals is calculated; and 
then the degree of consolidation values for the time intervals 
are compared with the theoretical one-dimensional consolida- 
tion curves for the appropriate boundary conditions . Thus 
the values of time factor for consolidation are determined 
at any time interval starting from the post-construction 
Beciod tds ythe “initial times» Making «ause-cof »the respective 
time factors, 1, the time interval t oe the length of the 
drainage path H, the value of Cy (= i) isicaley kateds. eihe 
value of c, thus obtained is on the assumption of one-dimen- 
sional dissipation. This value gives an idea of the order 
of magnitude of the coefficient of consolidation for two-_ 
dimensional dissipation. However, the magnitude of Cy 
Obtained is not the correct value nor can it be used for the 
duration of the construction period. 

By the above procedure a rough approximation of the 


Order of magnitude of Cy is obtained and it is suitably modi- 
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fied to arrive at an appropriate value. Some of the published 
literature, however, yields the laboratory (or field) com- 
pression characteristics together with the permeability tests. 
This information directly leads to the determination of the 
parameter oF 

(d) B, the ratio of pore pressure developed to the 
principal stress, may be varied. The variation depends mainly 
Pune moisture content: of»the soil’ going: into the damuas 
referred to the optimum moisture content determined in the 
laboratory compaction tests. In a saturated soil the com- 
pressibility of the soil skeleton is almost infinitely greater 
than that of the pore water, and thus essentially all of a 
Stress increment applied to a saturated soil is carried by 
the pore fluid and hence B = 1. Ina dry soil the compressi- 
Dility of the pore air is almost’ infinitely greater than 
the compressibility of the soil skeleton, and thus essenti- 
ely all of the increment! in total stress applied to the dry 
soil element is carried by the soil skeleton, i.e., B= 0. 

In partly saturated soils the very high compressibility of 
air relative to those of water and the soil skeleton results 
in values for the parameter B somewhere between 0 and 1 until 
the percent saturation approaches 100%. 

If the field moisture is above the optimum moisture 
content determined in the laboratory, then the soil is very 
nearly saturated and a value for B may vary between 0.8 and 
1.0. On the contrary, if the material has gone into the fill 


dry of optimum, the value of B may be taken as varying between 
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0.6 and 0.7. The response to a variation in B is reflected 
in the different results obtained and may be compared with 


ecttial performance in the field. 


Cee CNS EY HIS TORDES 


S'Tx ee histories have been studied and analyzed with 
respect to pore pressure development. The pore pressure data 
used in this study consisted of measurements in the embank- 
Hewes during and/or “after construction, as indicated by 
piezometers installed prior to the completion of the embank- 
ments. The dams analyzed together with observational data 
are\listed in Table 6.1. Materials and construction procedures 
varied widely for the dams, as is natural, reflecting the 
various countries and tneir organizations. In the following 
Seevions , a Peeagiied’ analy and the comparison between cal- 
culated and observed results are given for each of the cases 


considered, 


6.4.A OTTER BROOK DAM 


The cross-section of the Otter Brook Dam (New Hampshire, 
U.S.A.) with the piezometer locations is shown in Figure 6.2. 
Otter Brook dam is of rolled fill construction. The height 
is about 133 feet. It is of. homogeneous impervious cross- 
section, except for the pervious fill drainage blanket and 


chimney, and rockfill and gravel slope protection. Both tne 
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upstream and downstream slopes are 1 on 2 1/2. 

On the right abutment the foundation is rock consisting 
of mica schist with some granite and gneiss. The left abut- 
ment consists of a deep deposit of glacial till. Impervious 
Fill was placed directly on bedrock on the right abutment. 

The glacial till embankment material was obtained from 
a borrow area opened in the left abutment above the elevation 
of the top of the dam. It consisted of well-graded gravelly 
clayey sand having 10-20% gravel and boulder-size material. 
The material was placed slightly dry of optimum moisture con- 
tent and the average dry unit weight exceeded that at labora- 
Pory-optimum.. Explorations. during design of the dam had 
indicated that the borrow material would be 2 to 3% wet of 
optimum. However, these explorations were made in winter and 
Spring when the borrow area would have been at its wettest. 
Hettally.+ the. summer-of: 1957- proved toxbe one of) therdriest 
emerecord,.andisubstantial dryingiof the ,material*o¢ecurred 
both in the borrow pit and on the embankment. As a result, 
the material as actually placed averaged slightly dry of 
optimum moisture content and the average dry unit weight ex- 
ceeded that at laboratory optimum. 

Six closed-type piezometers were installed in the im- 
pervious fill during construction. After the upstream slope 
movements (bulging by about 3 feet), an additional five open- 
type piezometers were installed. At completion of the embank- 
ment, pore pressures ranged from 30 to 65% and 38. (toe Spop On 


the pressure of the overlying fill as indicated by the closed- 
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type and open-type piezometers respectively. Since completion 
of the dam, all piezometers have shown decreasing pore pressures 

The actual rate of construction is as shown in Figure 
6.3. This figure has been redrawn from the original given 
by Linell and Shea (1960). The assumed rate of construction 
is also shown in the same figure by dotted lines. The assumed 
ioees Or CONS truction are very close to the actual rates. The 
minor variations in actual rate of construction are smooth 
emeenis idéalization. This is'done so as to avoid programming 
and computational difficulties. 

The assumed cross~section of the dam considered for 
analysis is as shown in Figure 6.4. That portion of the dam 
beyond the chimney drain was not considered. For this 
assumed cross-section, the upstream and downstream slopes 
bear an integer (=4) relationship; and the idealized cross- 
Seetionof the dam2is "shown in™Figure 604 (b)V Pi Thus <it may 
be seen, the idealized cross-section does not differ to any 
appreciable degree from the original cross-section. The 
height of the dam remains the same; while the length of the 
drainage path in the horizontal direction at the upper reaches 
is affected; it has been found, however, the length increase 
is very small. The section shown in Figure 6.4(b) is used 
for the calculation of pore pressures in conjunction with 
the assumed rate of construction. 

The construction pore pressures in Otter Brook dam as 
observed are shown in Figure 6.5. A value for cy Was estimated 


using the method outlined in section 6.3 Cfablerouc) «1. rom 
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maese calculations the value of Cy obtained was 1200 feet 
Ssauare per year. This value is for the after-construction 
period and is overconservative. As such a value of 1500 
Peet Square per year for the coefficient of consolidation 
c, was chosen for the duration of the construction and is 
made use of in the analysis of the estimation of construc- 
tion pore pressures. 

The value of B was varied from 0.5 to 0.7 since the core 
material went in dry of optimum. 

Using these values, the analysis was carried out with 
the assumed cross-section. The pore pressures computed are 
expressed as a percentage of the existing overburden and 
plotted against time (Figure 6.5). The computed values have 
been plotted only for the piezometer location 3A. Also the 
ratio of the computed pore pressure to the observed pore 
pressure is plotted against time. If the computed and the 
observed pore pressures are the same over the entire period 
of observation, the resulting curve should be parallel to 
the time axis and have a value of 1. The obtained results 
fore various values of B equal to 0.5, 0.6, and 0.7 are shown 
in Figure 6.6 for the piezometer location 3A. As may be 
inferred, the B value, which not only depends on the moisture 
content of the material but also on the stress increment 


ratio, operating during construction was between 0.5 aitcm Ui Or. 
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6.4.B TOOMA DAM 


Tooma Dam (Pinkerton and McConnell, 1964) is an earth 
and rockfill structure 68 meters high, situated at an eleva- 
tion of 1220 meters on the Tooma river in the Snowy mountains 
of south eastern Australia. The dam forms a reservoir which 
essentially operates to provide flood retention and diversion 
of the waters of Tooma river. 

The tae Site is located in a V-shaped section of the 
valley with steep abutments rising about 80 meters above the 
river bed. The foundations for the dam consist of biotite 
granite, except for the upper part of the left abutment, the 
spillway and the downstream toe where the bedrock is granitic 
gneiss. In the river bed and lower parts of the abutments 
the dam is founded on sound rock. 

The zoning of the embankment as constructed is shown 
mmrigure 6./. .A description of the materials in each zone 


ise as’ follows: 


(1) IMPERMEABLE ZONE 


The material used in the impermeable zone was a resi- 
dual soil resulting from complete weathering of the biotite 
granite. The density of the material in place was required 
to be not less than 98% of the laboratory maximum dry density. 
The moisture limits specified were. 1% on either side of the 
optimum moisture content, although for the upper reaches of 


the embankment these limits extended to 2% below optimum and 
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1.3% above optimum. During handling and compaction of this 
material, considerable breakdown occurred and Substantially 
different properties were obtained for the material in the 
fill as compared with those obtained from samples tested for 


design purposes. 


(2) FILTER AND DRAINAGE ZONES 


The material placed in the filter and drainage zones 
consisted of selected quarry fines blended as necessary with 
crushed sand to give the required gradation for satisfactory 


filter and drainage properties. 


Pe ROCKET LEE ZONES 


The rockfill was dense, fine grained quartzite which 
was hard and resistant to weathering. It consisted of a 
free-draining, well-graded mixture of rock fragments with 
maximum size about 15 centimeters. 

The embankment was constructed in two distinct seasons 
embracing the 1959-1960 and 1960-1961 summer and autumn 
periods. Moisture conditions in the borrow area for imperme- 
able material were markedly influenced by air temperature 
and humidity and it was only during these periods that the 
moisture content was low enough for placement. 

The complete instrumentation is given by 


Pinkerton and McConnell op. cit.) and installations for 
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measuring pore pressures are shown in Figure 6.7. Typical 
pore pressures recorded within the impermeable zone by 
vibrating wire type gauges during construction are 

given by Pinkerton and McConnell op. cit.). These show 
a characteristic build-up during periods when construction 
is in progress and dissipation during the off-season and — 
eueer completion of construction. 

The Cross-Section. ot Toonacdam sis asmsioWwne in bl gune 
6.7. The assumed cross-section of the core for purposes of 
this analysis is given in Figure 6.8(a). The upstream and 
downstream slopes for the assumed section respectively are 
ievertical,tov0<.225 horizontal and } vertical to 0.9 horizon- 
tal. The downstream slope is a multiple '4' of the upstream 
Slope. The upstream slope is approximated by steps as shown 
in Figure 6.8(b). The number of layers in the vertical direc- 
tion can always be so chosen that the number is always a whole 
himbers multiple of '4'. Under such conditions the upstream 
Slope can be approximated by steps. The idealized section 
(Figure 6.8(b)) is made use of in the analysis. 

To obtain a representative value for Cy post-construc- 
tion pore pressure values are very useful (section 6.3). In 
the case of Tooma dam such data is not available. However, 
"the authors have cited the in-place coefficient of perme- 
moat ylas 0.5" x 10° centimeter per second . They have also 
presented the compression characteristics from which a value 
for the coefficient of volume compressibility is obtained as 


varying between 0.002 and 0.0025 centimeters square per kilo- 
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gram. From these values Cy has been calculated and it ranges 
from 0.25 to 0.35 centimeters square per second. In other 
words, the value of c, ranges from about 8000 to 12,000 feet 
Square per year. These values are conservative since they 
wereobtained from one-dimensional dissipation. 

foscarrysout tohesanelysts tavchotce hasrTto be tmade®for 
the value of c,. The measured pore pressures indicate that 
the pore pressures developed during the first stage of con- 
struction are quite low, and they almost remain at those 
values during the construction stoppage. This behavior in- 
dicates tnat the core material was able to AP Ge are the 
pore pressures quickly during the construction phase and dur- 
ing stoppage and fice cain was slower. Hence it is con- 
cluded that the Cy had a certain value during construction 
and it was reduced during the subsequent period of construc- 
tion stoppage. With this in view and the values for ee cal- 
culated from field data (a value of 12000 feet square per 
year during dona rucdion and 8000 feet square per year during 
shutdown) are chosen for the analysis. 

The rate of construction made use of in the analysis is 
acusitown iin igure 629) This’.is redrawn from the rate of con- 
peruetton supplied in the original contribution. 

The in-place moisture content of the fill was 17.9%. 
The optimum moisture content determined in the laboratory from 
soil samples for design was on the average 18.42. USE ibr aS 
inferred that the fill went in dry of optimum. Hence, the 


Value of B is varied through 0.6, 0.7, to 0.8. 
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The analysis was performed for the idealized section 
and for the material properties and the rate of construction © 
noted above. The computed results were compared with the 
observed values of pore water pressure. There was a wide 
variation in the computed values and in the recorded pore 
Paessures.eas live authors attributed the difference ain) theore- 
tical curves (one-dimensional consolidation) and the recorded 
wanes) dWeiito erratic behavior of the piezometers. Flushing 
of piezometers at and soon after installation may have con- 
tributed to the relatively high initial pore pressures recor- 
ded at some points. | 

The computed results are expressed in terms of pore 
pressure ratio, ns The computed and observed pore pressure 
ratios are then compared (Figuresm64:10"and> Grlils)mat Onecie two 
piezometer locations Py. and Pao: Figures GeclOsfanidasors slp eie Ui 
Strate the variation of the pore pressure ratio with time. 
The computed values are in agreement with the observed vaiues 
moni practical kimtts. “Both the figuves eindicatemetnat 
durcangr initial stages of«construction a value of B between 
feoeand 0.7 may bevoperating. eDuring the nextoconstruction 


Season, an average value of 0.50 for B seems to be most 


appropriate. 


pe4eG NSE LTENOLKEA, DAM 


The Seitenoikea Dam was built in 1960 across the Ema 


River, Finland (Arhippainen, 1964). The main section of the 
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dam crossing the river is 900 meters long, has a maximum 
meightsor. 35 meters. It has a composite earth and rock-fill 
section with a central rolled earth-fill core surrounded 
by transition zones of graded filter material and supported 
by upstream and downstream rockfill zones (Figure 6.12). 
The foundation material at the site was oe ae ie oneo ind ated 
coarse silt and fine sand. Underlying this sediment was 
glacial moraine and bedrock. The silt was excavated under 
the highest part of the dam and the core was founded either 
On moraine or on bedrock. 

The embankment core material which was glacial moraine 
has the following characteristic properties. 

(i) Average dry density in place was 2.01 tonnes/ 
meter?. The fill dry density was expected to 
reach at least 95% of the Standard Proctor dry 
density. 

(ii) Average moisture content in place was 10.22, 
while the average optimum moisture content ob- 
tained in the laboratory was 8.8%. That means 
that the fill went in wet of optimum. 

(iii) The degree of saturation of the placed material 
was in the range of 64% to 82% with an average 
value of 74%. 
The location of the three pore pressure cells is given 
in Figure 6.12. The pore pressures recorded during the con- 


struction period and after the construction are given by 


Arnippainen, op. cit. 
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ineraterof cons truceton, as reported .o has three differ- 


ela Stages: 


July 15 - Aug 10 1960 0.22 m/day 
Aug 11 - Aug 28 1960 0.10 m/day 
Aug 29 - Sept 20 1960 0.41 m/day. 


iomioedata 1S plotted in Figure 6.13. 

The assumed cross-section of the Seitenoikea dam for 
analysis is shown in Figure 6.14(a). The approximation is 
Made so as to result in side slopes bearing a definite inte- 
ger relation to each other. The downstream and upstream slopes 
chosen are 0.8 on 1 and 1.2 on 1 respectively. The idealized 
cross-section of the dam (Figure 6.14(b)) is then worked out 
on the basis of the assumed cross-section. As can be seen, 
the height of the idealized dam section remains the same as 
that of the assumed section while the lengths of seepage in 
the horizontal direction are somewhat longer than the actual 
lengths (of seepage). To investigate if these lengthened 
Seepage lines have any effect on the pore pressure dissipa- 
won. another idealized section for the assumed cross-section 
Wasechosen (Figure 6.14(c)). Pore pressures were calculated 
for both the sections (Table 6.3) and it has been found that 
the differences in values at the points of comparison (pore 
pressure cell locations) are very small. Thus the former 
idealization (Figure 6.14(a)) though it lengthens the seepage 
path to a certain extent, in no way affects the ultimate 
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The recorded pore pressures are quite difficult to 
reproduce. Arhippainen (op. cit.) has quoted the pore pres- 
sure ratios at the end of construction. The pore pressure 
ratio r, at a point (Bishop and Morgenstem, 1960) is defined 
as the ratio of the pore pressureto the overburden weight 


above that point. The values are 


ear non .4 re tO 7, 


u 
cetl no. 5 Ge O20 
cells no. -6 ry 220 029% 


These values are made use of for the purposes of comparison. 
Also the values for pore pressure ratio during the construction 
Stages are retrieved from the published results and are 
presented in Table 6.4. 

To obtain a representative value of the coefficient of 
consolidation Cy for the core material, the procedure outlined 
in section 6.3 was used. The values of Cy obtained (Table 
6.5) for various piezometer locations were examined. The 
Walues thus obtained are for the after construction period and 
are conservative since only one-dimensional consolidation 
was Supposed to have taken place. As such the value for Cy 
was varied from 5.0 to 6.0 meters squared per day. 

The material for the core construction was placed wet 
of optimum and as such the value of B will be in the range 
mie0.,o0 to 0.90. 

Making use of the above values together with the Con 


struction sequence noted above, a detailed comparison with 
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field observation may be made. The pore pressure ratios, 

ry, are computed at several time intervals and compared with 

the observed values, Table 6.4 (see Figure 6:15). From 

Figure 6.15 it may be concluded that a value of 6 meter Squared 
per day for Cy and a value between 0.85 and 0.90 for B give values 
ofppores pressure very close to the observed: values. © The 

values of pore pressure ratios calculated are within prac- 


tacals limits: of accuracy. 


6.4.D MIRAFLORES DAM 


The Miraflores Dam (Colombia), Figure 6.16 1s located 
in an area of a deep seated intrusive formation of igneous 
rock known as 'Antiogueno Batholith'. The bedrock is quartz- 
diorite and is generally found Sine Rally 60 feet or more 
below the ground level. 

The soil varies from a pinkish silt having an average 
in-situ density of 80 pounds per cubic foot to a gray silty 
sand (locally termed ‘decomposed rock') of average in-situ 
@ensityiof about 105 poundssper cubic Toots “The Miraflores 
dam was built with predominantly sandy silt because there 
was not enough 'decomposed rock! readily available. 

For measuring pore pressures, the twin-tube piezometers 
of USBR type were installed (31 numbers). Almost all the 
piezometers recorded consistent and valid readings. It is 
interesting to note that the piezometers were found to be so 


sensitive that it was often possible to tell from their read- 
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ings the fill placement activities such as Starting or stopp- 
ing the fill placing operations and the thickness of the layer 
placed. 

The construction pore pressures developed are given by 
Li (1967) who reported that very high pore pressures 
developed in Miraflores dam. In general, the ratio of pore 
pecssiire to vertical stresses (assuming that vertical stress 
equals the vertical load of the fill) is very high during the 
initial stages of construction and decreases gradually as the 
height of fill increases, as shown by the convex shape of the 
observed pore pressure curves. 

The major part of the earth fill was placed in one dry 
season (3 1/2 months - mid-December to beginning of April) 
Bound the clock.’ "The rate of construction of the dam 71s 
Siowm in Figurev6.17. \The assumed rate of construction is 
aomcirownr in Figure .6: 1/7)" bhevassumed* yates are very “close 
momenesactial rates of “constructton. * The local minor varia- 
tions in’the actual rate of construction are smoothened by 
the assumed idealization. 

Miraflores dam (Figure 6.16) is non-symmetrical in sec- 
tion and is provided with a chimney filter. For purposes of 
this study the portion upstream from the chimney filter is 
considered and the assumed cross-section is shown in Figure 
6.18(a). The upstream and downstream slopes bear an integer 
relationship. The idealized cross-section of the dam (Figure 
6.18(b)) is very close to the assumed cross-section. The 


percentage increase in length of seepage tines is very, sma lil. 
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The average in-place moisture content was 21.5% and is 
greater than the Proctor optimum moisture content, Li 
(1967) which ranges from 19-20%. It is inferred, hence, 
tnat the Fill went in wet of optimum. The value of B may 
petveried from 0.8°to 029. 

Pavabueyhoynthescoethict ent«of sconsolLidation Gy is 
Beeainedion’the; lines outlined in section 6.3. rThe.post- 
construction recorded pore pressures are used and the values 
for Cy obtained are given in Table 6.6. 

Also the core material characteristics and behavior 
are reported by Li (1967). The coefficient of permeability 
waecne Till, material in Miraflores dam is of the order of 
ex In. centimeter per second as determined by laboratory 
fests. ihe compression characteristics of the fill material 
are given in Figure 4 (Li, 1967). From these curves an 
average value of 0.G11 feet square per ton for the coeffi- 
cient of compressibility is obtained. Using these values 
the Beecicdentaot consolidation Cy is calculated and Cy is 
950 x 10°" centimeters Square per second. 

ba ope itenacadculatedpa value for Cy using the 50% 
time factor to compare the calculated (one-dimensional dissi- 
pation) and actual measured pore pressure values. The appro- 
4 


ximate value for Cy is. 1000 x*lOsfi centimeters, squares per 


second, 


Comparing the values of Cy obtained by above methods, 
it is seen that they agree very well in their numerical values. 
For purposes of computation a representative value of 1000 x 


1074 centimeters square per second for Ces 1S Chosen: 
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Using the above values for the pertinent parameters 
the construction pore pressures were determined. The com- 
puted results at the end of construction are used to draw 
pore pressure contours which are compared with the measured 
values (Figure 6.19). As can be inferred from DERG UPeLG a9 
the computed values are lower than the measured values but 
are within 10% to 15% of the measured values. Figure 6.20 
Gepicts the ratio of pore pressure to vertical load of ‘the 
mils! i) -In. general, these ratios iare very high at the initial 
Stages of construction and decrease gradually as the height 


Oia id hhincreases:. 


6.4.E JARI DAM 


ihe. Imdus, Basin Progect GPakistan) sis cthesl argesteisdngle 
water development in the world and its two storage dams, 
Mangla and Tarbela are both exceedingly large. The substi- 
ome ion of Jari. Dam (Binnie, et ales. 1967), in place.of Mirpur 
dyke increased the capacity of Mangla reservoir without any 
increase in the owen tat level. 

The Jari dam (Figure 6.21) is so positioned that the 
axis of the dam is parallel to the strike of the bedrock. 
The core, which is composed of rolled silt, covers a sand- 
Stone bed Oy» which is contiuous across the site. The higher 
Sandstone bed oF is either covered by the core or outcrops 


in the upstream side of the core trench. Thus the continuous 


Clay bed between Oy and bo forms a cutoff connected to the 
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core. It was at Jari that shear zones were first discovered, 
and it was during the redesign that different design para- 
meters for clay along and across the bedding were first used. 
| The instrumentation on the Mangla dam project has been 
designed to assist the detailed study of the construction 
dnd performance of all the major structures. Selected cross- 
sections of the dams have been instrumented in detail; Jari 
dam was instrumented in cross-sections at chainages 161 + 50 
and 171 + 50. Pore pressure values are given for selected 
Boepointsein Figuve TO3M(Biniieretvals } Sop .scit® )< 
theteross-section of Jari dam'is as°shown in Figurec6.21. 
Only the portion of the dam between the upstream filters and 
the downstream drained silt (type B) is considered. The 
assumed cross-section for analysis is as shown in Figure 
6.22(a). The upstream and downstream slopes are 1 vertical 
fonOne horizontal and T*verticaleto OT Si horizontal *respec= 
tively. The upstream and downstream slopes bear a relation 
Of 5:2°and hence an idealized cross-section is made up as 
shown in Figure 6.22(b). By this idealization a rectangular 
mesh is obtained which is suitable for numerical analysis. 
Also as can be seen the ideal cross-section is very nearly 
the same as the actual section except for negligible lengthen- 
ing (or occasional shortening) of the seepage lines. Care, 
however, has been taken not to unduly lengthen the seepage 
lengths. The height of the dam remains the same as the actual 
Section. 


The rate of construction of the dam at chainage 161 + 50 
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is obtained from Binnie, et al., COD. oCEE) 
end 1s reproduced in Figure 6:23.~° The assumed rave of con- 
Bande tion is as depicted in Figure 6.23. As can be seen the 
assumed rates adhere to the actual rates very closely. 

To obtain a representative value for the coefficient 
of consolidation of the silt type A, the detailed pore pres- 
Sure readings recorded were obtained from Little (1969). 
The calculated values (obtained on the lines mentioned in 
meGcton 6.3) of the coefficient af consolidation Cy showed 
a Wide irregularity and no consistent results could be ob- 
maimed (able 6.7). This is thought to be due to a number 
of construction shutdowns in between the construction seasons. 
The values of Cy obtained seem to be small for the silty type 
Of material that makes the core of the Jari dam. On observ- 
ing the measuredpore water pressures at various locations, 
it was found that the development of construction pore pres- 
sures at the initial stages tsi small, thereby sindicating that 
G@issipation is rapid. «At later stages of constructionmit 
has been observed that the pore pressure values maintain 
almost at the same level -as the walues “at ithevinittial stages: 
mies econcludes «that as construction proceededseither the 
value of B is decreasing or the value of c, has decreased 
considerably. In the absence of any other data, the value 
of Cy is decreased in this analysis starting from a value of 
Cc. that yields comparable pore pressure results to the measured 


V 
values. The average values of Cy used in a set of calculations 
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during first 3 months of constructio 
during next 3 months of shutdown 
during next 7 months of construction 
otieung next. 1 Wy7omontise or 


shutdown 

during next 7 1/2 months of 
construction 

dieeinds NEXts | dVizenomuns aon 
shutdown 


duringenext- 3 imogtnsmot construction 


during remaining post-construction 


Season; 


Theresis notrecord avaltlableefor a valde of Bamade tin 


the laboratory. B determinations for the Jari core fill 


material were not made on site. The original calculations 


for the estimated construction pore water pressures refer 


Onmy tO ai pore pressure ratio of 0.4, °and*not to anysspecific 


Pavalue (Little, 1969). 


The idealized cross-section of the dam was analyzed for 


the specified rates of construction and the material proper- 


ties mentioned above. 


Comparisons were made of the computed 


results with the observed data. To fit the computed values 


to the recorded data, a variation in set of values for Cy 1s 


envisaged. 
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The measured and calculated values of pore pressure 
are shown in Figure 6.24 for the piezometer locations 27 and 
28. The figures indicate that the computed results are in 
close agreement with observed values within practical limits 
eeeraccurdacy.- “It way be inferred that a value of 0°76 °for B 
15 Operative for a set of values for Cy 4000 - 3500 - 3000 - 


2000 - 1000 feet squared per year, 


6.4.F USK DAM 


The earth dam situated at the upperreaches of the Usk 
River (U.K.) is 1575 feet long at crest level, 109 feet high 
with a cut off trench of maximum depth 77 feet. The volume 
Of water impounded is 2700 million gallons. 

The dam is made of boulder-clay fill carefully com- 
pacted in 12 inch layers sloping slightly away from tne center 
Mine so-as to avoid*ponding. In’ the center, there is a diaph- 
ragm of puddle clay 6 feet wide at the top and increasing to 
lomfeet at 7the deepest part where it joins the. “spear@=nead” 
@rmtne 6 feet concrete filling “in the’ cutoff otrencn. 

Embankment construction began in April 1951 and proceeded 
rather slowly at first (Sheppard and Aylen, 1957). Weather 
Was an important factor controlling the use of equipment 
and placing of fill was possible for only part of the avail- 
able time. With the comparatively high rainfall of about 60 


inches per annum - half of which can fall during the construc- 
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#100 Season -/it was realized: that. the fill was being com- 
pacted in a very wet state and a watch was kept on the pore 
water pressures. 

Pressure cells were placed in the embankment ifill.as 
the work proceeded. These cells showed that the pressures 
resulting from the placing of the first 30 feet of embanking 
during the summer dissipated very slowly during the following 
winter season and that either the rate of construction would 
have to be reduced or some method found of speeding the rate 
Or dissipation, of pore»pressures. The investigation of) this 
problem, and its eventual solution by incorporating horizon- 
tal drainage blankets in the body of the dam, was the most 
interesting and unusual feature of the scheme (Sheppard and 
meen, @0p. cits). The. subsequent pore, pressure: readings.have 
supported the theoretical anticipations. A detailed compari- 
Son between calculated and observed pore pressures has shown 
that towards the end of construction there was not good agree- 
ment when the lateral flow of pore water became important 
(Gibson, 1958). 

The drainage blankets, which stopped 20 feet from the 
puddle core, consisted of 12 inches of river gravel placed 
on the prepared surface of the previous season's fill; then 
6 to 9 inches of broken stone sized from 1/2 to 3 inches, 
followed by 18 inches of river gravel. The new season's fill 
was directly placed on this. 

The actual rate of construction and the adopted Gfoy 


Calculations) rate of construc tion are given in Figure 
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6.26. The difference in the two rates is because of the thick- 
ness of the drainage layers (two in number at elevations 948 
and 980). In the numerical calculations, the drainage layers 
are assumed to be Chee eee (af mo thickness) nandsealsovare 
assumed not to contribute substantially to the overburden 
weight. | 

The actual cross-section of the dam at chainage 800 + 00 
Hseshown in-Figure 6:25. The cross-section assumed for 
analysis is given in Figure 6.27. Only the downstream half 
of the dam cross-section is considered. The boundary along 
the center line of the puddle core is considered impervious. 
The bottom of the assumed cross-section is on a drainage 
mattress and hence it is considered a Per vious boundary. 

A value for the coefficient of consolidation Ge for the 
fill material is quoted as equal to 11 feet square per month 
(Gibson, 1958). This was a value determined in the labora- 
PORV. «Thissvalue of Cy, is made use in the present analysis. 
Gibson (op. cit.) quotes average values of B and y respec- 
tively as equal to 0.85 and 142 pounds per cubic foot obtained 
from laboratory and field tests. Later field evidence showed 
that B associated with further load increments (after first 
construction season) decreased with increasing consolidation 
(Skempton, 1957). This phenomenon has been discussed in 
detail by Bishop (1957). The value of B used in this analy- 
sis varies from 0.8, 0.85, to 0.90. This is so because the 
field moisture content was 2 to 3% more than the optimum 


moisture content determined in the laboratory (the wet site 
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conditions resulted in the material being placed on the wet 
Side of the Proctor optimum moisture content). 

The pore water pressures at the end of the October 1953 
construction season are computed and drawn (Figure 6.28(a)) 
fora Vatue of B equal to 0.85.. The pore pressure contours 
are mainly parallel to each other in the body of the dam. 
They are nearly horizontal with a close spacing near the 
drainage blankets. Thus it can be easily inferred that the 
Puincipal dissipation of pore pressures*is inithe vertical 
direction. Similar pore pressure contours are drawn at the 
end of the completion of the dam, September , 1954 (Figure 
6.28(b)). These figures also show that the main dissipation 
is in the vertical direction. This compares very well with 
the actual pore water pressures measured. Hence the con- 
tention of Gibson (op. cit.) that ‘there was not good agree- 
ment between the calculated and tne observed pore water pres- 
Sures towards the end of construction when the lateral flow 
of pore water became more important! does not seem to hold 
Booed enGibson (opr cit.) predicted the pore pressure values 
based on one-dimensional dissipation in the vertical direction. 
Figure 6.28 show the calculated as well as the measured pore 
water pressures. The observed and computed values agree very 
well. There were very few pore pressure cells installed near 
the downstream side of the cross-section and as such pro- 
bably the pore pressure contours were drawn the way they 
were presented. It may be concluded that a two-dimensional 


study of the effect of dissipation on construction pore pressures 
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yields an assessment of the pore pressures within acceptable 


accuracy. 


bros DISCUSSTON 


The method of solution presented in this study to pre- 
dict pore pressures during and after construction is found 
mompe™ very useful as the’ results ofthe six case histories 
Suggest. As remarked earlier, the method is applicable to 
two-dimensional problems. The results obtained yield repre- 
sentative values on which basis a rational design of dams 
and embankments can be made. Also the method suggests an 
approach to control the rates of construction and the pro- 
vision of drainage features in a dam or an embankment. 

The ADI method is very versatile and can handle any 
changes in input values of material properties as time pro- 
e@resses. [t also can mad ofany tynev or cross -Seetron (with 
any form of drainage features incorporated) provided the cross- 
Section 1s idealized. 

The one important property of the embankment material 
that should be made known is the coefficient of consolidation, 
c,. The method (section 6:3) "adopted inthis study to 
determine a representative value of Cy has proved to be very 
Good. For example, in the case of Miraflores Dam (section 
Oe4-0) sa value of 950 x 107" centimeters square per second 


for c. was calculated from the available laboratory and field 


data. The value for Cy obtained using the method mentioned 
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miesecttonm 623 ths 1000 .x 107 centimeters square per second 
Thus the values for Cc, compare very well. However, it may 
be noted here that the coefficient of consolidation Cy has 
a different definition (section 2.4) depending on whether 
the consolidation is one- or two-dimensional. As such the 
value obtained for Cy by this process is only approximate. 
All the same this approach yields representative values to 
Give an idea.of the order of magnitude for Cy 

There are a number of limitations tothe analysis pre- 
sented here..-In the process’ of idealization (to discretize 
an embankment cross-section), the resulting embankment used 
for analysis will be slightly different from the actual 
cross-section. Also the paths of seepage, especially in the 
horizontal direction, are lengthened (or shortened, as the 
case may be) at certain sections. By this means, the pore 
pressure dissipation is delayed and the time taken is longer. 
With judicious arrangement of the mesh size this discrepancy 
Can be minimized but cannot be overcome altogether. 

As thesembankment.is. being constructed... fresh layerssof 
soil are being added at the top. At any vertical section of 
the embankment, it was assumed that the increase in the maxi- 
mum principal stress (in the vertical direction) at any point 
is the same and is equal to the weight of the added soil 
layer. Thus the increase in principal stress at all points 
along the depth is assumed uniform. This is at variance with 
reality. Any increase in load at a particular height is not 


felt uniformly throughout the soil depth; thes load. increase 
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is maximum at the top and decreases along the depth. Ata 
certain depth and below it, any variation int the ubodd: atatihe 
top does not have any effect. However, the actual distri- 
bution of the load along a vertical section is hard to get, 
especially in a moving boundary problem where the thickness 
Msachanging. with times. For purposes of simplicity, it is 
assumed that all points (along the depth) in a vertical sec- 
tion are affected to the same extent ete to the weight of 
the soil layer added), i.e., there is a uniform distribution 
Of added stress throughout the depth. 

ie to the addition of load (because of the growing 
Reayens) cat ithe «top, the total principal ‘stress at ahh points 
in the central portion of the dam increases by the magnitude 
Of the added load. Because of dissipation of pore pressure 
Mere will ibe«changes im effective stress in both s vertical 
and horizontal directions. What proportion of the load con- 
eeibutes to the horizontal stress: is not ‘known. 

While developing the computer program for the solution 
of non-symmetrical sections of dams, it was assumed that the 
weight added at the top is operative on all nodal points then 
existing irrespective of their location. Thus the algorithm 
yields Br hn uid results and are on the safe side. 

Thus the algorithm developed in this study solves dam 
cross-sections in a grid form. All the grid points are sub- 
jected to the same load irrespective of its position; at the 
Same time, the lines of seepage are almost the same as in the 


original cross-section(s). In other words, the analysis is 
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concerned with rectangular cross-section embankments as far 
as the load is concerned, while the solution actually treats 
real drainage paths. This approximation could be eliminated 
in subsequent work. 

The observations made from the examples presented in 
section 6.4 are intended to present a new approach for pre- 
dicting pore pressures developed during and after the con- 
struction of embankments. The following observations may be 
made in.this connection. 

There are a large number of factors which influence the 
build up of construction pore water pressures. Because of 
their interdependence, it is difficult to isolate the influ- 
Sice,Ofpany one .factor on the development of construction 
pore pressures. The influencing factors may be summarized 
erence. Of construction. including, cons tructiongs toppages,, 
nature of the core material, the location of drainage fea- 
tures, the length of the drainage path, the overburden weight, 
and the placement-water content. 

The rate of construction definitely influences the 
pore pressures developed during construction. Increasing 
the rate of construction results in an increasing rate of 
pore pressures with overburden weight. Bishop (1957) has 
Shown that if dissipation of pore pressures occurs during 
construction stoppages the rate of inerease in pore pressure 
with subsequent increases in overburden pressure is lessened. 


This has been found to be true in the case of Jari and Mira- 


mores. 
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The higher the value of the coefficient of consolida- 
tion oy Of a core material the quicker the dissipation of 
Bore pressure. [hus the soil type has a certain influence 
One tne DUuTTd up-oF construction pore pressures. 

Where internal drainage is provided adjacent to or 
oatnin the core, dissipation of pore pressure occurs (e.g., 
Usk and Miraflores). This helps in lowering subsequent 
increases in pore pressures with increasing fill height. 

‘The fill height influences the length of the drainage 
path and the overburden weight. Thus the dissipation of 
excess pore pressures developed are influenced to a certain 
extent. 

Sherard et al. (1963) stated that ‘water content at 
which the embankment is constructed has the largest influence 
on the magnitude of the pore pressures which develop’. This 
study, in general, supports this conclusion. For placement 
water contents on the wet side of optimum water content, 
pore pressures rapidly increase with increasing water content. 

At this stage some of the anomalies occurring in the 
field may be mentioned. All values of water content, density, 
etc. reported are usually the weighted mean values found by 
various sampling techniques. Some random variation of these 
quantities is expected. However, in several cases because of 
the erratic nature of the weather during placement, or the 
borrow material significant deviations may be observed. 

Also the type of piezometer used has an influence on 


the measured value of the pore pressure. Piezometer designs 
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have been modified and made more efficient. Bishop et al. 
(1964) have pointed out the difficulties in measuring pore 
pressure in partly saturated compacted soil and conclude 
that much of the Standard equipment currently used for pore 
Paessuire; measurements in rolled fills is not in fact suit- 
able for this purpose. The accurate measurement of pore 
Water pressure in partly saturated soils requires properly 
designed piezometers if the difference between pore air and 
pore water pressure is significant. The importance of the 
error due to inadvertent measurement of air pressure will 
depend on the soil type, the placement water content and the 


height of fill above the piezometer. 
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able 6.35 Values of pore pressure at cel | 
no. 5, Seitenoikea Dam 


B=0.90 c, = 6.0 meter/day 


Time (days) Pore Pressure Values (t/m)* 
Cross-Section J Cross-Section 2 
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14 Oraliez eri09 
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ASSUMED CROSS-SECTION OF OTTER BROOK DAM 
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CHAPTER VII 


FINITE STRAIN ONE-DIMENSIONAL CONSOLIDATION 
SEDIMENTATION 


7.1 GENERAL 


This study attempts a comprehensive solution for 
the estimation of excess pore pressure in sedimented layers 
during one-dimensional consolidation. The original theory 
Of consolidation due to Terzaghi (1923) assumes that several 
SOil parameters such as compressibility, permeability, and 
coefficient of consolidation remain constant throughout the 
paecess Of Consolidation: —In*most of the theoretical vresults 
published so far, it has been assumed that the soil behaves 
linearly, both with respect to permeability, and with res- 
Peetetor?tus«stress strain properties. That iss tovsay, 1t 
has been assumed that the coefficient of permeability is 
constant throughout the layer and remains unchanged as con- 
Solidation proceeds, and that the compressibility or modulus 
bmeclasticitysorethe soil eds aisomconstantey tnekactethesper - 
meability of real soil decreases with decreasing void ratio 
and hence with increasing effective compressive stress. 
Thus in a deep bed of clay the initial permeability will 
be appreciably less at the bottom than at the top due to 
overburden effects and the permeability will decrease during 


consolidation as the stresses transfer to the soil skeleton 
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and compress it. The compressibility of the soil skeleton 
will always vary with depth and time in a similar manner. 

ite iseeertaii’thatyvobserved soilibehaviaor) in requir- 
ing a large number of linear elements to describe it, is in 
fact nonlinear. The treatment in terms of linear elements 
has been necessary because exact mathematical treatment is 
then possible. However, the mathematics become more diffi 
cult as the number of parameters increases. For solutions 
to be of more practical use, they must involve a minimum of 
parameters. It is perhaps timely to investigate the possi- 
bility of introducing nonlinear terms involving the minimum 
number of parameters necessary to describe the essential 
physics of the problem, and to produce sufficiently accurate 
numerical»ssolutions. on-ascomputer. 

The major sources of nonlinear behavior for clayey 
SQOils are (Barden, 1968): 

fake finitesstrain: 

(b) varying permeability; 

(c) varying compressibility; 

(d) compressible pore fluid in partly saturated soil; 
and he) pesecondavywcreep Seqtecustegnonlipedn structural 

V 1ISCOS NEY a 

It tsenassibié boeincludéetall theseavauious nonlineari- 
ties into a single general treatment. However, it is usual 
to consider various common types of soil and to include only 
the dominant nonlinearities relevant to each case. 


Theoretical analysis of the effects of nonlinear 
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behavior of a real soil could be carried out for general 
conditions but so far most results have been published for 
one-dimensional conditions. Davis and Raymond (1965) gave 
the results of a theory which assumes the coefficient of 
consolidation Cy to remain constant, while still permitting 
k and Mm, to vary. Their results only encompassed the varia- 
tion with time, not with depth, so their results are applic- 
able to oedometer tests and thin layers in the field. Barden 
and Berry (1965) analyzed the case of a soil whose permea- 
bility at any instant is expressed as a polynomial approxi- 
mation of the excess pore pressure at that instant. Gibson 
et al. (1967) have produced general equations assuming one- 
dimensional finite strain and any arbitrary variation of 
compressibility and permeability. An extension of the work 
of Davis and Raymond (op. cit.) to include the depth effect 
is presented by Davis and Lee (1969) and by Raymond (1969). 
Haretsky et al. (1969) derive a basic system of differential 
equations for three-dimensional consolidation of a soil. 
Nonlinear relations between stresses and strains in the 
soil skeleton, the dependence of permeability on water Sa- 
uration and the effect of the rheological properties of the 
soil skeleton are taken into account. 

In the case of a compressible saturated clay, the 
assumptions of nonlinearity may involve the following empiri- 
cal relationships: 

(a) void ratio - logarithm of permeability is a 


straight line (Raymond, 1966; Raymond and 


: avin (Amel). a ase Me 2h eB) ef po Tes | 
a ee iohiesg Viite a te a hee mer toa *. 08, vt oa its] 
iz MAT higy-sove At laced gee y At e 4 sis oan ont : 
eee tara: ee Er jw cyt abel 
fees! ¢ loots y silt selma satnantd | 
ROTI SH Gas PR, pee men $9 ateee NOB | “s ne 
| | ‘css Ua | 

ano Gaiimere snatsaupe..1 ptadng veiibomy “SNH (02th 
ae HOTS TEN RAGTS DONG Mabe ete Whey sarin Ea ok 


noeie vtuewerl! Bony. VuReane ene () 


hate 
De AGoW ene We aotenedixe WA Cede nae Bei 
. pea. Wy op ue age Gere. (COUP ‘bles tye rr bos 

Ls P hier sy Vt): tn aor eye steeds 60 Ptah (oogsy | ie 3 
- i : 7 ? ot tran Vo weboebt is ; ” oui i terol 2ndwtb-sareilt oh 
| va irae ih 2 hetde.Was-.qumas’ 22 awonbied: enotastes Me 
‘4 ; ‘ id re, ci twinw ed P ideancay Fo ab ae Vb bare oils: poral 


ae ve fe: 
, ai. 


5 aap “the Laebtocgane | rE steals me etsy one 39° J o0tts ious 


ee * 


a Pd reine) AOD a ogat asie9 iss c4y eae | 
oan aly 


‘ 
. e ott fate teen alas eastqians 8. io! oe9:bay oo “ge 
q ptwot tor oat weviornt NOM a : 4 


ee mS ve > i? 


_ 


ef itt dooming 2 snisrnega i 
' ity aes wank <b 


2 por 
ae iy ae . 


200 


Azzouz, 1969); and 

(b) void ratio - logarithm of effective stress is a 
Straight line (Raymond, op. cit.; Raymond and 
AZ Z20UIZ, TOPUSet ty) % 

Lai (1968) noted these empirical relationships and 
generalized the assumptions so as to achieve any arbitrary 
variation in permeability and compressibility with effective 
stress. the coefficient of permeability k and ’the coeffi- 
cient of compressibility m, are related to effective stress 
o'(z,t) by power laws. These relations in turn mean that 
the coefficient of consolidation oy TSPattuncenonsots the 
effective stress. In all these cases, the results show that 
the rate of pore pressure dissipation in an oedometer test 
is affected by the load increment ratio used, a property 


not predicted by the classical theory. 


7.2 MATHEMATICAL FORMULATION 


[twas mentionedSearliervin sectiong7s li that changes 
in layer thickness with time occur due to the finite strains 
involved in the process of consolidation. This implies that 
the quantity (l+te), where e is the current void ratio, re- 
presenting the thickness of the sample at a given depth and 
time can no longer be considered a constant during consoli- 
dation. 

It is customary, when deriving the equation governing 


a physical process which involves a region of space and time, 
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to fix attention on the physical phenonema which occur at 

any single point of space (Gibson, et al., 1967). More 
correctly the neighbourhood of this point is considered and 
mrs Velads’ to the’ discussion of the events taking place in 

an element of volume. We can arbitrarily fix our attention 
on an element of space through which the medium moves and 

to consider the sequence of events in this element. This 

Pee Euler's approach. Alternatively, an element of mass which 
always encloses the same material particles may be chosen and 
we describe the events taking place in this moving and dis- 
torting element as time progresses. This is the method 
adopted by Lagrange. 

Since in a finite strain consolidation problem the 
thickness of the soil sample is continuously changing, the 
top boundary is always moving. Thus at the top, the boundary 
fouttme dependent.- Since this condition 1s both the answer 
we are seeking and the condition required to obtain the 
answer, this type of problem is normally solved by trial and 
error. This is an impossible task with a partial differen- 
tial equation of second order with variable coefficients 
which governs the consolidation process (Gibson, et al., op. 
cit.). This difficulty may completely be overcome if the 
problem is presented in terms of Lagrangean coordinates. 

In this approach the boundary is always identified and the 
boundary conditions on it introduced into the analysis, | 
although we are anon of its exact location. 


In the Lagrangean system of coordinates, an element 
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rie considered to enclose always the same soil grains, and 
hence the element does not always embrace the same pore 
water. (By the same token, an element enclosing always the 
Same pore water may be chosen). In this system the element 
moves and its dimensions change as time progresses. On the 
other hand, in Euler's approach attention is fixed on an 
element which is assumed not to distort. It may be noted 
that when the displacements and strains are small the distinc- 
tion between the two approaches becomes negligible. Thus 
in small strain theories it is of no consequence which sys- 
tem of coordinates is employed. However, in finite strain 


theory, an element has to be defined rigorously. 


An element may be located at a distance ‘a’ from any 


prescribed datum plane at a given time (Figure 7.1). Let 


the element be at a distance 'x' from the same datum, at any 


Subsequent time 't'. Then, 
xX = cleat I ped 


which indicates that x depends on time as well. Thus the 
element is fixed no more in space as in the case of Euler's 
approach and the dimension dx varies with time. 


(a) The vertical equilibrium of the solids and fluid 


currently occupying the element dx yields 
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where o denotes the total stress 
e denotes the void ratio 
ve denotes the unit weight of water 


and sie denotes the unit weight of solids 


Since the chosen element always contains the same 
Weight of solids this leads to the equation of continuity 


GresOlids as, 


TET otal ae ee eg 7.3 
Ys Ite Ys hrecaet) 2254 : 


where iA is considered to remain constant. 


(b) Darcy's law is generally given as 
v = ki 


where v is the discharge velocity of the water and i is the 
hydraulic gradient. A more general form (Scheidegger, 1960) 
may be used relating the law as function of the velocity of 


the water relative to that of the soil skeleton. Thus, 


Mee Wh grou! pei 
a ae Yw 2X 


where n denotes porosity of the soil skeleton (= e/Ite), 


and u denotes the excess pore water pressure, 


ee denotes the velocity of pore water, 


vy. denotes the velocity of solid particles. 
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In the case of steady seepage, ves 0 leading to the 
familiar form. If the entire continuum moved bodily, ee ee 
and the equation then correctly predicts that this movement 
is not associated with the development of a hydraulic gradient. 


Equation 7.4a may be expressed as 
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(c) Continuity of flow of water demands that the 
rate of weight inflow of water must equal the rate of change 


of weight of water in the element. 


Rate of weight inflow of the water is 
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Rate of change of weight of the water in the element is 
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It 1s convenient io introduce. ati this stage (Gibson 
et al., op. cit.) a new independent variable z to replace a, 


such that 


a . 
z(a) = J oCtosmn da 7G 


ivs implies that a point of the soil) skeleton is: identified 
now by the volume of solids z in a prism of unit cross- 
sectional area lying between the datum and the point (McNabb, 
1960). Clearly, z is independent of time. All the relations 
so far derived may now be expressed in terms of the new 
Vomwapie z.. Thus, 


Equation 7.2 becomes 
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Equation 7.10 becomes 


dz/da from equation 7.6 is 


dz/da = 1/1+e, 


Therefore equation 7.3 is expressed as 
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Equation 7.4b, expressed in the new variable Ze es 
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assuming that the unit weight of water does not alter during 


the process of consolidation. 


The governing equation can now be developed using the 
pardons /./i, 17.85 7.9, and 7.10. Substituting equation 


veo in the equation 7.10, yields 
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at a given time is expressed as 
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where x is considered positive when measured against gravity. 
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pore water pressure. Differentiating with respect to t, 


gives: 


=» d(ho) 
corer teary 7.1lle 


Q 
cr 


Posituting equations 7.lJic, 7.11d, and 7.i11e7in-/.1la yields 


2) ie. de d 
sere =] + [ap (ac) - Shay, f saz] =0 © 7.12 


Equation 7.12 is the governing general relationship 
of the excess pore water pressure as a function of time t 
and space z. | 

As a check, the familiar Terzaghi equation may easily 
be recovered. On the basis of small strains (i.e., void 
ratio e is constant), Terzaghi made the rational assumption 
enat the permeability k, the compressibility de/do*, the 
@uantity 1/lt+e are constant during consolidation. therefore, 


eqiiation /.12 reduces: €0 
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which is the equation derived by Gibson (1958) except for the 
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term 1/(14e9)°. The term 1/(1#e,)¢ arises because the z 
reference plane is with respect to the volume of the solids 


as explained following equation 7.6 which may be expressed as 


dz/da = 1/1te, 


Where ‘a' refers to the soil as a whole from any arbitrarily 


chosen datum. The equation derived by Gibson (op. cit.) is 


Meche respect to ‘a Equation 7.13a expressed in terms of 


the variable 'a' yields 
C —s F ae - ao (Ao) (Pa Ae 


Which is exactly the form given by Gibson (op. cit.). 
If it is assumed, further, that there is no change 


in superincumbent load with respect time, 7.13b yields 


which is the classical Terzaghi equation for one-dimensional 


dissipation of excess pore water pressure. 


BOUNDARY CONDITIONS 


The solution of equation 7.t2 will be unique depending 
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on the boundary conditions. At the top, the layer of =—so11 
iseeaiwayS free to drain and 1s open to the atmosphere. As 
Such excess pore pressures will dissipate almost instantan- 
eously. Thus the top boundary of the soil will have zero 
excess pore pressure. That is, a free draining boundary is 
one for which the excess pore pressure is always zero. There 
is a possibility that the bottom boundary may also be free 
draining. 

If the soil layer is in contact with an impermeable 
rock stratum it will be impossible for water to drain across 
the boundary. Also no movement of solid Dahticles; occurs 


across the boundary. Therefore 


substituting the abeve equations in the equation 7.9 yields 


Q 


<u = 0 7.15 
Z 
Equation 7.15 implies that the flow velocity induced by the 
excess pore pressure is zero. That is, an impervious boundary 
is one for which the excess pore pressure gradient is zero. 
The impervious boundary condition can also be expressed 
in terms of the void ratio e at any given instant of time. 
Equation 7.15 is made use of in. conjunction with equations 
fey 8 sand 7.11b to arrive at the desired result. 


Differentiating the equation 7.11b with respect to z, 


he ¥ 
be se ae ‘tt ant we i 
2 A oy Bg zomby ait: ad Ht age et bn, 
Sdubnwdent deomt Pe rrr ph beak: 12549, | bi i 
lahkes aved 1H ines ont inoue) os ‘oes ‘aul wee 
sratnaed pv RO BOT A bt sae sweats ’ “i , 
wait he *, ‘al at GNWECS ofa eeadat ade ‘As 
i" gaa od © “ Be see gti : ea itn of “tadt. on pe 
cid 7 gy A oe, a gh yal 
, tied Wt we oxi fas 
tog a i satan YH? Lcnlaal ae. it wah 
hd sh sy Same ee were Fae a7 oft oeth) 
» kale vonodt -eatnay oa ‘ 
7 or G a 
ee 
ilahac,) Rees ae As 0 ora saa a 


% x 
7" ‘ Bas ; 7) wie 
esr. cs t ’ pei" 7 . i} 
1 ( Md 

7 > i? : his ‘5 

roll . < fi oe 1 

. ry " 4 =. i 

‘ ‘ r 
; 


ett Vd Beowlint: ar? ati welt aay an sotto! ‘asa 
Mehebn yd 1 gpk iam ng y2t dent ova at urease | 
lunes at Sal beng niberena aie reaps Poke dof 
oe aes | I wel'e ia nplithaca viehawed shan 

emi m ‘gasdcnt Weehe ene Te ottey bio af? 105 

engitovps Aho noo nae ny Mae! 
oly” PRACT binifeeb mi ta eye beard 


a iene Ro hw dot. jy oer) wit 


Mg 


gives: 


Using equation 7.15, the above reduces to 
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Equation 7.15b thus gives an expression for an imper- 
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vious boundary condition in terms of void ated) erty bra oee oo I 
prelds the variation of void ratio with depth at the boundary 
across which no movement of pore water or solids takes place. - 
ivemexpression for de/oz indicates ‘that, at the impervious 
boundary, the variation s; void ratio with depth depends on 
such soil characteristics as specific gravity and the current 
Compressibility. fe 

Asin thin tayersy, If the self weight..of the sorl is 


not considered correctly and the buoyant weight is assumed 


to be unity, i.e., G@ = 1, equation 7.15b reduces to 
Hs 7.15¢ 
OZ | 


at the impervious boundary. 
The free draining and impervious boundary conditions 
form the two extremes of a real situation in which there is 


partial drainage across the boundary. 


INITIAL CONIDITION 


In the case, of sedimentation. etiemsO | ala Vvers lane 
fmadually built up. Since layers are growing with time, 
at the initial time (i.e., say when t = 0) the initial excess 
pore pressure is zero. Also at that initial time, there will 
be no clay layer to start sedimentation with, which means 


that the void ratio at the beginning is zero. A zero void 
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ratio always relates to an effective stress which is not 
accountable. To surmount this predicament, the sedimenta- 
£10n shall] start (at t = 0) with a negligibly thin layer 
Seman initial void ratio Cy» and a corresponding effective 
stress OG: At the next time step (t = 0 + At) a layer of 


thickness Az shall be added at a given void ratio e This 


0° 
will generate excess pore water pressure, which dissipates 
depending on the boundary conditions. 

Miesinittal “conditions thus fare: 

(1) a soil layer of no measureable thickness with 

a definite void ratio ey and a corresponding 
effective stress 9 is provided; and 

(2) there is no excess pore pressure, i.e., the 

excess pore pressure in tne thin layer of void 
ratio Ta is’ zero, 

Equation 7.12can be solved for given initial condi- 
tion and a set of prescribed boundary conditions. Unfortuna- 
tely it is impossible to give an analytical solution to the 
equation under these conditions. However, an approximate 
numerical solution can be arrived at (to determine Gerace 
function of z and t) using an implicit finite difference 
technique. Once the value of the excess pore pressure is 
known, the quantities such as effective stress, void ratio, 
the height of the sediment (which is not the sum total of the 


layers added because of finite strains) and the density of 


the sediment may easily be obtained. 
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7.3 PERMEABILITY, COMPRESSIBILITY, VOID RATIO AND 
EP PEC}IVESPRES S REL ATLONSHIPS 


Kieenis SeetlOns ducts cussion. in general, of the 
possible relationships existing between compressibility, 
permeability, and effective stress is given. 

Janbu (1963) concluded, after an extensive study of 
rocks and clays, that an effective stress-strain law may be 


expressed as 


Where the subscript '0' refers to an arbitrary reference state 
and the exponential p iS a non-dimensional constant. The 
value of p ranges from -1] to 0. Equation 7.16 in its pre- 
sent form is the one adopted by Lai (1968); but is identical 
mimeli respects to that given by Janbu-. When p. = 05 de/do. 

1s a constant and hence the effective stress-strain curve 15s 

a straight line representing a linearly elastic material be- 
Mavdor. When p = -1, de/do' varies inversely as the effective 
stress and a plot of e versus log o' yields a straight line; 
thus representing a normally consolidated material. When 

foe 0.5. a plot Of eaversucac: yields a parabolic curve; 
which may not in particular represent any one material. But 
this value of p = -0.5 represents a state which is somewhat 


in between the states represented by p = 0 (a linearly elas- 


tic material) and p = -1 (a normally consolidated SOnuiemecne 
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value of p = 0 roughly represents the case of an overconsoli- 
dated soil in a limited stress range. Thus p = -0.5 possibly 
represents the middle of a spectrum of soils, one end of 
which is represented by normally consolidated soils and the 
other end being represented by overconsolidated soils. 
momever, Tor. purposes of this study. equation. 7.16 4s.made 
use Of with a value of p = -1 because special attention is 
devoted here to normally consolidated soins. 

Reva luc Of Dac =e yi el dss a. plot. ofseuversis wig) ' 


as a Straight line; the equation of which may be written as 
e2'C i= D-logsa Teall 


where C and D (both positive) depend on the soil characteris- 
mess such as liquid timit, plasticity index etc. Such an 
expiacit relation.as given by equation 7.17 1s necessary for 
solving the nonlinear consolidation equation 7.12 as will be 
Shown later. 

An attempt was made to find a relationship between 
ciempermeability K and the effective stress Go. based On 
experimental evidence. Strangely enough there is bare cae 
published evidence available for such a relationship to be 
evaluated. According to Taylor (1948), 'a plot of the void 
ratio to a natural scale against the coefficient of perme- 
ability to a logarthmic scale approximates a straight line 
for any soil'. Schmid (1957), based on a study of several 


published results, reached the same conclusion. 
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Figure 7.2 gives the variation in permeability caused 
by consolidation. The quantity log (k/lte) bears a Strath t 
immerrelationship with log o' for the soils studied. At low 
consolidation pressures, the curves are not strictly straight 
Tines; though they are perfect straight lines at higher effec- 
myisuresses. For purposes of this study the plot of k/i+te 
Peysus 0° 1S Considered a straight line on a log - 109 scale. 
(The curves are plotted from the data of Normand, 1964.) 

Thus it may be assumed that an experimental curve 
relating the effective stress o' and the ratio k/Ite can be 


represented by 


where q is a non-dimensional constant to be determined by 
emevemcitting., AS in the equation 7.16)the subscript 0” 
merers to the same arbitrary Strained state. 

From Figure 7.2 it may be concluded that a given 
change in void ratio e causes a much larger change in k when 
e is large than when e is small change in the ratio k/Ite 
for a given change in e must approach zero as e becomes 
Smaller. Also the rate at which the ratio k/lte changes 
decreases continuously. Equation 7.18 possesses these pro- 


perties. The first cerivative with respect to o' yields 
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derivative 
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also approaches zero with o' at a much faster rate. 


7.4 ONE-DIMENSIONAL CONSOLIDATION EQUATION 


Equations 7.16 and 7.18 may be substituted in equation 


fei to yield 
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Equation 7.20 is the most general equation describing 
the one-dimensional consolidation process. In this study the 
Ovantities p,q, Kkg/1tegs Ty's and (S27) 5 are considered to 
Deweconstants.. The time factor T is defined in the same 
manner as Terzaghi's except that this time factor T is not 
divided by the square of the length of the drainage path. 

Tne time factor T has a dimension of (length)¢, 

As mentioned earlier in connection with equation 7.17, 
the governing equation 7.19 involves the effective stress oa' 
and as such an explicit relation between void ratio e and 
pmrcetive stress o° such @s 7.17 is essential. 

hercan be shown that fromthe: equation 7.19, stnen form 
of the governing equation as presented by Davis and Raymond 
(1965) for a normally consolidated soil may be derived. As 
assumed by Davis and Raymond in their formulation, the weight 
of a thin soil layer (for an oedometer case)may be taken such 
that 5 Sn aie xe Ac(t) = 0, since no additional’ layers are 
Being added. Since it is a small strain theory 28 ="Q. 
Again, for a normally consolidated soil p= -1. They further 
assumed that the decrease in permeability k is proportional 
to the decrease in compressibility cee hence q = -l. All 


these parameters are substituted in the equation 7.20 to yield 
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which is of the same form as given by Davis and Raymond 

(op. cit.). The presence of the term (1/ite)? may be ex- 
plained on the same basis aS was done for the equation / 13a, 
Equation 7.2la may be expressed to the independent variable 


oe to yield 
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which is Davis and Raymond's equation for a normally con- 
SOlidated soil with the value of the coefficient of consoli- 
dation Cy maintained constant. 

Also the form of the equation governing one-dimensional 
consolidation with finite strains as derived by Lai (1968) 
can be extracted from the governing equation presented here. 
It would be much simpler and convenient to proceed from the 
equation 7.]la instead of the equation 7.12 or equation 7.20. 
Combining the equation 7.1lla with equation 7.16 and 7.18 


fas carried out by Lai), yields 
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which is Lai's equation for one-dimensional finite strain 


eonsolidation. 


7.9 THICKNESS OF THE SEDIMENT 


During sedimentation, fresh layers are deposited on 
the existing layers. As the layers (say) of certain thickness 
are being deposited, the pore water dissipates whereby the 
thickness of the layer changes. The change in thickness 
depends, among other things, on the overlying sediment and 
the everchanging physical properties such as compressibility, 
permeability etc. These changing properties can be related 
by a single variable, namely the void ratio e. The value 
of the void ratio depends on the position of the point under 
consideration and also on time. A change in void ratio under 
the changing overburden reflects the changes in compressi- 
Drlity, the permeability and associated properties. thus 
mae change in thickness of a sediment, or the thickness of 
a soil layer may be designated in terms of the void ratio. 

Because of the dissipation of the excess pore water 
Pieessures, the effective stresses change. In the (case of 
finite strains, the thickness of a layer changes depending 
on its location. Let the thickness of the soil deposited be 
dz at a height z from the reference datum. Further, let the 
oad yatio of this layersat any instant of time be e(z,t)" 


From the continuity of solids (equation 7.8), we have 
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where dx is the current thickness of the soil layer. There- 


fore, 


The value of the void ratio e(z,t) can be evaluated using 
the governing equation for a specified initial condition and 
eosumed boundary conditions. Equation 7.22 may be utilized 
in computing the current height of a soil layer of given 


iiitial thickness. 


feo DENSTARY .OFs THE SEDAN ENT 


Since the void ratio © Varies, from position to posi- 
tion within ia seu and also from time to time, the unit 
weight at any given point in the sediment will change with 
time. Once the value of the void ratio e(z,t) is known, the 
density of the sediment at any depth z and time t may be 


obtained as 
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the sediment is required, the individual values of density 
at various chosen intervals may be integrated by Simpson's 
miberand averaged. This-leads to calculating the total 


weight of the sediment deposited. 


77) SPENT DE DIERERENCE REPRESENTATION 


EXact Solutions Of equation 7.20 for specified boun- 
Momyeand initial conditions are difficult to obtain. tne most 
general method for solving this equation is by finite differ- 


ences. Equation 7.20 is reproduced below as 
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To write the above equation in a simpler form the 


following equivalences are made. Let 
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The governing equation in central differences may now 


be expressed as 


or 


where 


which is 


equation. 
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When L = ] 
-U1(1)*U(O,N+1)+(14+2*U1(1))*U(1,N+1)-U1(1) *U(2,N41) = D(1) 


Since at L = 1 an impervious boundary exists, there is no 


hydraulic gradient; and as such 


U(O,N+1) = U(2,N+1) 
imenefore (14+2*UT(T))*U(1,N+1)-U1(1)*2*U(2, N41) = D(1) 


Thus the difference equations may be written for all values 
Cae except, when L = Ns that is; at the top of*the current 
sediment. 

When. L = N, 1.e., at the top of the sediment at the 
instant t = N, the gradient au equals (U(N+1,N)-U(N-1,N))/2*DZ 
where the (N+1,N) point (Figure 7.3) does not exist at the 
mimemc. ilo obviate this difficulty, 1t, 1s, safe to, assume 
that the hydraulic gradient excess is linear between the mesh 
Bonincs (N-1.N) and (N,N), 1.e. , instead of taking central 


Girnerences, backward differences are taken. © —[neretfore: 


-U1(N)*U(N-1,N+1)+(14+2*U1(N))*U(N,N4+1)-U1(N)*U(N+1, N41) 
= D(N) 
where U(N+1,N+1) = 0, always by assumption | 
and D(N) = U(N,N)+DH+DET(N)*DT-(q*U1(N)/Sigma(N))*((UCN-1,N)- 
U(N,N))**2+DZ*UWW*(G.-1)*(U(N-1,N)-UC(N,N))) 


where DET(N) = DET(N-1) 


Thus the governing partial differential equation is 
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converted into a set of simultaneous algebraic equations 
using finite differences. These Simultaneous equations may 
then be solved in the usual manner. 

Once the values of U(L,N) are known, the values of 
U(L,N+1) can be computed. This pROcedirerors narchingrout:. 


the solution yields results at subsequent time intervals. 


7.8 SUMMARY 


A study has been made in order to ascertain the excess 
pore pressures developed in an environment, when sedimenta- 
tion takes place. The treatment has been on the basis of 
general nonlinear behavior of the soil where the soil pro- 
perties were assumed to vary with depth and time. To account 
for large changes in the void ratio (i.e., finite strains 
were allowed to develop) the continuity equation for one- 
dimensional fluid flow was derived in a general manner. 

Since large changes in void ratio usually entail large changes 
iaesoib properties of which theecompressibilitytand’ the per- 
meability are the main factors, these were taken into account 
by assuming that both the coefficients of compressibility 

and permeability are related to effective stress by power 
laws. These assumptions mean in turn that the coefficient 

of consolidation is also related to the effective stress by 

a power law. By combining these with other assumptions made 
in classical Terzaghi theory together with the vertical 


equilibrium of soil water system including the self weight 
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of the soil, the equation governing the process of consolida- 
tion was derived. 

The consolidation equation obtained in this SUV een 
general, is nonlinear. From this general nonlinear equation 
otner Auster equations governing linear as well as non- 
linear consolidation behavior of soils can be derived (as 
demonstrated on previous pages) provided the appropriate 


assumptions that have been made are incorporated. 


eon MANGA SE SHES TORY 


Subject to some assumptions, the observations of Fisk 
and McClelland (1959) on the deltaic deposits of the Continen- 
tal Shelf off Louisiana can be reproduced using the theory 
described in the previous sections. A brief description of 
the area is presented below. 

The search for offshore oil was being conducted in 
shallow waters of the continental shelf along the entire 
Bouisiana Coast (Fisk, 1956). Data from a number of test 
borings, from the gulf floor samples, indicate that the shelf 
surface throughout the region is underlain by relatively 
unconsolidated deposits of the Quaternary age. Laboratory 
investigations of the shelf sediments penetrated by those 
borings have been directed primarily toward the measurement 
of strength characteristics rather than consolidation pro- 
perties. Consolidation tests have been performed for only 


a small percentage of the locations investigated (McClelland, 
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1967), Morgenstern (1967) calculated the degree of consoli- 
dation for a range of rates of sedimentation and coefficients 
of consolidation when the layer gradually grows upon an 
impermeable base. The results revealed that underconsolida- 
tion (average degree of consolidation was of the order of 
Peto 0.50) is significant for silty clays and clays 
eeposited at), deltaic rates of sedimentation. 

The general nonlinear one-dimensional consolidation 
theory Aerenio ved in an earlier section may be applied to 
calculate the void ratio and effective stresses of any depth 
ietcne seaiment. Some basic information pertaining to rates 
Of sedimentation, type of sediment etc is obtained from 
Nectelland (op. cit.) 

The topstratum deposits of the continental shelf off 
Louisiana vary in thickness to a maximum of approximately 
600 feet in the central area south of New Orleans. These 
recent deposits, together with the relatively firm soil 
layer upon which they rest locally, control. foundation condi- 
Puonsaat drilling mint rocmecn neces The nature and distribution 
of the uppermost layers of the sediments on the continental 
slope and in the deeper parts of the gulf have been dis- 
cussed by Greenman and LeBlanc (1956) and appear to be largely 
Clay deposits. 

The principal control on the character of the top- 
stratum deposits has been the Mississippi River, and their 
thickness depends mainly on the distance from sites of river- 
mouth sedimentation in the Aone plain. dhe stream dis- 


charaes more than one million tons per day of sediment through 
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the mouths of its birdfoot delta distributaries. The load 
of sediments carried to the gulf is more than 75 ber cent 
Silts and clays and less than 25 per cent very fine sands. 
The sands are deposited close to the mouths of the distri- 
butaries in river-mouth bars, whereas the silts and clays 
are-carried by gulf currents for considerable distances from 
the delta. 

The nature and distribution of Recent and latePleisto- 
cene deposits and the physiographic features of the shelf 
provide evidence of events which have occurred in the region 
in the immediate geological past. The presence of buried 
valley systems across the shelf, the submarine canyon at the 
shelf margin, and the soil zone which mantles a buried ero- 
sional surface indicate that the shelf was exposed and was 
deeply eroded by stream systems in late Pleistocene time. 

The burial of the erosional surface by the sequence of recent 
Beltaic  depositswiseproof that thessea Tevelesubsequently 
rose and flooded the shelf. The gradational sequence of 
recent deposits from a coarse substratum of sands and gravels 
through a fine-grained topstratum of silts and clays provides 
evidence of the gradual lowering the gradients of streams 
entering the gulf and the reduction of their carrying capa- 
city while sea level was rising. A cycle of sea level change 
is represented by the erosion of the late Pleistocene valley 
system on the shelf area and by the subsequent flooding of 
the shelf by gulf waters. This cycle can logically be 


explained by appeal to the climatic cycle responsible for 
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the growth and retreat of glaciers during the closing stage 
of the Pleistocene. Radiocarbon dates (Fisk et al., 1954) 
Suggest that the modern birdfoot delta began to develop 450 
years ago. The deposition of the leaf-like mass of the bird- 
foot delta has proceeded uniformly and has resulted in the 
seaward elongation of the distributaries and has given rise 
to a pattern of diverging finger-like bar deposits of sands 
and sandy silts separated by wedges of silty clay. These 
sedimentary units, together with an underlying thin layer 
Gieprodelta marine clay, form the framework of the delta 
platform. 

Fisk and McClelland (1959) have provided data for 
three locations of similar composition, but of different 
degrees of consolidation (Table 7.1). 

To apply the nonlinear theory to the above cases the 
following information must be known. 

(a) The material input values, such as the initial 

Youd) fat oy, ey and the corresponding initial 
effective stress Og’: 

(b) The mass rate at which deposition over the years 

took place so as to yield the present thickness. 
{c) The compressibility and permeability characteris - 

fics chery wiht taal esa at eo and the manner 

in which they vary. 

To determine the mass rate of deposition, the total 
weight and the original thickness of the sediment are cal- 


culated. From the data of McClelland (1967) the void ratio 
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and unit weight variation with depth are obtained. The 

total weight of the sediment is calculated. The average void 
ratio of the sediment is determined. It may be noted that 
the thickness of the sediment is proportional to the quan- 


tity (l+e). Therefore 


ae ; (1+e5) 
Original thickness = Present thickness (mmr 
The value thus obtained gives a rough approximation of the 
actual rate at which sedimentation took place. The mass 
rate of deposition is then obtained by dividing the total 
weight with the original thickness. 

Thesinpwt intial word ratio @g may be termed as that 
void ratio of the sediment a few centimeters below the mud- 
line. Also the value of @g may be varied witn the water con- 
tent, w% (e, = W Go» where G. is the specific gravity of the 
Sioulei disi)s, 

To arrive at the input value for the effective stress, 
Ty's a relation between e and o' must be known. The sediment 
at Eugene Island Block 188 is fully consolidated and behaves 


as a normally consolidated soil (Figure 7.4). Thus the soil 


behaviour is expressed as 
elise 1.048 -90:.79 1og19 ) 
This relation is made use of for the sediment at the other 


two locations 


To determine the variation of permeability, a relation 
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between m/1+te (m, the coefficient of volume compressibility) 
and ag' is drawn (Figure 7.5). If the assumption is made that 
the value of the coefficient of consolidation, Cy is a con- 


Spant. the following relation is obtained. 
nies Ee) aoa 


An estimate of the value of Ko corresponding to Tg’ 
Can then be made. The value of Ky was varied together with 
imervalue for q-(in the equation 7.18) for the underconsoli- 
dated sediment at Grand Isle Block 23 to fit the observed 
foromon void ratio with depth (Figure 7.8). This 1s carried 
out since the rate of deposition for Eugene Island Block 188 
Sediment is so slow that a variation of c, and/or q does not 
have any influence on the degree of consolidation. Figure 
7.6 yields the variation of te (i.e., the variation of kg) 


with q. <A decrease in the value of q is compensated by an 


increase in a ia-tierre lation 
0 
i ]+ 
(lus \2 ie ) q 


V Vo Ite) 9 


so as to yield the same void ratio-depth relationship as the 
observed (Figure 7.8). It has been observed from Figure 7.2 


that the values of q for the various soils tested by Normand 


(1964) are: 
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Weald Clay GR=t=r0 767 
Avonmouth Clay q = -20..60 
Kaolinite Ge -a-80263% 


Also from Figure 7.5 for the fully consolidated sediment 

at Eugene Island Block 188 the value of q obtained is -0.814 
(when Cy is assumed constant). Keeping these values in view, 
for purposes of this study, the following set of values are 


chosen: 


qrewi-foezdmee & LOMnT oeyRN 


(poe hates Poe 50 FT@/YR. 


A value of 1 x 107° em/Sec stom on x ere cm/sec for kp is 
obtained, which is quite reasonable for the type of sediment 
dealt with. The properties of the sediments used in the 


analysis are: 


era leo4 = 0.75 10946 On 


Figures* /.7,9788, and /699 llistrace the relationship 
between the void ratio and depth for the three locations. 
The observed data is represented by points while the calculated 
results are shown by full lines. The calculated values agree 
very well with the observed data. Figure 7.7 yields the re- 
lationship between void ratio and depth for two input void 
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ence in the curves lies only in upper few feet of the sedi- 
ment. A similar trend was observed at the other two locations 
(this is not shown in Figures 7.8 and 7.9). 

Begures /.10, / all and 7.02 show the manner an which 
permeability decreases with depth at the end of deposition. 
ire Variationvis typical and it “decreases “from a value of 
ee Oeror K/K to a value of the order of 0.04. However, it 
is noted that the variation of permeability with depth for 
the South Pass Block 20 sediment (Figure 7.12) is very small 
in the upper few feet and then gradually the decrease grows 
rapidly. This is because the sediment is very underconsoli- 
dated and the upper few feet are almost at the input void 
ratio ey (and effective stress Og'). As the sediment is 
being built/up the effective stress increases and hence <a 
decrease in the permeability occurs. 

migures 72132 72 14 ynande/lom Llustravestnetvaniation 
of the coefficient of consolidation with depth for the input 


values of 
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The variation (increase or décrease) of Cylcy with depth is 
mainly dependent on the value of q. The smaller the value 

of q (smaller than -1.0) the value of c, decreases with depth. 
However, a value of q greater than -1,0- makes» ithe, value of 


c. increase with depth although the increase is very small. 
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This demonstrates the manner in which the magnitude of q 


atfects both k and Cy with depth. 


Wee 0: UES CUS STON 


Waa OSes Udy, Cie sSta vel of -consolidation wat several 
Hocations on the continental shelf off Louisiana is studied. 
The one-dimensional consolidation theory adopted for analysis 
is the most general accounting for a mass input at a pre- 
scribed initial void ratio. The results calculated are in 
good agreement with the observed field values. The rate of 
deposition.and the soil characteristics such as the compres- 
sion index, aieian void ratio, and coefficient of permeabi- 
lity control the pore pressure dissipation (and hence the 
effective stress). 

A close study of the equation 7.20 reveals that a 
value for the exponential q has a significant influence on 
tneaultimate results...A value of zero. for q reduces: the 


equation, 7.20 to 


The above equation does not embody any more the term con- 

taining G.. In other words, if by any chance the ratio k/It+te 
S 

remains constant (i.e., q = 0) which is most unlikely, during 


consolidation the specific gravity G. (or the self weight) of 
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the consolidating material does not contribute anything to 
the progress of consolidation. 
Also from equation 7.20 the expression it Veeder 


contains both the hydraulic gradient ae and Gatisucire 


Specific gravity of solids. If either one of these terms 


is zero then the expression au my CGS 1) vanishes. The hy- 
draulic gradient (su) is zero at an impervious boundary where 
there is no exchange of water across it. At such a boundary, 
mueeseems , the influence of G. tLsanovere teat choughmated| 
other sections G. has its influence. By the same token the 
influence of G. is the maximum where the hydraulic gradient 
jmsythe highest. The hydraulic gradient is the highest at 

the top of the sediment and hence the influence of G. is 
maximum in the upper layers of the sediment. 

To assess the influence of G. analytical sorrttons 
were obtained for values of G. Secepednd G. = 1.0 (buoyant 
weight). It has been observed that the computed results 
of void ratio (with depth) and the degree of consolidation 
do not differ substantially. It may be concluded that G. 
(i.e., self weight of the sediment) is not too important in 
Quantitative results. 

In summary the following may be expressed. 

The deltaic sediment samples (very under consolidated 
to fully consolidated) reported by McClelland (1967) were 
examined in the light of the theory presented in section 7.4, 


The rate of deposition is believed to be the main factor 


most closely associated with the consolidation characteristics 
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observed. The rates of sedimentation and related properties 
are available in sedimentology literature. In Ge OWOiss On 
high rate of sedimentation, there will be a lag between the 
accumulation of the material and the consolidation associated 
with it. This gives rise to an excess pore pressure and the 
under consolidated material is prone to slumping. The pre- 
diction of excess pore pressure can be made to a high degree 
Of accuracy by the analysis presented in section 7.4. The 
State of consolidation of the clay mass is very important 

in connection with foundation studies fOr Opnonore structures, 
for the strength of these deposits is a function of pre- 
consolidation pressure. 

Other parameters that affect the build up pore pressures 
during sedimentation are the material properties such as 

(a) the consolidation characteristics of the sediment; 

(b) the variation of compressibility and permeability; 
and (c) the environment in which deposition takes place. 

It has been demonstrated by Skempton (1970) that the 
behavior of (fully consolidated) deltaic sediments is mostly 
normally consolidated. And hence the value of p in equation 
7.16 may be taken as equal to -1.0 for all practical con- 
Siderations. 

The sedimented void ratio e) (and hence the corres- 
ponding effective stress) depends on the type of environment 
in which the sedimentation takes place. Depending on the 
sedimented material properties such as eg and the rate of 


sedimentation the behaviour of the sediment at any stage of 


Deu eet iy otal i + i 6 e ne gay Br bet i) 


| ae a 
is aig ' eo i haat | aT uotod an 
i ba pe r Bis Pe ot 4 a or - 1 i et f r om) ee - : a a " : 


a.” Wert. sie ro tee Seana wee, Dhb er ahi ‘id ie 


in. elt r BS nawals ee | scale 
hr a i nie ee hi 40 atts) oe” 


a B a ag si: coe i +e . | 
ath’ | bvgeteigees aie rut ane outa: 4d (ee 


if 
hs we nor shit tora 

Le a 
i elt oe i2Nre ‘ny - ad) bake nova 0% fey die 595 si 
arr" | oe dint 
‘hVi-z vieont a9 ni gaara 


a 


eintoatn skcg oy bf loud ahi do0hpe tans eae smote yi ys : 


aweReng, nor 


oh ies ® wong forestey ore iotaeseaqnhee 


Aswah eee wit Mo aie! ete np 4 anion 


Varin doe yar Te dpersigneg) Me not yelnay, sig, ia | 
| be tt ae srouna tno 6 
| foil CHOY): hols L bate wajeiupiey nave @ pote 
¥i de i uy e9R i fuse ote ifob (ha: +a Sioa sg turiere 


ri 


Hobie No 70 Chiay. off Mone bal op Sta 


= 
¥ > 
B fd 
en 
‘ 
4 


A. 
sa 


wo) (hereaehg 118 AGT O. 1- @2 Veoge a oon oy. een a 


wrist Raed Dee) Yd otdad hhen - bs mat bee i 


ay 
~ 


he 
¥ ; ns 
saa fi 


‘avune Lived Dy evs he io eeneyeh (azhade ovtsvarte 
Oe oy. po Me gee a riedsd hotdedanntiing i ‘om 
Aig wher nt bee yet on "doe std 7EQONG fatnedan 


to no eta qa 3 potiet tne at) ye. welsh 


23.8 


accumulation can be predicted provided the value of sedi- 

mented permeability Ko (i.e., the input coefficient of con- 

SOolidation Cy ) and its behaviour during sedimentation are 
0 


known. The magnitude of the values of k, and q of the sedi- 


0 
ments off Louisiana are precisely not known. It has been 
smown that q has a-significant influence on the variation 
of the coefficient of consolidation and the consolidation 


Giaracteristics of the sediment. 
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FIG, 7.8 RELATIONSHIP BETWEEN VOID RATIO AND DEPTH, 
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FIG. 7.12 VARIATION OF PERMEABILITY WITH DEPTH, 
SOUTH PASS BLOCK, 20. 
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WITH DEPTH, EUGENE ISLAND BLOCK, 188. 
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CHAPTER VIII 


EROSION AND SWELLING IN VALLEY FORMATION 


Gal - GENERAL 


This chapter deals with the pore water pressure-time 
relationship for a fully saturated clay layer subjected to 
a uniformly-distributed eroding load. The governing equa- 
tion is developed on the basis of the Terzaghi's classical 
one-dimensional consolidation theory and it is weakly non- 
linear in character. With this nonlinearity and the boundary 
Conditions encountered, it is difficult to obtain a closed 
form solution. A numerical solution based on finite differ- 
ence approach is presented. The study is of practical value 
in the assessment of pore pressure equalization of an eroded 
soul masse. Further, the study aimsirat corre latingigeotogica | 


and physical parameters involved with the theoretical solution. 


8.2 FORMULATION 


The problem considered is that of a semi-infinite mass 
Gest) eroded ata: presornibed arbitraryimate at the top. The 
properties of a real soil are time and space dependent. For 
purposes of this study, however, the material properties are 
assumed constant (on the lines of Terzaghi's classical theory), 
i.e., the physical properties of the soil such as the coeffi- 


cient of swelling (consolidation), the coefficient of perme- 
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abial 1 bysnetCG. vane constant. 
The problem may be considered to Fepres cio canny n0,f 
the following categories: | 
(a) the eroded material removed instantaneously; and 
(b) the usual case of finite rate of removal of 
eroded material. 
Thus the problem is one of a moving boundary condition. 
Generally, two types of moving boundary problems are encoun- 
tered: 
(i) motion of the boundary is due to the rigid body 
motions of the entire body (Carslaw and Jaeger, 
1959); 
(ii) motion of the boundary is due to-local conditions 
near the boundary. 
Further, in tyye (ii) two cases can be distinguished: 
(a) problems in which the motion of the boundary is 
prescribed; and 
(b) problem.in which the boundary must be determined 
as a part of the solution of the erosion problem. 
The problem of erosion is developed along the lines 
given by Gibson (1958) for sedimentation. In the case of 
sedimentation fresh layers of soil are added at the top 
whereby the thickness of the sediment increases with time. 
In the erosion study, the layers are removed gradually at 
the top and hence the thickness of sediment decreases with 
time. Thus erosion may be looked upon as similar to sedimen- 


tation except that instead of soil layers being added at the 
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top, the soil layers are removed Gmadually, dice kconsitant 
Batev?’ Thus §the quantity Ac in Gibson's (1958) equation for 
sedimentation (equation 7.13b) may be changed to -Ao; then 


equation 8.1 will describe a process of erosion. 


d 
= i (Ac) 8.1 


where Cy denotes the coefficient of swelling of the material 
u denotes the deficient pore pressure 


Ao denotes the weight of the soil removed by erosion. 


Equation 8.1 governs the dissipation of deficient 
pore water pressures in one dimension when erosion takes 
place ata prescribed rate of d(Ac)/dt. The expression Ao 
is equal to y dh/dt where y is the bulk unit weight of the 


soil and dh/dt is the rate of removal of the soil with time. 


BOUNDARY CONDITIONS 


A unique solution to equation 8.1 depends primarily 
on the boundary conditions. The top of the sediment is al- 
ways open to the atmosphere and is) free to drain; thespore 
pressures will dissipate almost instantaneously. Pius sae 
the top the excess (deficient) pore pressure will always be 
zero. But then the top boundary is 'moving' with time. The 


top boundary is to be located at any instant OTreL ime. 
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pt@the mins banteoretimente=00 4{ite.. anitialdyjecne 
top of the sediment is designated by z = 0. The rate of 
erosion dh/dt, is equal to m; the amount of sediment eroded 
Cleetsaz =-int after a time t. = hus thee top boundary at time 
te t exists at ‘a distance z = mt from the original refer- 
euce plane. It .issat this boundary the deficient pore pres- 
Sure is zero. Therefore the top boundary condition may now 


be defined as 
ULGZy tas0 sam Zhscm ty) 822710 8.2a 


ire Dotcom poundary COnd1.cion, condition mus talso 

Derspecitted. It-is tacitly assumed that the problem in- 
volves a semi-infinite homogeneous soil medium. It may be 
said with no loss of genrality that at a great depth (i.e., 
Zao). tere Will be no change an the deficient pore pres. 
sure (equal to the weight of soil removed) at any given time. 
Or, in other words, at infinite depth the soil will in effect 
act as an impermeable layer since across that layer no hy - 


draulic gradient will be set up. Thus the bottom boundary 


may be defined as 


cu = ( Fhe og ody, inaee 10) 8.2b 
cs ez 


The soil mass as a whole may be assumed to be in 


equilibrium with its surroundings before the process of 
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erosion sets in. As such it may be said that no deficient 
pore pressure exists throughout the soil mass at the instant 
When erosion “is about to set lin’ The dnitial condition, 


CAMS, 15 


“ez, ty) = 0 Way arucon t = 0 sae 


Making use of the boundary and initial conditions 
given by equations 8.2 a unique solution to equation 8.1 may 
be obtained. A closed form solution of the equation 8.1 for 
Ene nonlinear boundary conditions specified is difficult to 
obtain. Hence a numerical method is sought. 

The procedure for the numerical solution of equation 
8.1 involves replacing the equation by a finite difference 
fori. “Itdt- 1s accomplished by central differences using 
Crank-Nicolson's wmplicit method. The resulting set of 
Simultaneous linear equations is solved by the Gaussian 
elimination procedure. 

However, the nature of the bottom boundary condition 
specified may influence the accuracy of the results. Hence, 
the artificial boundary specified by z + » (equation 8.2b) 
was carefully studied, both with respect to its location and 
with respect to its assigned boundary condition, to determine 
its effect on the computed results. The conclusions were: 

Gus thatsa distance (Lon approximately 2000 feet is 


adequate in order to neglect the influence of 


the boundary; and 
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(b) that consideration of this boundary as impermeable 
ou 


yz = 0), or that at this boundary the pore pres- 
sure u(z,t) is always equal to the weight of the 
Soil eroded at that instant had no Son hiicant 
errectsshences tne solution was obtained for L = 
2000nTeset. ~Thetsollutionswasefound to be Stable, 
According to this method of solution, the behavior 
of the soil mass is determined at a finite number of discrete 
nodal points; if values of the parameter u are desired at 
points intermediate to the nodal points, linear interpola- 


tion may be employed. 


Oo DISCUSSION OF NUMERECAL (RESULTS 


Typical results are presented for the erosion of a 
low-swelling clay sediment under constant rate of erosion. 
The clay deposit may be considered as recently glaciated 
(quaternary age) and as being eroded by fluvial action. The 
Sediment issassumed tovbe fully*saturated* and» on the top of 
the eroded layer, it is assumed, water is always available. 
By this means, the deficient pore pressures developed are 
Satisfied at the boundary. Further it is assumed that uni- 
form pressure distribution takes place throughout the depth 
when a soil layer is removed by erosion. 

Eigquresecs|at0e8.deshowlgthe deficientu pore pressure 
isochrones at the end of erosion and also at the end of the 


duration of standstill. The depth of erosion is atbiteertly 
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chosen at 250 feet. The plots have been obtained for various 
values of the coefficient of consolidation (WO ae Phe seo— 
efficient of swelling). The values for the coefficient of 
Swelling C. Us Gd Hae Ol, Gh 50n) 1.0'0. “and hO0n0ahect Square 
wemnyear. hAlsoevarying periods. of oan and subsequent 
standstill are tried. The combinations are: 

(a) 5000 years of erosion and 15000 years of subse- 

quent swelling; 

(b) 10000 years of erosion and 10000 years of sub- 

sequent swelling; and 

(c) 15000 years of erosion and 9000 years of subsequent 

swelling. 
In each case the value of the coefficient of swelling is 
varied through the specified ranges. 

Since it was assumed that the sediment is fully satu- 
rated and also that the groundwater level is flush witn the 
top of the (eroded) sediment, prior to erosion there exists 
a hydrostatic pressure of water which would be recovered after 
equilization. Combining the deficient pore pressure with 
the hydrostatic pressure at each level the resulting pore 
water pressure isochrones, that are existing, are obtained. 
Biguvesecal tosS.4 ihlustrate thatyat thegend sok the period 
of erosion, substantial deficient pore pressures remain at 
the top within few feet although the sediment is: fuddyusatu- 
rated; thus there exist very high negative hydraulic gradients 
at the top. These high deficient pore pressures Shah) exist 


after a large number of years of standstill, when neither 
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deposition nor erosion takes place. 

These figures show that the dissipation of deficient 
pore pressure depends, largely upon the coefficient of swell- 
iatg!, C.. The higher the value of C. the quicker the dissipa- 
tion and the greater is the depth to which pore pressures 
are equated. Thus the behaviour of a soil sediment depends 
mainly on the permeability and compressibility characteristics 
(on which C, depends). The lower the permeability, the longer 
une =trmei needed for the dissipation of! deficient pore pres- 
Suir 2Sy 

For a given value of Cos the duration of erosion and 
Subsequent period of standstill have a profound influence on 
the magnitude of the developed deficient pore pressure. The 
smaller the duration of erosion (for a given erodable depth), 
the greater the effect on the deficient pore pressures and 
uve hydraulic gradient at upper Payers is’ very Targerwe The 
longer the duration of erosion (for the same eroded depth), 
comparatively the hydraulic gradients are smaller. To 
illustrate the above statements, ‘the total’ water ‘pressure 
at anylipoint (as itexisits at the end of the duration of 
standstill) is expressed as a percentage of the nydrostatic 
pressure at that point that would have existed had there been 
no swelling. The calculated per cent of thydros tat c pressure 
is plotted against depth for a given value of C.. Figures 
8.5 to 8.8 illustrate the effect of rates of erosion and 
subsequent standstill on the pore water pressure for a speci- 


fied soil mass (whose coefficient of swelling is known). The 
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values expressed, e.g., as 5000 + 10000 years on Figures 
8.5 to 8.8 denote 5000 years (first term) of erosion to 
Form a 250 feet valley and 10000 years (second term) of 
subsequent standstill with no erosion or deposition. 

Figures 8.5 and 8.6 suggest that whatever be the rates 
of erosion and subsequent swelling for soil masses with a 
value of C. OF etivetordertom 0.1 or 1.0 foot savare per year, 
the per cent hydrostatic pressure does not alter with depth. 
This means to say that in low swelling soils the dissipation 
oie pOnecmprossiures isevery slow .2r Figure 8l8 sitgqgesits® that for 
Soils with a high value of C. (e.g., 100 feet square per 
year), the rate of erosion and swelling have almost the same 
ehrecu. ln other words, tne hign swell ing soi lss adjust 
themselves to the surroundings very quickly. A value of 
10.0 feet square per year for C. (Figure 8.7) shows a marked 
immluence of the various*rates of erosion and subsequent 
swelling. “AIL the curves (Figures 8.5 to 8.8) become paral- 
liel@ebowtive: denith vaxis: ‘at; largeivallues: of deptn. 

| Essentially the same remarks (as for the rate of 

erosion) may be made for the duration of subsequent stand- 
stil 1 ethic length of the duration of swelling has a signi- 
ficant influence on the deficient pore pressures. The longer 
is this duration the more is the degree of dissipation. 

It is of interest to investigate if cases exist in 
nature that confirm to the type of behavior described earlier 
in this section. Large areas of Western Canada and the 


North Central States (USA) are underlain by Bearpaw Shale. 
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This material geologically is classed as bedrock. 
In the undisturbed state Bearpaw Shale is dense, 
homogeneous and impervious (Peterson, 1954). Some of the 


pertinent material properties are: 


Average Degree of Saturation 96.4% 
Coefficient of Permeability cm/sec HG sto on” 
Average Water Content MASSE 


Tests conducted at this University yield an average value 

of 0.10 feet square per year for the coefficient of swelling, 
Cy. Same order of magnitude for Cy is quoted in the published 
literature. 

The analysis presented in the previous section is rele- 
vant to the valley formation in Bearpaw Shale. South Saska- 
tchewan River drains most of the Prairie provinces and has 
cut deep valleys in the Bearpaw Shale formation. In the 
immediate vicinity of Gardiner Dam the valley is approximately 
-200 feet deep. Thus the river flows through a deep bed of 
homogeneous, low swelling, (almost) saturated sediment and 
has cut a valley 200 feet deep. Figures 8.1 and 8.5 are 
ideal for comparison for this case and the figures indicate 
that a substantial deficit pore pressures may exist. 

Although it has been assumed that large depths of 
homogeneous material are involved, the analysis can be applied 
to homogeneous materials of limited thickness underlain by 
relatively impermeable material such as most rocks. However, 


any inhomogenity in form of stratification in the material 


will yield substantially different result tS. 
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= Pore Pressure (units) + 
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FIG. 8.2 PORE PRESSURE ISOCHRONES 
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= Pore Pressure (units) + 
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FIG. 8.5 PERCENT OF HYDROSTATIC PRESSURE vs DEPTH 
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FIG. 8.6 PERCENT OF HYDROSTATIC PRESSURE vs DEPTH 
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FIG, 8.7 PERCENT OF HYDROSTATIC PRESSURE vs DEPTH 
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COMPUTER PROGRAM AND USAGE 
FOR TWO-DIMENSIONAL IMPEDED DRAINAGE 
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APPENDIX A 


ORGANIZATION OF COMPUTER PROGRAM 


The computer program described in this Appendix is 
based on the theory presented in Chapter V of this Study. 
The program is in Fortran IV language and may be used directly 
on computers of the type IBM 360/67. 

The program furnishes the excess pore pressure, and 
the average degree of consolidation of a rectangular wedge 
at any time. The program is applicable for a rectangular 
section of neight L and width 2H provided with side drains 
of thickness d on either side. The program is based on the’ 
ADI method of analysis presented in Chapter V. The program 
is presented for uniform distribution of initial excess pore 
pressure which may exist at the start of the program and 
which will dissipate with time. However, the algorithm can 
easily be modified to incorporate a continuous variation 
across the entire height or width. The variation may be 
different from one direction to the perpendicular direction. 
The variation may be defined by a polynomial interpolation 
of input data. 

The equation governing dissipation of excess pore 
pressure in two dimensions is given by 5.2. The extent of 
retardation to drainage is measured by a parameter called 
the Impedance Factor, A. The treatment of the problem is 


in a dimensionless form. 


The program reads all the input data regarding the 
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rectangular wedge and sets up the initial and boundary condi- 
mons. Also it establishes the impedance factor governing 

the retarded consolidation. Depending on the governing 
material and operating parameters, the program calculates 

the excess pore pressure at each nodal point at presicribed 
time intervals. ye average degree of consolidation both 

by two-dimensional and one-dimensional (along a specified 
section) are calculated using the pore pressure values 
obtained earlier. The Simpson's rule is made use of in 
obtaining the degree of consolidation. The calculated results 


are printed out as per the output formats. 


COMPUTER PROGRAM USAGE 


INPUT DATA 


The first step in the analysis of retarded consoli- 
dation for a rectangular wedge with inefficient side drains 
is to determine 

(a) the height to half width ratio ‘Beta’ of the 

wedge; 

(b) the existing boundary and initial conditions; 

(c) the extent of impedance to drainage as measured 

by a numerical value of the impedance factor i; and 

(d) the space and time intervals to suit the computer 

memory requirements. 

The input data cards reads the space steps in both the 


@irections, the time (factor) steps and the ratio Beta. 
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OUTPUT INFORMATION 


The following information is computed and printed 
by the program at specified time intervals: 
(a} the pore pressure excess at each nodal point; and 
(be tne ssduare Leousof time factor, thestine factor, 
the degree of consolidation (by two-dimensional 
integration) for the entire wedge, and the degree 
of consolidation (by one-dimensional integration) 


along a chosen section. 
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APPENDIX B 


ORGANIZATION OF COMPUTER PROGRAM 


The computer program presented in this Appendix 
describes the construction sequence of JARI Dam (section 
6.4.5) and calculates the pore water pressures at several 
nodal points. The program is in Fortran IV language and 
may be used directly on computers of the type IBM 360/67. 

The algorithm developed typically predicts the influ- 
ences of dissipation on construction pore pressures. Of the 
Six case nistories studied, the program developed and written 
for Jari Dam is presented as an illustrative example. 

Tne equation governing the dissipation of construction 
pore pressures is given by 6.4. The method of analysis is 
Oyeine Die technique described in Chapter V7 ine pore prese 
sure generation term is DH (= y B*DX where DX is the thick- 
ness of the added soil layer) and this is equated to zero 
when there is a construction stoppage. 

The program furnishes the excess pore pressure at the 
nodal points of the idealized cross-section of the core of 
thewdart Dam. The first step in the analysis 1s to deter— 
mine the construction sequence and the rates of construction 
and work stoppage. The material properties of the core such 
as unit weight and the coefficient of consolidation Cy 
(alternatively the permeability avid compress ipl 1 hly -onaracter- 
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that the algorithm can handle any variation in the values 
of oT OL B at any Stage of construction or work Stoppage. 
The program is furnished for a uniform distribution of the 
load added at the top (due to increase in layer thickness) 
on all the nodal points irrespective of their position and 
Pepin ihis Ws justitied to an extent. because proper dis- 
tribution of load along the depth and along the lateral 
eirection fOr a moving boundary problem is difficult to assess . 

the entire program consists of a main and fourteen 
Subroutines. The main program reads in the input data and 
Sercmup acne: initial andi boundary conditions. “lt establishes 
the physical constants such as core dimensions and the 
material properties. The main program also regulates the 
Dranching to the Subroutines by setting up appropriate 
conditional statements. 

Subroutine KOP generates the boundary conditions for 
the first construction stoppage of three months. Subroutines 
KOPPUL and KOPPU set up the boundary conditions for the 
second and third construction stoppages of one and half 
months each. Subroutine NM] computes the boundary condi- 
tions for the period after construction for the entire core 
section. 

The subsequent subroutines calculate the excess pore 
pressures generated (or dissipated) at the current nodal 
points by the ADI technique. The computed quantities in all 
the subroutines are returned to the main program. Tne main 


program prints the desired information according to the out- 
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COMPUTER PROGRAM USAGE 
INPUT DATA 


The first step in the analysis of the influence of 
dissipation on construction pore pressures is to decide on 
the cross-section of the core in the case of a dam to be 
built or to note the dimensions of an existing dam under 
examination. The core section is idealized on the basis 
described in Chapter VI to get a rectangular mesh amenable 
to finite difference approach. The boundary conditions at 
the bottom and sides are also noted. 

The second step is to note the rate of construction 
together with the durations of construction stoppage, if any. 
When necessary the actual rate of construction is modified 
to enable easy computation on the lines suggested in section 
Gr. BD. 

The most important parameters in calculating excess 
pore pressures are the coefficient of consolidation Cy and 
the pore pressure parameter B. These values could, however, 
be varied throughout the construction operation. 

The next step is to discretize the core of the dam 
into a rectangular grid and decide on the time interval so 


Eiatethe space. and time steps suit the computer memory re- 
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The input data card reads the space steps in both 
directions and the time steps. 


The READ variables are: 
WT sr 


The FORMAT is, 


OUTPUT INFORMATION 


The program computes and prints out the following 


information: 
(ae Brenprint of input data; 
(b) pore pressures at nodal points for specified 
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GO TO 674 
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CCATINUE 

GG Th 674 
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OO) GIL tn = WIN Se 
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SULLA) a TRACE SUR IE INSURE LLIN SINT) AEC WG a2 Go eSROESUNUE OIL IN INI IBA 
CUR/SA®*2)*R(LA-1L)I/01L.4+02. FIRS SA¥X%2) D-(RISA**®2) 4 A(LA-1) ) 
CONTINUE 

Jy = IGS 
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CUR/SAX= #2) XB ULA-1))/0 1.412. 3(R/SA%*2) J-(R/SAF*2)KALLA-1)) 

CONTINUE 

UCL» JB Nel) = B(JB) 

Leet e Qk JA =. 2 

JE 2 dey = hs 

OO) Gy Wan = wy whe 

MB = JB-VM 

Otis Uitte) = BINB)FA(MB) «UCL ,MBtl Nel) 

CONTINUE 

CONT TINUE 

RETURN 

END 
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APPENDIX C 


COMPUTER PROGRAM AND USAGE FOR 
ONE-DIMENSIONAL FINITE STRAIN CONSOLIDATION 
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ORGANIZATION OF COMPUTER PROGRAM 


The computer program described in this Appendix is 
based on the theory presented in Chap teneVLh of sth hs siudy. 
The program is in Fortran IV language and may be used on 
computers of the type IBM 360/67. 

The program deals with the one-dimensional, non- 
linear consolidation behaviour of a sedimenting soil under- 
going finite strains. The program furnishes the excess pore 
pressure and the void ratio throughout the current depth at 
chosen intervals for any time step. The program also supplies 
the current height, the total added height, and the average 
current density at each time step. Finally the program cal- 
culates the average degree of consolidation of the sedimented 
thickness at the end of the duration of deposition. 

The program is for a moving boundary problem where 
the top boundary is in perpetual motion. The program is 
applicable for the case of a normally consolidating soil 
whose compressibility and consolidation characteristics are 
known. The program is presented for a uniform distribution 
of the added weight on the top (due to sedimenting soil) 
throughout the current depth of deposit. However, the 
algorithm can easily be modified to incorporate a continuous 


Variation of the added weight across the depth. 


The equation governing one-dimensional, finite strain 
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consolidation behaviour of a Sedimenting soil is given by 
7.20. The average degree of consolidation is defined as 

the ratio of the difference between the total added weight 
and the total excess pore pressure to the total added weight. 

The main program reads all the input data concerning 
the sedimenting soil and sets up the initial and boundary 
conditions. Also it establishes the space and time intervals 
and the initial effective stress of the very thin layer over 
which sedimentation takes place. 

Subroutine FS] calculates the relevant quantities 
depending on the material property values supplied by the 
main program when there is a steady growth of soil layers. 

{t also calculates the current void ratio for each of the 
nodal points at a specified time. The calculated quantities 
are then returned to the main program. 

Subroutine FS2 calculates the required dient when 
Cuemdenostt 1S ata standstill and there is mo Input th form 
of growth of layers. The dissipation is allowed to take 
place at a constant load. The computed quantities (same as 
in subroutine FS1) are returned to the main program for 
furchner calculation. 

The main program makes use of the quantities obtained 
from the subroutine(s) to calculate the excess pore pressure 
for each of the existing nodal points at prescribed time 
intervals. Further the program calculates the current height 
of the sediment (by trapezoidal rule), the average current 


density (of the sediment thickness), and the average degree 
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Of consolidation at the end of the deposition period by 
Simpson's rule. The computed results are printed as per 


the output format. 


COMPUTER PROGRAM USAGE 


—.__-. 


INPUT DATA 


ne TiyStestep in the analysis of one-dimensional. 
ees train consolidation of a sedimenting soil is to 
determine the rate of mass deposition and its duration. 
Next, the compressibility and permeability characteristics 
of the soil are to be determined. A choice is to be made 
as to the input void ratio and corresponding effective stress 
at which the soil is being deposited. 
The required input data cards are as follows in the 
order presented. 
(ay) First card. Mhe first inbutedata arc. cads inc 
specific gravity, the input void ratio and the 
input coefficient of consolidation. 


The READ variables are: 
GSO IE0;, They 
The FORMAT is, 
Sh L0ne 


(b) Second card: The next data card reads the numeri- 


cal values of C and D in the equation 
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See Cle cee.) Togi9 Ge 
and the value of Q in the equation 
k a(g!)9 
The READ variables are: 
Cee Dee 0 
The FORMAT is, 
Bie NOs} 


(c) Third card: The third input data card reads the 
time steps required for the entire deposition 
ana’ standstitl operation and” tnat fom on wivene 


duration of deposition. 
The READ variables are: 
Kal 
The FORMAT is, 
215 


The space step (DDZ), the time sitep (DDT), and the 
unit weignt of water (UWW) are built in the main program. 


The appropriate values are to be substituted for these 


quantities. 
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OUTPUT INFORMATION 


The following information is printed by the program 


at prescribed time intervals: 


(a) 


Che, initial (input) properties of the sediment- 
ing soil together with the time and space steps, 
the weight of the added layer thickness and the 
unit weight of water; 

the current height and total added thickness; 
the average current density; 

the excess pore pressure and void ratio at 
specified space intervals; and 

the average degree of consolidation of the soil 


deposit at the end of the deposition period. 
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(6, RAE LK A IE ole RKC Oke drake ak tere ts ak teak dete ok ok fc He Ae Ae Ke eke He Ye te Xe ake Xt He ke ek Ke 
C PROGRAM FOR THE ANALYSIS UF ONE DIMENSTONAL NONLINCAPR CONSOLIDA 
C -TION SFOIMENTATION FOK A PERICO & THEPEAFTER DESSIP. 
2 THE SOIL CONS IDFRFD 1S NORMALLY CONSOLIDATED 
C RR LO KOR A Kok RO a tr tok tle eek A atk tek IO ROI LK AO LAC te oleate 
c 
ON ALOU TIVO), UCIO2,102), EC 102,102), SIGMAILOL),XAC1G2), 
Pee Os ULCIOVY) RVCIOLL: S101), AN(IO1), PETLYOL), X10101), 
Oks Kl, Ls Ny GS, Q, UWW, D7, OT, 2, OH, SIGMO, EQ, C, D 
4 | . pee Cae kes 
(ORR Eg Rokk tate rota lak Ate fe teak tl tok ae tak fe x: 
C 
a PEAD AND PRINT MATERIAL PROPERTIES 
C 
os OK OR Rr a 3K Re Xe oe EK REX a ae Xe ke oe ae He KE ake 
: ‘ 
By FORMAT. (3710.3) 
E05 5166S 060, CV 
BEAD (54 o)-Cy 0, 0 
g 
REAN (5,10) Ky Kl 
10 FORMAT (215) 
C 


Perie ois) GS, FO, CV, Cy 05 0 
15 FORMAT (25X, "PROPERTIES OF SEDIMENTING SOIL'///, 25X, *GS=", 
Meee oN fF ERH), FB. 3/5) 25%), 1CVE (1 68a 3s) 25%s (Cals FS.S/s 
Mee eal, Fb.3/, 29%, 'Q=!, FR. 3/171), 
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THEREFORE OT = CV*DOT 
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N00 40 N=1,K 
NA = N 
IF (N.GT.K1) NA=K1 
NQ 30 L=1,NA 
peat Cee Oeste) Los GOs NA). GO, 11) 20 
GO. 10,30 
oo SUA, N) = Os 
E(L~N) = EO 
30) CONTINUE 
Oo) CONT LNUE 
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APPENDIX D 


COMPUTER PROGRAM AND USAGE FOR 
EROSION AND SWELLING IN VALLEY FORMATION 
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APPENDIX D 


ORGANIZATION OF THE COMPUTER PROGRAM 


The computer program presented in the Appendix is basec 
on the theory offered in Chapter VIII of this thesis. The 
program is in Fortran IV language and may be used directly 
On computers of the type IBM 360/67. 

The program determines the deficient pore pressures 
For a low swelling deposit eroded at a constant rate over 
a period of time and for an assumed standstill thereafter. 
Essentially the program is based on Gibson's (1958) theory 
for consolidation of sedimenting layers. The deposit is 
assumed to extend to infinity in its depth. The top of the 
deposit is eroded to a depth of 250 feel. Por ali uract)- 
cal purposes the infinite depth is achieved at about 2000 
feet from the uneroded surface. The algorithm developed is 
for a constant rate of erosion and for an uniform distribu- 
tion of load (removed) throughout the depth. However both 
these assumptions can be altered and the program modified 
for any sort of erosion schedule (which can be later appro- 
ximated into a suitable number of linear rates of erosion) 
and to include a continuous variation of the load removed 
across the depth. 

The governing equation for the dissipation of deficient 
pore pressures for an eroding sediment is given by 8.1. The 
analysis can be carried out for any value of the coefficient 


of consolidation and for several rates of erosion (or ec5U 
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feet depth of erosion) and subsequent swelling (same dura- 
tion for all rates of erosion). 

The program has all the input data built in and sets 
up the boundary and initial conditions. It establishes the 
rate of erosion and the space interval for the finite differ- 
ence grid. It then determines the time interval for the 
above space interval and the same time interval is used for 
the duration of subsequent swelling. Thus the total space 
steps, time steps for erosion and subsequent standstill are 
established. Depending on the governing material parameters, 
the program calculates the deficient pore pressure at each 
nodal point at prescribed time intervals. The calculated 
results are printed out as per the output formats. 

Ines programeaiso plots the deficient pore pressures 
adainst the depth: /This is carried out by the Calcomp 
plotter which operates on the principle that digital commands 
from a tape activate step motors to produce a plot. Further 
for various nodal points, the program calculates the hydro- 
static pressure of the eroded sediment and the existing pore 
water pressure at the end of erosion and at subsequent times. 
The existing pore water pressure is the sum of the deficient 
pore pressure and hydrostatic pressure. The prograi plots 


both these quantities against depth at prescribed time 


intervals. 
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QUTPUT INFORMATION 


The following information is computed and printed 


by the program: 


(1) 
a) 


(3) 


Peper ince fen put da ta: 
pore pressures (deficient at nodal points at 
Specified time intervals; 


plots the deficient pore pressure isochrones 


and existing water pressures at the end of 


erosion and subsequent time intervals and the 


hydrostatic pressure of the eroded sediment. 
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